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XX. BAND ZWEITES HEFT 1952 


Zur Theorie der turbulenten Strémung von Gasen 
stark veranderlicher Dichte *. 


Von W. Szablewski. 


1, Einleitung. Fundamental wichtige Phanomene in der Theorie der turbulenten Strémungen, 
wie die Ausbreitung heiBer Gasstrahlen in Luft oder der Warmeiibergang an umstrémten heiBen 
Kérpern, sind mit gréBeren Dichteunterschieden verbunden. Wir stellen zur Untersuchung dieser 
Phanomene in Verallgemeinerung der Betrachtungsweise von 0. Reynolds und L. Prandtl ein 
Gleichungssystem auf, das die turbulenten Mischungsvorgange in Feldern mit stark verander- 


licher Dichte beschreibt. 


2. Der turbulente Mischungsvorgang in statistischer Betrachtungsweise. Wir erinnern kurz 
an die klassischen Ansatze von O. Reynolds!. Bei Beschrankung auf Felder von konstanter oder 
annihernd konstanter Dichte betrachtet O. Reynolds Strémungen, die im zeitlichen Mittel sta- 
tionar sind. Fiir diese fiihrt er eine statistische Betrachtungsweise durch, indem er bildet 


u=utw, v=v+’, G=H1L. 
Uberstrichene GréBen bedeuten zeitliche Mittelwerte, gestrichene GréBen momentane Schwan- 
kungen. Ferner bedeuten x die Hauptstrémungsrichtung, y die Querrichtung, u, v die Geschwin- 
digkeitskomponenten in x- und y-Richtung, @ die Dichte, 7? die Temperatur, T die absolute Tem- 
peratur, c, die spezifische Warme bei konstantem Druck. Vom Standpunkt der im zeitlichen 
Mittel stationaren Strémung ergibt sich dann eine turbulente Schubspannung 


t=—o0u'v’ (la) 
bzw. eine turbulente Warmeleitung 
Q0=—c,00'v’. (1b) 
Wir gehen nun dazu iiber, die statistische Betrachtungsweise auf Felder mit stark verander- 
licher Dichte zu verallgemeinern, und betrachten hierzu einen ebenen Strémungsvorgang, der ge- 
kennzeichnet sei durch eine Hauptstrémungsrichtung mit turbulentem Austausch quer zu der- 
selben. Die Strémung soll im zeitlichen Mittel stationar sein. 
Wir kénnen den Strémungsvorgang, wenn wir molekulare Wirkungen (Diffusion, Reibung, 


Warmeleitung) vernachlassigen und der Einfachheit halber konstanten Druck annehmen®, durch 
folgende Kontinuititsgleichungen beschreiben: 


dou dou 
Masse: A ape (ee 
_ OQuu , OovU _ D) 
Impuls: - Fae re OF (2) 


Energie: — Tae — 
Zeitliche Mittelwerte sind wieder durch Uberstreichen gekennzeichnet. 


me Auszug aus der Géttinger Dissertation 1947: W. Szablewski, Die Ausbreitung eines Heifluftstrahls 
in bewegter Luft. 

1Q, Reynolds, Phil. Trans. Roy. Soc. 186 (1895). ate 

2 Bei Phanomenen der freien Turbulenz entspricht die Annahme konstanten Drucks den itblichen 
Voraussetzungen. Bei Grenzschichtproblemen haben wir, wenn ein Druckgradient vorhanden ist, auf der 


dp oe as 
rechten Seite der Impulsgleichung (2) das Druckglied — a auf der rechten Seite der Energiegleichung (2) 


tritt das Glied a auf, demgegeniiber die von den Schwankungen herriihrenden Glieder u as ae = 
x ; i 


als vernachlassigbar anzusehen sind. 
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Wir beschriinken uns dabei auf die Betrachtung der Warmeenergie als Reprasentanten der 
substantiellen Eigenschaften. Wiirden wir auch noch Beimengungen irgendwelcher Art der Kon- 
zentration ¢ zulassen, so hatten wir als weitere Gleichung 


deren Betrachtung ganz analog der Energiegleichung verlauft. 
Unser Ziel ist es, die stationaren Werte u, v, 0, cp T’ zu berechnen. Zu diesem Zweck stellen wir 
U, Uv, 0,¢, 0 durch ihre zeitlichen Mittelwerte und Schwankungen dar: 


u=itw, w=tty, e=o+e', |T=oT+(p7). 
Ks ist dann ou =pu ep ue. OU =pv+ 0'v’, 
und aus der Kontinuitatsgleichung der Masse wird 
Bey BS (ov) +507). (3a) 
Beziiglich der Impulsgleichung erhalten wir 
guu—ou- w+ (ou)'w —puu+orw ut (ou), 
Quv—pu-v + (guy =guv+o'w-v + (guy; 


das ergibt 
Cg eg OR VES Oe pul = (omy al no as a one (3b) 
Ox oy 0 . 9 oy : 
Entsprechend kommt fiir die Warmeenergie 
O0UCpT O0UCpT rs) rH ee TER aL OW Ca y ; i ee, 
oS = z.|-¢'u - oT — (Quy TY] + 5 l-e' TY -& (omy e']. (3c) 


Vom Standpunkt der im zeitlichen Mittel stationdren Strémung aus werden also durch den 
Mechanismus der turbulenten Vermischung turbulente Stréme der Masse, des Impulses und der 
Warme erzeugt. 

Wir betrachten nun die Gré®enordnung der durch die turbulenten Schwankungen induzierten 
Stréme. Bei der Kontinuitatsgleichung der Masse 

A(outew’) , Agvu+e'e) 4 
Seem sone, “= 
Ox oy 


handelt es sich um die GréSenordnung der Glieder o’u’lou und 9’v’/ov. 

Nach bisherigen Erfahrungswerten, die sich allerdings nur auf maBige Dichteunterschiede be- 
ziehen', iiberschreiten die relativen Schwankungen der Geschwindigkeit in der Hauptbewegungs- 
richtung u’/u kaum die GréBenordnung 0,1 und sind gegen 1 vernachlassigbar. Dasselbe gilt fiir 
die substantiellen Eigenschaften 0’/o, T’/T. Dagegen nimmt die relative Schwankung der Quer- 
geschwindigkeit v’/v erhebliche Werte an, da hier den kleinen Schwankungen v’ die ebenfalls 
kleine Querkomponente v gegeniibersteht. 


Wir vernachlassigen entsprechend o’u’ low und erhalten fiir die Kontinuitatsgleichung der 
Masse 

Te tet aa 

Be 4 Be 8 (ov) (4a) 


Bei der Impulsgleichung handelt es sich um die Glieder 


o’wjow, (ou) wjguu und (gu)'v’Jouv. 
Wir stellen die Schwankungen des Produktes (gu)’ durch die Schwankungen der Faktoren dar: 
(gu)’ = ou’ + wor. (5) 
Das folgt so: 
gu=ou + (ou), 


ferner 


| ou=(6+0')(utw)= p+ ow bie tow, 
1 H. Reichardt, Z. angew. Math. Mech. 24 (1944), S. 268. 


XX. Band 1952. Szablewski: Zur Theorie der turbulenten Strémung von Gasen usw. 69 


daher 


a eu (ou) = ou-+ ou’ + ug’+o'u’; 
hier ist pu—ou-0’u’, so daB wir schlieBlich erhalten 
(Qu)’ = ou’ +up’+ o'w—o'w’, 

0 


) 
car 


(ou)’ ~ou’+ uo’. 
Qu)'u we Te 
Fu MepyES erhalten wir dann — -+ ot 
u u? ou 
oe QU ry’ uv ‘yy 
fiir ( Ee ebenso ——- + oa 
Guv ut Ov 


Unter Vernachlassigung der Glieder kleinerer GréBenordnung reduziert sich also die Impuls- 


gleichung auf 


douu dovu OMe aR 
as aE = a (Qu) v ) ; (4b) 
Entsprechend ergibt sich fiir die Warmeenergie 
doucpT . dovepT a paca 
ae eee ay = By (= (oc,T) v Me (4c) 


Wir haben damit fiir den turbulenten Mischungsvorgang in Feldern mit stark veranderlicher 
Dichte erhalten 


M=~— 0’v’ fiir den turbulenten Massenstrom. (6a) 
T =—(ou)/v’ mit (ou) ~ou’ + uo’ fiir die turbulente Schubspannung, (6b) 
Q=—(0¢,T) v’ mit (o¢,T) ~o(c,T) + ¢,To’ fir die turbulente Warmeleitung. —_ (6c) 


Dabei ist zu beachten, daB die hier abgeleiteten turbulenten AustauschgréBen sich auf ,,Kon- 
tinuitatsgleichungen‘‘ beziehen, vgl. (4a) bis (4c). 

Bei auftretenden Konzentrationsunterschieden einer Beimengung tritt noch D=— (@c)’v’ (6d) 
fiir die turbulente Diffusion hinzu. Diesen turbulenten AustauschgréBen gegeniiber sind die 
molekularen, die wir vernachlassigt haben, meist bedeutend kleiner!. 

3. Die Prandtlsche Hypothese des Mischungsweges. Fiir Felder mit annahernd konstanter 


Dichte 9 hatte O. Reynolds auf statistischer Grundlage 


T=— 0u'v’, J —_ 000 


erhalten. Daneben stand langere Zeit unvermittelt der auf phanomenologischer Grundlage be- 
ruhende Ansatz von Boussinesq? 
ou 
SA a 
Hier ist ¢ die sog. scheinbare kinematische Zahigkeit. Das Produkt g¢ wird oft zu einer GriBe 
A,, dem sog. Austausch, zusammengefaBt : 


ee 


T dy . 
A, hat die Dimension einer Zahigkeit und ist meist um Zehnerpotenzen gréfer als die laminare 
Zahigkeit yu. 

Entsprechend gilt fiir die turbulente Warmeleitung 


Oo 
eI 


wo ¢, Ao die sog. turbulente Warmeleitzahl ist. 

Als ein Hauptproblem der Theorie der Turbulenz kann man die Aufgabe bezeichnen, die 
SchwankungsgréBen durch die zeitlichen Mittelwerte darzustellen. Ein wichtiger Schritt in 
dieser Richtung erfolgte durch die Prandtlsche Hypothese des Mischungsweges*, die gleich- 


zeitig eine Vermittlung zwischen dem statistischen und phinomenologischen Ansatz schuf. Nach 


1 H. Schlichting, Grenzschichttheorie, S. 338. Karlsruhe 1951. 
2 Boussinesq, Théorie de l’écoulement tourbillant, Paris 1897. 
3 L. Prandtl, Z. angew. Math. Mech. 5 (1925), S. 136. 
6* 
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dieser Hypothese entstehen die Schwankungen aus der Bewegung von Turbulenzballen, die fiir 
die Dauer ihrer Eigenexistenz ihre Geschwindigkeit und Temperatur von einer Schicht zur an- 
deren transportieren. In dieser Auffassung ergeben sich die mittleren Schwankungen der Ge- 
schwindigkeit in der Hauptstrémungsrichtung bzw. der Temperatur als proportional dem Gra- 
dienten der Geschwindigkeits- bzw. Temperaturverteilung : 


aor du al a8 | 
77 |% ey py eel if 
wi=hl|. | | 5 (7) 
Fir |v’| wird der plausible Ansatz gemacht 
|v'|~|u’|. (8) 
Mit den Korrelationen oe ae 
u’y’=k|u’| - |v’ (9) 
Ha = KF 
erhalt man dann unter Beachtung des richtigen Vorzeichens 
t=0F Ee fiir die turbulente Schubspannung 1, (10a) 
Oc or ae a6 fiir die turbulente Warmeleitung. (10b) 
TOCA 
Die scheinbare kinematische Zahigkeit wird demnach 
Ou 
ey fp pee 
emt lay 
und der Austausch 
du 
Seo 
AeA a 


In neuerer Zeit hat Prandtl? fiir Aufgaben aus dem Gebiet der freien Turbulenz vorgeschlagen, 
fiir den Austausch, der nach der obigen Formulierung itber dem Querschnitt variiert, einen Mittel- 
wert zu bilden. In der obigen Formulierung hat der Austausch differentiellen Charakter, was 
an Nullstellen der Ableitung dem Begriff des Mischungsweges nicht aquivalent ist und in ge- 
wissen Fallen denn auch den Tatsachen nicht gerecht wird (fiir die Ausbreitung von Luftstrahlen 
z. B. wird in Strahlmitte mit du/dy auch der Austausch gleich Null, was rechnerisch zu un- 
natiirlich zugespitzten Profilen fiihrt). Es wird infolgedessen der titber dem Querschnitt kon- 
stante Austausch in Vorschlag gebracht, der sich mit dem Mittelwert 

€(x) =x, 6, (x) | Umax — Umin| (11) 
ergibt, wobei b,(x) die Breite der Vermischungszone und x, eine empirische Konstante ist. Es ist 
ersichtlich, da durch diesen neuen Ansatz auch die Rechnungen erheblich vereinfacht werden. 

Wenden wir uns nun dem turbulenten Austausch in Feldern mit stark verinderlicher Dichte 
zu, so machen wir uns fiir den Mechanismus der turbulenten Vermischung die Prandtlsche Vor- 
stellung des Transportes durch Turbulenzballen zu eigen. Wir tragen jedoch dem Umstand Rech- 
nung, daf die experimentellen Beobachtungen (in Feldern mit annihernd konstanter Dichte) ein 
ungleiches Ubertragungsverhaltnis fiir den Impuls und die substantiellenEigenschaften anzeigen®. 
Setzen wir nach Prandtl 
| du 


eee 
|u’| 1 ay 


(12) 


1 Es moge bemerkt werden, dafs man den Prandtlschen Schubspannungsansatz auch in einfacher Weise 
aus der Hypothese erhalt, daB die turbulente Schubspannung bestimmt wird 
durchdenlokalenVerlauf dermittleren Geschwindigkeit u(y). Der lokale Ver- 


lauf von u(y) wird bei hinreichend kleiner Umgebung analytisch festgelegt durch = (== 0 vorausgesetzt). 
ask ly ; 

Der exakte Reynoldssche Ausdruck t = — gu’v’ ergibt dann mittels der genannten Hypothese aus di- 

du 

oe 


: 2 ee du 
mensionsanalytischen Griinden sofort t = @/? fe re Das Vorzeichen erhalt man mittels der Bemerkung, 


daB die durch den Mischungsvorgang erzeugte scheinbare kinematische Zahigkeit > 0 ist. Die Lange 1 in 
der Forme! ist als geometrische Funktion anzusehen. 

2 L. Prandtl, Z. angew. Math. Mech, 22 (1942), S. 242. 

3S. z. B. L. Prandtl, Strémungslehre, S. 114, Braunschweig 1949. 
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so schreiben wir demgemaB 


eee ee. see aT 
=) ae A — 
Q 1 dy ? Vee El, ay (13) 
Beziiglich der Querbewegung nehmen wir wieder an 
woke (v1=4,[w). 4 
Wir setzen : ‘ Be 
[| = kel fe 
eigen aan (15) 


(epTY v”| = b\(cpT)'| |v’ . | 
Der Korrelationskoeffizient k ist dabei fiir alle Formen des turbulenten Austausches derselbe. 


Dies ergibt sich daraus, daB8 der Austausch fiir alle Formen an den Transport durch denselben 


Trager gebunden ist. Fiir den turbulenten Massenstrom z. B. erhalt man, wenn L bzw. EL den 
momentanen Mischungsweg bedeutet, 


pe ete sare tes ee 
o’v’ = EL x v" == | 2 | |ELv’| . 


Mit der Korrelation |Lv’ = k|L|- |v’| erhalt man 
Na heenaes 75 eines 
oul = oy | REIL v'l==Klo’| |e. 


Entsprechend ergeben sich die anderen Korrelationen. Der gemeinsame Korrelationskoeffi- 
zient ist also definiert durch 


[Lo'| = IL [v7. (16) 


Mit den Ansatzen (12), (13) und den Korrelationen (14), (15) erhalten wir 
a) fiir den turbulenten Massenstrom 


M=—o= 2k [o'|[e|== Ah Fe 
du||0o du 00 
= thh, ER ||2\= +kh, ERS % 
leh ET aay £kh, El 
du 39 
= + El? dy Oy (k | el be I?) ° 


Offenbar ist das Vorzeichen von — 9’v’ identisch mit dem von 00/0 y; denn bei positivem v’ z. B. 
ist 9’ positiv fiir negatives 09/0y und umgekehrt. Das gibt 


du| do 
ae SF) | eee eens 
M—=El By Oy 
b) fir die turbulente Schubspannnung! 
t= — (guy =— Gu’ fag) =— 0 wv — ug’ 
= + ok|u’| v’| + uk|o’| |v’| = +ekk,|u’'| |u| tukk,|o’| |u’| 
J Te pee 190 | | au! it aa = ae ou 
aot PA hashes pCO Cbs 25. ery. fe as 2 
= +0okk,P By +ukk, El len +0 kali (55) & uhh, il By ay 
OU \Cu OOLOL: 
2 ce 
+ol (35) +tuEkl By Dad 


1 Auch hier ergeben sich die Ansatze fiir den turbulenten Austausch in einfacher Weise aus der H y po- 
these, daBderturbulente Austausch bestimmt wirddurchdenlokalenVer- 
laufdermittlerenVerteilungen. Betrachten wir beispielsweise die verallgemeinerte Reynolds- 
sche Formel t = — (ou’v’ + u 0/v’) fiir die turbulente Schubspannnng. Der lokale Verlauf der mittleren 

du do 
Verteilungen wird hier festgelegt durch und “°. Die obige Hypothese liefert dann aus dimensions- 


dy d 
analytischen Griinden sofort — ou’ =oP 
du| {_ du Ee dole. 
5a (eas + Eu dy{ fiihrt. 


do Le 
ive was mit 7 = E auf 


dy 


oeele 
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Das Vorzeichen von —u’v’ ist identisch mit dem von du/dy; das von —9Q’v’ mit dem von 
do/dy . Das gibt 
du) /_ du _ 00 
RICE om . 
taal ay (2 ay Bus, 


c) fiir die turbulente Warmeleitung 


00 cpT 
dy 
entsprechend bei etwaigen tian bometae care einer Beimengung fiir die turbulente 


Diffusion 


= EP); es 


‘du dec 
D= EP — ; 
|ay| ay 
Der Mischungsvorgang in Feldern mit stark veranderlicher Dichte wird also durch folgende 
AustauschgréBen beschrieben: 


ou, do 
a 2 17 
Mes ay (17a) 
ee mee 17b 
=P (55 is + Ewe). (17b) 
du| dOcpT iL 
pop 28 eben (176) 
C= EP \ay| ay 
p= Er || 2¢ (17a) 
vO 
Fiir kleine Dichteunterschiede bzw. 0 ~ Feat erhalt man wieder die bekannten Formeln 
— zo |u| du 
EES 
Tol By By 
und (abgesehen vom Faktor E) 
~ 79 |u| OT 
—=He..0 1 — 
NEG RN 


Bei gréBeren Dichteunterschieden sind also die Zusatzglieder 


EP a5 poe fiir die turbulente Schubspannung, 
und 
El? 35% C, re fiir die turbulente Warmeleitung 


zu beriicksichtigen. 
Die oben abgeleiteten turbulenten AustauschgréBen beziehen sich, wie noch einmal hervor- 
gehoben werden soll, auf ,,Kontinuitatsgleichungen“, vgl. (4a) bis (4c). 


4, Zusammenfassung, Die statistische Betrachtungsweise nach 0. Reynolds wird auf Felder 
mit stark veranderlicher Dichte verallgemeinert. Auf der Grundlage der Prandtlschen Konzep- 
tion des Mischungsweges wird dann ein Gleichungssystem zur Berechnung der Mischungsvor- 
gange aufgestellt. 


(Eingegangen am 3. Oktober 1951.) 
Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, JagerstraBe 22/23. 
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Turbulente Vermischung zweier ebener Luftstrahle von fast gleicher 
Geschwindigkeit und stark unterschiedlicher Temperatur. 


Von W. Szablewski. 


1, Feldgleichungen fiir die turbulente Vermischung zweier ebener Luftstrahle stark unterschied- 
licher Temperatur, Wir setzen im folgenden den Druck als konstant voraus. Nach Gleichung (17) 
der vorangehenden Arbeit des Verfassers wird dann ein ebenes turbulentes Mischungsfeld mit 
der Hauptstrémungsrichtung x durch folgendes Gleichungssystem bheschrieben: 


_douu  douv O20 ys 0 
Impuls: una é(x) ay (2 By 7 Eu oy) ; (1) 
- . OBcpTu , OepTH GGcnl 
Energie: oe ay SED) ak , (2) 
"Ogu Oor 2. 070 
Masse: sao + ina E e(x) ay?” (3) 


(Hierbei bedeuten « die Hauptstrémungsrichtung, y die Querachse, u und v die Geschwindigkeits- 
komponenten in x- und y-Richtung, g die Dichte, T die absolute Temperatur, c, die spezifische 
Warme bei konstantem Druck, E eine empirische Zahl (Ubertragungsverhaltnis von substan- 
tiellen Eigenschaften und Geschwindigkeit beim turbulenten Austausch). Das Uberstreichen 
deutet an, da es sich um die zeitlichen Mittelwerte einer stationéren turbulenten Strémung 
handelt. Im folgenden werden die Querstriche fortgelassen.) Dabei haben wir, einem Vor- 


schlag von L. Prandtl! gemaB, die scheinbare kinematische Zahigkeite = /? = dure Wee 
Mittelwertbildung . 

€ (x)= x, by (x) (Up— uy) (Up > Uy) . (4) 
ersetzt, wo b,(x) die Mischbreite des Feldes der Geschwindigkeit, x, eine dimensionslose empirische 


Konstante*, uy und u, die Geschwindigkeiten der ungestérten Strahle 0 und | bedeuten (Abb. 1). 
Die GréBe c, wird als Konstante behandelt. Wir fithren noch die Ubertemperatur 


= T—T, 
ein und erhalten dann fiir (2), wenn wir (3) beriicksichtigen, 
dgué , dovd 008 
af 


= E¢(x) (2’) lt 


Ox 0 Oy? | 
Mittels (3) kénnen wir unser Gleichungssystem (1), = 
(2’), (3) auch schreiben : 


3 a au du do/dy = : 
wget gy el) [a A te a pale (5) Sa 


ao ao 029 08 do/dy aaa 
so le ee SSE elle , 6 
u Ox v oy é(x) ly | dy 0 | ( ) SA 
-s 
dou . dov 0"0 0 
<— a = = Ee(x) dy?" (7) Abb. 1. 
Mit der Zustandsgleichung fiir ideale Gase bei konstantem Druck 
o [—konat: (8) 
ersetzen wir 9 durch #: 
konst. _ konst. il (9) 


CS Ace ee Te OTs 


1 L. Prandtl, Bemerkungen zur Theorie der freien Turbulenz, Z. angew. Math. Mech, 22 (1942), S. 242. 
2 Der empirische Koeffizient x, wird fiir das Phanomen der Vermischung ebener Luftstrahle als Kon- 


stante angesehen. Unter der Voraussetzung, dab |—— nur wenig variiert, deckt sich die obige Annahme 
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und erhalten dann 


ao/T, 
. ~ eu Ou tu Qu oy |, 10 
Impuls: Hae Ae Pa &(x) dy? (E 1) dy 1+ 9/T, : { 
LAS 
DP PTOS Oo ad aa Dade: ll 
Energie : Ua, cle! ag eee) dy? Ov 14+- WT, Oy 


Fiir die Masse ergibt sich zunachst aus (7) 


| ee | 
du , dv 1 Oo/T, | ad/T, | Oy el tye 1078) Tai 
Sas Gy, ita, |" ye vie Oyhag BHO a ay, Oy? J- 


Hier ist aber nach (11) die rechte Seite der Gleichung identisch Null, so da8 wir erhalten 

du , Ov 

ee ee eee 12 
pos, 0. (12) 


Infolge des Fehlens einer kennzeichnenden Linge ist es naheliegend, Affinitat der Profile iiber 
der Mischbreite anzunehmen. Aus dimensionsanalytischen Grinden ergibt sich dann fiir die 
Mischbreiten 

b~ x, b,~x, (13) 


wo b, die Mischbreite des Temperaturfeldes ist. Wir transformieren auf die Ahnlichkeitskoordinate 


¥y . 
i foes | 


mit (14) 
UN nodes ed Re Cie alee 
(a5) bate die (55), hee le 


Das gibt, wenn wir gleichzeitig mit ug, 09(—= T)— T,) dimensionslos machen, 


d9/T, 
_ @uluy , du/uy (: 2a | x | di a 
Impuls: dn? - dn \ uo 1] Up &(x) ( | 1) 1+ 0/T, == > (15) 
d9/T, 
., BO,  dd/O| 1 (= DNB dye ike. 
Energie: iy? ‘e Fis a a! | MO L407, 0, (16) 
‘ du/uy dv/uy __ 
Masse: a a ae =) yh, (17) 
Hier ist mit a = ¢, nach (13) 
€,(x) Up —U 
at = Ky, Cy Fi uy 
Fihren wir nach dem Vorbild von H. Gériler! in die Impuls- und Massegleichung die Variable 
il i 8 
oS Gta oe ee ee 
Oy smit-o°—— 2, 0, SG) (18) 
Uo 
und in die Energiegleichung die Variable 
1 1 
it es — 1 
0x2 mit o4= > ou, e, on (19) 
Uo 
ein, so lautet unser Gleichungssystem schlieBlich 
: Puy du/uy uae v E-- 1 @/T,| _ 
Le erat tk eps 2 e oti a5.) 1+9/T, don | wees a 
. 23/9, , dO [> /u » 2 dT 
E > 4 4 ° Se) Ae \— 
pide d(o,n)? | doy7 2 f ol ack zl 1+6/T, do, 1 wh cu 
; du]ug dov/uy 
Masse: Co ae f dan (22) 


1 H. Gértler, Z. angew. Math, Mech. 22 (1942), S. 244. 
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Das sind die Gleichungen, welche die turbulente Vermischung zweier ebener Luftstrahle stark 
unterschiedlicher Temperatur beschreiben. Gegeniiber dem Feld konstanter Dichte treten in den 
Differentialgleichungen die unterstrichenen Glieder zusatzlich auf. 


2. Berechnung des Mischfeldes fiir den Fall (uj—u,)/uy) <1. Wir beschrinken uns im folgenden 
auf die Behandlung des Falles, daB, bei beliebig groBen Temperaturunterschieden, die Strahle 
fast gleiche Geschwindigkeit besitzen: 

pean 


0. (23) 


Up 
Offenbar kénnen wir dann in den Gleichungen (20) und (21) 

ual — a y (24) 
setzen. Denn mit (u)— u)/u)—> 0 strebt nach(4)auch die scheinbare kinematische Zahigkeit e— 0. 
(Zu denselben Gleichungen gelangt man, wenn man fiir u/up, v/uy und 3/3, einen nach Potenzen 
von (U)— u;)/Uy fortschreitenden Reihenansatz macht und die héheren Potenzen von (up— u,)/Up 
vernachlassigt.) 

Wir erhalten dann das vereinfachte Gleichungssystem 


_ d®u/up du/ug (E+1) d3/T,] __ 
Ce deat sadn, 2 1+0/T, don as (25) 
.. &/9,  dd/9y | 2 do/T;) 
Energie: dc)? eek 204% TFO/T, down ==); (26) 
: dujuy , dov/uy _ 5 
Masse: oY og utdan —= (5 (27) 
Die Integration dieses Gleichungssystems ]aBt sich bis auf Quadraturen durchfiihren. 
Fiir die Temperaturverteilung #/3) mit den Randbedingungen 
aes ==€© 
OO 5> i fal amome 1 a 
erhalt man aus (26) 
0/8, = — : 
| 0 &/T, |aP(o,n) +B 
mit der Fehlerintegralfunktion 
Ox7 
2 ieee oy)? A 
P(o, 4) = ix | e—(%%")* doy 4 (28) 
0 
und den Integrationskonstanten 
ie eat ee iA) 
{go jiaey wale: Bes zi 120 O/T, 
Fiihren wir in (25) zweckmaBigerweise die Variable o,7 statt oy ein, so kommt 
Cele 
BUG sie eg 8 Tay gy = 0 (29) 
du/ug 14+0/T, doy 
do, 
mit den Randbedingungen | ae 
Ca athe 47 : 
Ug— Uy Vo ate! | 
| Man erhalt 
u u a ene d ir) 
SE ey, - aim Oy? 
Ug— Uy } | (x B(o,.) +B)" ** eat 
0 
: 30 
mit den Konstanten hg (30) 
oo 
a) a) 
So Viel 
0 0 0 0 
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und schlieBlich die Querkomponente 


—E(o4n)* , 
0% v/Uo fe 28 G Ege do, x + konst “4 (31) 


LY Oneemeg 

onal 
Die Integrationskonstante in (31) kénnte man nach einer Bemerkung von v. Kdrmdn' aus der 
Gleichung 

QU Yo t 01% = 9 

(v, bzw. v, asymptotischer Betrag der Querkomponente des Strahls 0 bzw. 1 im Mischfeld) er- 
mitteln, die besagt, daB der Impuls in der Hauptstrémungsrichtung infolge des Fehlens auBerer 
Krafte keine Anderung erfahrt. 

Eine andere Bestimmungsgleichung ergibt sich auf Grund folgender Bemerkung: Bei der 
Ausbreitung eines aus einer Diise tretenden Luftstrahls in bewegter Luft tritt das Phanomen der 
Vermischung zweier ebener Luftstrahle in unmittelbarer Nahe des Diisenrandes auf. In dem 
durch die Vermischung noch nicht berithrten Teil des Diisenstrahls, dem Strahlkern, ist das Ver- 
schwinden der Querkomponente anzunehmen. Legen wir dem Diisenstrah] die Anfangsdaten 
Up, T) bei, so haben wir damit die Randbedingung 

ae ee eet OyN>— oO. (33) 
Up— Uy 

Se 
In Hinsicht auf spatere Untersuchungen iiber die Ausbreitung von Heifluftstrahlen legen wir 
(33) der weiteren Rechnung zugrunde. Das gibt 


Ox 7] s —oo 
OxP/Uo i PSA RAMETU pe [ | (34) 

ee cn a | 

ELO 
Fiir kleinere J)/T, ~0 vereinfachen sich die Formeln (28), (30) und (34) erheblich: 
1 i 
9]8=— 4 (qn) +4. (35) 
Ti 1 1 

ene cece oe Ol) stow: (36) 


Die Verteilungen der Temperatur und Geschwindigkeit werden hier durch die Fehlerintegral- 
funktion dargestellt. Da nach (19) die Streckungsfaktoren o und o,, im Verhaltnis < == jE 


stehen, so ergibt sich das Mischbreitenverh4ltnis 


b, = 
a =] E (37) 
Die Querkomponente wird hier 
ovjup 1 1 — (on? 
est 2 \/x : 3 (38) 


en 
Zur vollstindigen Festlegung des Strémungsfeldes ist noch die Bestimmung der Streckungs- 
faktoren o und o,, erforderlich. Das kann so geschehen: Aus der berechneten Geschwindigkeits- 
verteilung (30) laBt sich offenbar o,,c, ablesen (wo nach (13) cj=6,/xist). Es folgt dann aus (19) 


1 il 
SPR ee (39) 


Hier tritt der Faktor x, auf. Dieser stellt einen empirischen Koeffizienten der Theorie dar und ist 
aus Messungen zu ermitteln. Wir gehen an dieser Stelle nicht weiter darauf ein. Es verbleibt 
noch der empirische Faktor E, der das UbertragungsverhAltnis von substantiellen Eigenschaften 
und Geschwindigkeit beim turbulenten Austausch angibt?. Nach L. Prandtl? setzen wir E—2. 


“1A, Kuethe, Journ, appl. Mech. 2 (1935), S. 87. 
2 Vel. W. Szablewski, vorangehende Arbeit. 
3 L, Prandtl, Fiihrer durch die Strémungslehre, S. 114. Braunschweig 1949, 
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Tabelle 1. 


0,055 
0,2 0,345 | 0,185 03123 0,057 0,028 —0,2 0,655 0,452 0,343 0,188 0,102 
0,4 OFZ 0.102 0,065 0,028 0,014 —(,4 0,788 0,611 0,500 0,307 0,181 
OSOmmia OT US 0,050 0,031 0,013 0,006 —(),6 0,885 0,756 0,663 0,463 0,300 
0,8 0,055 0,022 0,013 0,005 | 0,003 0,8 0,945 0,866 0,805 0,636 0,460 
1,0 0,023 0,009 0,005 0,002 0,001 —1,0 0,977 0,944 0,905 0,794 0,643 
Ne: 0,008 0,003 0,002 0,001 | 0 —1,2 0,992 0,978 0,962 0,905 0,807 
1,4 0,003 0,001 | 0,001 0 —l,4 0,997 0,993 0,987 0,964 0,918 
1,6 0,001 0 0 —1,6 0,999 | 0,998 0,996 0,989 0,972 
1,8 |} 0 —1,8 al 0,999 0,999 0,997 0,992 
2,0 —2,0 1 1 0,999 0,999 
23 =2)2 | rei 1 
2,4 =— 24 
2,6 | —2,6 | | 
2,8 | | 32.971 | 
3.0 | | 3 | 
8/9, 
on |8/T,=0| 1 | on iam 10 | oxy ,/T, = 0 1 oe las 10 
— SSS 
0 | 0,500 0,333 0,250 0,143 | 0,083 
0,2 0,389 0,241 0,175 0,096 0,055 —(),2 0,611 0,440 0,344 0,208 0,125 
0,4 0,286 0,167 0,118 0,064 0,035 —0,4 0,714 0,555 0,454 0,294 0,185 
0,6 0,198 0,110 0,076 0,040 0,022 —0,6 0,802 0,669 0,574 0,403 0,269 
0,8 0,129 0,069 0,047 0,024 0,013 —0,8 0,871 Osa 72 0,692 0,530 0,380 
1,0 0,079 0,041 0,028 0,014 0,008 —1,0 0,921 0,859 0,796 0,661 0,516 
ip? 0,045 0,023 0,015 | 0,008 0,004 —1,2 0,955 0,919 0,876 0,780 0,659 
1,4 0,024 0,012 0,008 | 0,004 0,002 —],4 0,976 0,953 0,932 0,872 0,788 
1,6 0,012 0,006 0,004 0,002 | 0,001 —1,6 0,987 0,977 0,961 0,933 0,884 
1,8 0,005 0,003 0,002 0,001 § 0,001 —1,8 0,995 0,989 0,985 0,968 0,943 
2,0 0,002 0,001 0,001 0 0 —2,0 0,998 0,995 0,993 0,986 0,975 
2 0,001 0 0 | —2,2 0,999 | 0,998 0,997 0,994 0,990 
2,4 0 —2,4 iL | 0,999 0,999 0,998 0.996 
2.6 —2,6 1 1 0,999 0,999 
2,8 | —2,8 1 1 
3,0 | 9 | 
Ox0/U, 
) * 
o4n =| 0,/T, = 0 | i 2 | Boy tem IO oy $,/T, =0 | 1 2 5. pd ae 
0 | 0,200 | 0,339 0,444 0,643 0,815 | 
0,2 | 0,184 | 0,329 0,436 0,638 0,812 —(),2 | 0,184 0,324 0,431 0,630 0,811 
0,4 0,145 0,304 0,419 0,630 0,808 —0,4 ; 0,145 0,270 0,384 0,594 0,787 
0,6 0,097 0,279 0,402 0,622 0,804 —0,6 0,097 0,204 0,302 0,516 0,727 
0,8 0,055 0,259 0,390 0,617 0,802 —0,8 0,055 0,128 0,203 0,394 0,614 
| 1,0 0,027 0,247 0,383 0,614 0,800 —1,0 0,027 0,067 0,114 0,253 0,450 
, ie 0,011 0,241 0,379 0,613 0,800 —1,2 0,011 0,030 0,053 0,132 0,269 
1,4 0,004 0,238 0,378 0,612 | 0,800 —1,4 0,004 0,011 0,020 0,055 0,127 
1,6 0,001 | 0,237 ORT 0,612 0,799 —1,6 0,001 0,003 0,006 0,018 0,046 
1,8 0 0,237 0,377 0,612 0,799 —1,8 0 0,001 0,002 0,005 0,013 
a0) 0,237 0,377 0,612 0,799 —2,0 | 0 0 0,001 0,003 
2.2 —2,2 0 0 
2,4 244 | 
2.6 —2,6 
2,8 | —2,8 
3.0 5.0, | 


Derselbe Wert ergibt sich auch aus den Messungen von 0. Pabst! an Heifluftstrahlen, wie in 
einer spateren Arbeit gezeigt werden soll. 


10. Pabst, U. u. M. 8004 (1944) und U. u. M. 8007 (1944). 
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An Hand der Formeln (28), (30) und (34) haben wir die Geschwindigkeits- und Temperatur- 
verteilungen fiir die Parameterwerte J/T,;—0, 1, 2, 5, 10 berechnet. Die vorkommenden In- 
tegrale wurden dabei approximativ berechnet. bole 

In Abb. 2 bis 4 haben wir die Verteilungsfunktionen iiber o,7 aufgetragen. Die in Abb.2 
und 3 gestrichelt eingetragene Kurve verlauft parallel zur Verteilungskurve fir 3,/T,~ 0, die 


nach (35) und (36) durch 
die Fehlerintegralfunktion 
dargestellt wird. Man sieht, 
daf die bei gréBeren Dichte- 
unterschieden auftretenden 
Verteilungen nur wenig von 
der Fehlerintegralfunktion 
abweichen. Abb. 4 zeigt den 
Verlauf der Querkompo- 
nente der Geschwindigkeit. 
Entsprechend dem Um- 
stand, da wir im Mi- 
schungsraum bei fast glei- 
cher Geschwindigkeit eine 
geringere Dichte haben als 
im Strahl1, tritt in diesem 
ein Abstrom vom Mi- 
schungsraum weg auf. 

Es zeigt sich demnach 
aus Abb. 2 und 3, daB die 
Mischbreiten sich nur wenig 
mit dem Dichteunterschied 
andern}. Dagegen tritt, wie 
aus den Abbildungen er- 
sichtlich, folgender charak- 
teristischer Effekt auf: 
Wahrend bei 3,/T,~0 der 
Winkelraum der Mischzone 
symmetrisch zur x— Achse, 
also symmetrisch zu beiden 
wane Strahlen, gelegen ist, dreht 
geen er sich mit wachsendem 
; Temperaturunterschied in 
a markanter Weise zum hei- 
a 18 I] =10 Beren Strahl hin. Man kann 
| auch sagen, daB bei der 
S =06|—— —— = Ausbreitung von Heifluft- 
strahlen (die aus einer Diise 
ee ay Ca Sd EE a: treten) in Luft von fast 
gleicher Geschwindigkeit 
Q2 ORES evar pessoa der Innenwinkel der Misch- 
| zone mit der Temperatur- 
0 0,4 differenz wachst, waihrend 
wa ihae! Fook bani iad “OS 0 Os 10 15 20. der Aufenwinkel kleiner 

Abb. 4. Verteilung der Querkomponente der Geschwindigkeit. wird. Wir haben diese Ver- 

haltnisse noch besonders in 

Abb. 5 dargestellt, welche die Winkelraume der Mischzonen fiir das Geschwindigkeits- und 

Temperaturfeld zeigt. Die Grenzen der Mischzonen wurden dabei durch (u—uw,)/Uy bzw. 0/4 
gleich 0,99 und 0,01 festgelegt. 

Aus dem geschilderten Verhalten kann man den Schluf ziehen, da® bei der turbulenten Aus- 
breitung von HeiBluftstrahlen in Luft von fast gleicher Geschwindigkeit ((up—u,)/uy ~ 0) heiBere 


Abb. 2. Verteilung der Lingskomponente der Geschwindigkeit. 


HY, 


10 


Abb. 3. Temperaturverteilung. 


* Es ergibt sich zunaichst aus den Abb. 1 und 2, da® o,c; sich nur wenig mit dem Dichtenunterschied 
aindert. Nach (39) trifft diese Feststellung dann auch fiir c, und damit nach (19) auch fiir c,=b,/x zu. 
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Strahle eine kiirzere Reichweite (gemessen durch die Lange des Strahlkerns) haben werden. Man 
kann annehmen, daB dieser Effekt ganz allgemein (fiir beliebiges (uj—u,)/uy) bei der Ausbreitung 
von HeiSluftstrahlen in bewegter Luft in Erscheinung treten wird. Diese Vermutung findet sich 
durch die Messungen von 0Q. Pabst! bestatigt, wie in einer spateren Arbeit gezeigt werden soll. 


3. Priifung der Annahme konstanten Drucks, Wir haben unsere Rechnungen unter der An- 
nahme konstanten Drucks durchgefiihrt. Nach dem Vorgang von W. Tollmien? wollen wir diese 
Annahme an Hand der zweiten Bewegungsgleichung nachpriifen. 

a) Die Bewegungsgleichung in der Querrichtung lautet 


dgvu , Oovv _— ap 4 at | Qoy 
dx | dy dy ' Ox ' Oy’ 
wo (40) 
Ou 00 
t=c(x) ( ay +E 5) 
Ug-u 
two 
9g Oy/, T,=0 
Ony Geschwindigkertsteld %/Ty=-0 Oxy — emperaturreld 
7 7 : 
Wy) iy 70 
a 
10 
cD, | & 
Y he 05 10 
-7 =i 
Ug» lp 0 Uo, Ip 
(Y=) (W-1-;) 

| 5 0 

| 10 
fy =2 5 
70 


Abb. 5. Winkelraum der Mischzonen. 


und og, die scheinbare Druckspannung in der y-Richtung ist. Beachtet man (7), so kénnen wir 
(40) auch schreiben 


du OS OO a du | 00 Ties ae Ooy ; 4] 
Qua Cay dy I ax | €) Q dy I aay B(x) OF: ate dy ( ) 
Transformiert man auf 7 = y/x (s. (14)) und macht dimensionslos, so wird 
@ dvjuy /u B\ 1 dp 1 [Q du/uy L (E- 1) “leo du/ Ug LE u a 
Qo dy \Uo Up Qouy dy 207/09 dn® dy dn Up dy (42) 
Ev @Q/e | I do, | 
2o%u, dry? | ou2 dy 


Ersetzt man in dieser Gleichung nach (9) 9 durch #, eliminiert die erste Ableitung von v/Up 
mittels (22) und die zweiten Ableitungen von u/u) bzw. 0/9, mittels(20) baw. (21), so erhalt man 


1 OQ. Pabst, a. a. O. 
2 W. Tollmien, Z. angew. Math. Mech. 6 (1926), S. 468. 
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bei Einfithrung der Variablen o,7 (o?= Eo) schlieBlich 
1 OE Lp aD Oe du/ug ( u Maree 
4 (0,/2)u2 1+ 0/T, do,n 14+ 8/T, do, 7 \uo 1] 2 


1 1 do/T, u \2 . y \2 il OR ite apse Z doy | 
2 (1+0/T;) do.n fe 541) * ug) |! © (e9/2) ug 1+ A/T, don’ 
wo nach (18) o, = —— = ae 


(43) 


1 1 


22), 


b) Beschranken wir die Betrachtung auf Strahle fast gleicher Geschwindigkeit, so haben wir 
nach (42) 


v 
Ox —Ox rr 0x7] - 


Damit vereinfacht sich (43) zu 
1 o% 1 dp _—(Gm)*_ du/ug (647)? d0/Ty | 1 _ % : Be (44) 


1 
A (09/2) u3 109] T, dogn  1+0/T,doyn 2 (1+ 9/T,)? doyn | 4 (go/2) uj 1+ %0/T; doar 
Nach (30) ist 


—E(o,)? 
du/uy — Uy— Wy e (Gal 


do, Up y (x D(a n)+B)E+} é 


nach (28) 
=1 dd/T, 2 


(1+8/T,)? do. x 


e—(Fxn)* 


Damit erhalten wir 


fo. 3) 


Pe eres aN ema Gem) ons | 90 (45) 
(9/2) ug Oe = ACL Pol Ta) | (a Deyn TBE OT, 
Fir 0/T,=0 ergibt (45) 
Piosica Meroe ise Se res aad 46 
(oo/2)u2 ° uo me 


also Uberdruck im Innern des sich in bewegter Luft (Strahl 1) ausbreitenden Strahles (Strahl 0). 
Diese Formel, die iibrigens, wie man zeigen kann, fir beliebige (uy—uj)/Uy gilt, hat bereits 
H., Gortler? abgeleitet. 


Fiir nicht zu kleine %)/T, itberwiegt wegen a 0 das zweite Glied in (45): 
0 


Po P-« 9 ~ Do 47 
(eo/2) ug *  , * a 
Fiir nicht zu kleine Temperaturdifferenzen haben wir also einen der Temperaturdifferenz propor- 
tionalen Unterdruck im Innern des sich in Luft fast gleicher Geschwindigkeit (Strahl 1) ausbrei- 
tenden Strahls (Strahl 0) zu erwarten. 
Wegen der groBen Betrage von o? baw.o2 (im Fall Jy/ T,= 0, (up — uy)/Up—=1,0 z. B.ergibt sich 
nach Gértler ae ai ——0,0055; im iibrigen s. Formel (39)) sind die auftretenden Unter- 
= C0 0 
bzw. Uberdriicke sehr klein. Thre Vernachlassigung erscheint also gerechtfertigt. 


4, Zusammenfassung, Es wird das durch die turbulente Vermischung zweier ebener Luft- 
strahle von fast gleicher Geschwindigkeit und stark unterschiedlicher Temperatur erzeugte Feld 
der Geschwindigkeit und Temperatur berechnet. Als charakteristischer Effekt ergibt sich eine 


mit der Temperaturdifferenz wachsende Drehung des Winkelraumes der Mischzone zum heiBeren 
Strahl hin. 


(Eingegangen am 3. Oktober 1951.) 


Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, JagerstraBe 22/23. 


Deutsche Akademie der Wissenschaften, Institut fiir Mathematik, 
Abteilung ,,Angewandte Mathematik‘‘, Berlin W 8. 
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Das polarisationsoptische Verhalten der Kunststoffe Celluloid und Dekorit“. 
Von A, Schafer. 


1, Einleitung. Fiir polarisationsoptische Versuche, die nicht im rein elastischen Gebiet 
durchgefiihrt werden, ist die Kenntnis der sich im Inneren eines belasteten Modellwerkstoffes 
abwickelnden Vorgange wichtig. Eine Deutung dieser Vorgiinge kann aber erst erfolgen, nach- 
dem eine Erklarung fiir die unter der Belastung auftretende Doppelbrechung gegeben ist. Es 
wird danach mit diesen Kenntnissen Dekoritmaterial verschiedener Hartungsgrade untersucht. 
Folgende Fragen werden beantwortet: 

1. Was bewirkt die Hartung im Innern des Materials ? 

2. Was ist bei der Erstarrung raumlicher Spannungszustinde! zu beachten ? 

Die Ergebnisse der vorliegenden Arbeit werden aus zahl- 
reichen Versuchsreihen mit Celluloid und Dekorit gewonnen, 
um Fehlschliisse durch nicht vermeidbare MeBungenauigkeiten 
auszuschalten. 


2, Allgemeines zu nachfolgenden Versuchen, Bei der Doppel- 
brechung durchsichtiger, amorpher Stoffe unter Spannung gilt, 
wie bekannt? 


d= = (0, — 0) 4. (1) 


Hierbei ist 6 der Gangunterschied zwischen auberordentlichem 
und ordentlichem Strahl, C eine spannungsoptische Konstante, 
A die Wellenlange des verwendeten Lichtes, (¢, — 0,) die Haupt- 
spannungsdifferenz, d die Dicke der Scheibe. Die folgenden 
Versuche werden an Zugstaben durchgefihrt. Fiir den dann 
an allen Stellen des Modells gleich groBeneinachsigen Spannungs- 
zustand lautet (1) 

b= od. (2) 
Besser ist es, wie noch gezeigt wird, mit der Dehnung, statt 
mit der Spannung zu rechnen. Mit o = ¢,E wird (2) 


6= 5 Fed, (3) 


wobei ¢, die Dehnung des rein elastischen Zustandes — auch 
Lastdehnung genannt — ist. In den folgenden Versuchen treten ' 
noch Zeitdehnungsanteile, die Kriechdehnung und die plastische Abb. eee See teh 
Dehnung auf. Mit Kriechdehnung soll die nach der Entlastung 

nicht bleibende Dehnung, die sich beim Kriechen des Materials ergibt, bezeichnet werden. Die 
plastische Dehnung ist ein Dehnungsanteil, der bei sehr hohen Spannungen auftritt und eine 
bleibende Deformation zur Folge hat. 


3. Versuche mit dem Modellwerkstoff Celluloid. a) Der Aufbau und die Durchfithrung 
der Versuche. Ein Zugstab von 5 mm Starke hangt mit seinen breit ausgefiihrten Képfen im 
Versuchsrahmen. Die Belastung erfolgt durch Auflegen von Gewichten an einer Laststange. Die 
Dehnung wird mit zwei Huggenbergtensometern gemessen. Der Versuchskérper steht zwischen 
gekreuzten Polarisatoren; als Lichtquelle dient eine Natriumdampflampe 3, Die Abb.1 zeigt den 


* Die Untersuchung wurde im Mechanisch-technischen Laboratorium der Technischen Hochschule, 
Miinchen durchgefiihrt. — Dem Leiter des Institutes, Herrn o. Prof. Dr. L. Féppl, habe ich fiir die stete 
Foérderung der Untersuchungen zu danken. 

1 E. Ménch, Ing.-Arch. 16 (1948), S. 267. A 

2 L. Féppl und H. Neubert, Festigkeitslehre mittels Spannungsoptik. Miinchen 1935. 

3 L. Féppl und R. Hiltscher, Bauing. 20 (1939), 5. 231. 
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Versuchsrahmen mit Zugstab, der durch Gewichte belastet wird; zu sehen ist noch der Polari- 
sator und ein Teil des Lichtkastens als weiBer Streifen. 

Im unbelasteten Zustand erscheint der zwischen den gekreuzten Polarisatoren stehende Stab 
dunkel. Mit steigender Last tritt im mittleren Stabteil eine gleichmafige Aufhellung und Ver- 
dunklung ein. Der Ubergang von einer Verdunklung zur nachstfolgenden entspricht einer Iso- 
chromatenordnung. Die Zeitpunkte der starksten Abdunklungen lassen sich sehr genau fest- 
stellen. Ist die gréBte Abdunklung erreicht, so wird die Dehnung und die Last abgelesen. 


b) Die Ergebnisse der Celluloidversuche. Die Versuche zeigen, daB bei konstanter 
Last von einer bestimmten Spannung ab mit der Zeit ein merkliches Kriechen und damit eine 
merkliche Helligkeitsanderung des Stabes in der spannungsoptischen Einrichtung eintritt. Zu 
einer bestimmten Zunahme des spannungsoptischen Effektes bei konstanter Last gehért nicht 
mehr die Dehnungszunahme, die bei einer Laststeigerung zur Erreichung des gleichen optischen 
Effektes erforderlich ist. Zu diesem Ergebnis fiihren auch Versuche, die mit verschiedenen Be- 
lastungsgeschwindigkeiten durchgefiihrt werden. Als Beispiel zur Anderung der Belastungs- 
geschwindigkeiten seien hier die Versuche 4 und 5 der Zahlentafel 1 angefiithrt. Beim Versuch 
Nr. 4 betragt die Belastungszeit von der 3. zur 4. Isochromatenordnung 20 Minuten; beim Ver- 
such Nr. 5 nur 2 Minuten. Die Tabelle 1 zeigt, da®B der Isochromatensprung von einer Ord- 
nung durch eine bestimmte Last- und eine bestimmte Kriechdehnung erreicht wird. Die Deh- 
nungsanteile sind, da ihre optischen Wirkungen verschieden stark sind, mit entsprechenden Fak- 
toren zu versehen. Wird in Gleichung (2) fiir C/A a* und d= 1 gesetzt, so ist bei Beriicksich- 
tigung der Kriechdehnung ,,¢,‘* der optische Effekt 


6= a* og + be, (4) 
oder nach (3) bei eels: a 
6= aée,+ he, (5) 
In Zahlentafel 1 wird a= 1 gesetzt; dann ist b= 1/7,5. 


Die Summe aus beiden Dehnungen ist, wenn die Faktoren entsprechend bestimmt werden, 
fiir jeden Isochromatensprung konstant. 


Tabelle 1. Der EinflufB der Last- und Kriechdehnung auf den optischen Effekt. 


Dehnungen in 1/400 mm bei 20 mm MeBlange 


Vers. ton Pend : : u g 
Nr. matensprung sprung se eenO Cer eh asted Kriechdehnung =445 
i, Ordnung i. kg/em? Geek Gy, Ey €,,/ 7,5 & + 1/7,5 e, 
1 von 0 n. l 95 29 29 0 29 
2 pot eee os 91 36 28 aba 29 
3 EK ees 82 54 25 29/7,5 = 4 29 
4 5 ating, 58 101 18 if US) == i 29 
5 56 gn! 88 4,2 27 UY Obed 29 


4, Versuche mit dem Modellwerkstoff Dekorit, a) Allgemeines und die Grundlagen zur 
Versuchsdurchfiihrung. Die Versuche ergeben, da8B beim Kunststoff Dekorit 6 nicht kon- 
stant ist. Die Abnahme des 6-Wertes bei héheren Ordnungen deutet daraufhin, daB bei diesen 
Laststufen ein Dehnungsanteil vorhanden ist, der nur von einem geringen oder von keinem op- 
tischen Effekt begleitet wird ; zeitlich — bei konstanter Last — sind ebenfalls Schwankungen vorhan- 
den. Kin genauerer Aufschlu8 dariiber wird erhalten, wenn die Anderung des optischen Effektes 
mit gréBerer Genauigkeit und gréBerer Geschwindigkeit festgestellt werden kann. Diese For- 
derung erfiillt die Kompensation mit Hilfe eines //4-Plattchens, das zwischen Modell und Analy- 
sator zusatzlich angebracht wird’. Die Ausgangsgleichung, die fiir die Versuchsdurchfihrung 
maBgebend ist, ist Gleichung (5). An Stelle der Kriechdehnung «,, in (5) soll die Zeitdehnung « 
treten, zur Kennzeichnung, da folgend auBer der Kriechdehnung noch weitere Zeitdehuunce, 
anteile zu beriicksichtigen sind; fiir 6 ist dann b, zu setzen. Eine Zunahme des optischen Effektes 
Ad unter konstanter Last wird durch eine Zunahme der Zeitdehnung Ae, bewirkt: 


5+ A6= ae, +b, (e, + Ae,). (6) 


1G. Mesmer, Spannungsoptik, $.59. Berlin 1939. 
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Nach Gleichung (5) und (6) ist Adi== b, Ac 


die Zunahme des optischen Effektes wird in Winkelgraden am Analysator gemessen; deshalb soll 
fir 4d nun Aa gesetzt werden: 


. 
bah 


Ax = b, A 7 Deen 
x he ZANE oder he yer (7) 


b, ist die Zunahme der Ordnung, bezogen auf die Zeitdehnung in Prozenten, fiir einen Stab von 
lem Dicke. ' 

Ks 1a8t sich auch ein a fiir die Belastungsdehnung definieren. Eine Zunahme des optischen 
Effektes wird hierbei ausschlieBlich durch eine Zunahme der Belastungsdehnung Ae, hervor- 
gerufen. Mit Gleichung (5) ist 


6+ Ad=ae,+ ade, ; 


daraus folgt Aa = ade, oder = a P (8) 
LNEG 

Der Wert aist die Zunahme der Ordnung, bezogen auf die Lastdehnung in Prozenten, fiir einen 

Stab von lcm Dicke. 

Es werden hier nur solche Ergebnisse ausgewertet, bei denen die Dehnungen an den beiden an- 
gesetzten Tensometern immer die gleichen Wertzunahmen zeigen. Diese Bedingung setzt ein 
einwandfreies Material und eine sehr gute Bearbeitung voraus. 

b) Die Zeitdehnung in Abhangigkeit von der Belastung. Die Belastungshéhen 
werden so groB gewahlt, das ein merklicher Kriecheinflu® vorhanden ist. Die folgenden in Ta- 
belle 2 zusammengestellten Versuchsergebnisse zeigen, daB mit zunehmender Last die Anfangs- 
b,-Werte bis zu einem gewissen Grenzwert abklingen, um dann mit der Zeit alle einer gemein- 
samen tieferen Grenze (b~10), die durch die Kriechdehnung gegeben wird, zuzustreben. 


Tabelle 2. Die Abnahme des optischen Effektes — bezogen auf die Zeitdehnung — in 
Abhdngigkeit von der Spannung und von der Zeit. 


An ° 
o b, = Ae, %* 100 Bei konstanter Spannung nach ungefahr folgenden Zeiten in Minuten. 
in kg/em? 
3 | 6 9 12 15 18 21 

10,8 21,0 21-0205 — |= — — 

21,6 18,2 Lowen | eeelOss 16,2 == = = 

37,8 15,1 T3707, || 11258 ES. 11,2 10,6 10,4 

54,0 14,7 13,2 12,0 11,0 10,5 10,4 10,4 

81,0 14,3 12,8 11,7 10,8 10,4 10,3 10,3 
108,0 14,1 12,3 10,8 10,4 10,4 10,3 10,3 


Aus Tabelle 2 kann gefolgert werden, daf ein nach einer gewissen Zeit sich auskristalli- 
sierender Vorgang zu Anfang der Zeitmessung noch von einem oder mehreren Vorgangen, die eine 
b,- WertvergréBerung hervorrufen, iiberdeckt wird. 

Der sich bei der Erhéhung der Spannung ergebende a-Wert wird laufend kontrolliert; er ist 


bei allen Laststufen: a— Ax°/Ae,% 100 = 42 (8a) 


c) Der Zeitdehnungsverlauf bei einer bestimmten Belastung. Tabelle 2 sagt 
noch sehr wenig iiber die Vorgange im Innern des Materials aus. Um genauere Aussagen hieriiber 
machen zu kénnen, wird in kleinen Zeitabstanden (Bruchteile einer Minute) die Geschwindigkeit, 
mit der sich die Zeitdehnung und die Ordnung dndert, bestimmt. Die Versuchsergebnisse wurden 
durch zahlreiche Versuchsreihen erhartet und bestatigt. 

Untersucht werden Materialien im Anlieferungszustand, Materialien, die kurz vor der Messung 
einer Belastung unterworfen wurden — ermiidete Materialien — und Materialien, die beansprucht 
waren, aber durch langeres Lagern sich wieder erholt hatten; letztere Materialien standen nicht 
mehr unter den Nachwirkungen der vorherigen Belastung. Untersucht werden ferner Materialien, 
die vor der Messung einer bestimmten Spannung ausgesetzt waren und spater daraufhin bei einer 
Wiederbelastung bis zu dieser Spannungsgrenze die Merkmale der ermiideten Materialien und 
oberhalb der Grenzen die Merkmale der erholten Materialien zeigen. 

Das Ergebnis dieser Versuchsreihen ist folgendes: Bei entsprechend gewahlten konstanten 
Spannungen zeigt ein erholtes Material in den ersten Minuten des Versuchsablaufs folgendes 
Bild: Nimmt Ae, /At ab, so steigt zur gleichen Zeit Ao/At an und umgekehrt (Abb. 2). Der Kur- 


venzug Aa/Ae,— b, ist keine Gerade mit konstanter Neigung, sondern eine Kurve, die ihre Nei- 
u 
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gung in einem bestimmten Turnus andert. Der b.-Wert schwankt zwischen einem Mindestwert, 
der die GroBe des Beiwertes der Kriechdehnung hat und einem Héchstwert, der die GréBe von a 
erreicht (Abb. 3 Kurve 1). = ei 

Nach langerer Belastungszeit geht der Kurvenzug Ao|/Ae, in eine Gerade mit konstanter Nei- 
gung, die durch den Kriechvorgang gegeben wird, iiber (Abb. 3 Kurve 2 u. 3). 

Bei sehr groBen Belastungen, die in der Nahe der Bruchspannung liegen oder nicht so groBen 
Belastungen aber dann langeren Belastungszeiten ttberlagert sich noch ein dritter Dehnungs- 
anteil, der von keinem oder nur von einem sehr geringen optischen Effekt begleitet wird. 


5. Die physikalische Deutung des optischen Effektes, Bevor die Versuchsergebnisse be- 
sprochen werden, soll eine aus diesen sich ergebende physikalische Deutung des optischen Effektes 
beschrieben werden. Die Doppelbrechung kann an einem Stab, der Translationsschwingungen 
ausfiihrt, anschaulich erlautert werden. Es wird ein schwingender Stab von elliptischem Quer- 
schnitt betrachtet. Entsprechend seiner Steifigkeit in den Hauptrichtungen sind die Fre quenzen, 
in denen er schwingt. Wird er in einer beliebigen Richtung angestoBen, so tritt eine Zerlegung 
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Abb. 2. Der optische Effekt und die Zeitdehnung bei einem Abb. 3. Der optische Effekt in Abhangigkeit von der Dehnung 
erholten Material in Abhangigkeit von der Zeit. bei verschiedenen Belastungs- und Zeiteinfliissen. 


der Schwingungen nach der Richtung der kleinsten und gréBten Steifigkeit ein; beide Bewegungen 
erfolgen unabhangig voneinander, jede mit der ihr eigenen Frequenz. Der noch nicht mit einer 
Doppelbrechung behaftete, unbelastete Kunststoff kann wie ein durch Querschnitte aufgeteilter 
Stab behandelt werden, dessen elliptische Scheiben in vollkommener Unordnung zueinander 
liegen und somit keine ausgezeichnete Richtung mehr besitzt. Erst durch ein Aneinander- 
reihen der Scheiben, wobei auch die entsprechenden Hauptachsen aufeinander zu liegen kommen, 
also durch ein Ordnen, wiirde wieder ein dem elliptischen Stabe ahnliches Verhalten eintreten. 
Die Ordnung wird beim betrachteten Kunststoff durch eine Spannung hervorgerufen, die eine 
Doppelbrechung zur Folge hat. Die Doppelbrechung des Lichtes besteht darin, da® die ankom- 
mende Lichtbewegung in zwei zueinander senkrechte Komponenten zerlegt wird, die sich mit 
verschiedenen Geschwindigkeiten fortpflanzen. 


6, Diskussion der Dekorit-Versuchsergebnisse, Dic Analyse der Versuchsergebnisse lat auf 
drei sich unterscheidende Vorgange im Inneren eines belasteten Materials schlieBen, die sich am 
markantesten durch die Werte Ja/Ae unterscheiden und dessen GréBen unabhangig von der Be- 
lastung und der Belastungszeit sind. Durch Uberlagerung der einzelnen Vorgange ergibt sich bei 
einer nicht entsprechend schnell durchgefiihrten Belastung eine Erniedrigung des a-Wertes durch 
die Einwirkung von Zeiteinfliissen, wie Kriech- und plastische Dehnungsanteile. Bei den Zeit- 
dehnungsmessungen treten ebenfalls Uberlagerungen auf. Die Vorgange selbst lassen sich durch 
eine entsprechende Materialbehandlung, wie Harten und durch eine entsprechende Belastungs- 
héhe und Belastungszeit steuern. Alle Ergebnisse lassen sich erklaren, wenn der optische Effekt, 
wie es unter 5 ausgefiihrt wurde, als eine Ordnung von optischen Grundelementen, die durch die 
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oe dargestellt werden, gedeutet wird, Unterstiitzt wird diese Erklarung noch 
. 2 gendes ersuchsergebnis: Ein in der Belastungsrichtung durchstrahlter Stab zeigt in 

ieser trablrichtung bei der Belastung keine Anderung des optischen Effektes. Von einem sehr 
groBen optischen Effekt wird die Belastungsdehnung begleitet. Die optischen Grundélomonte 
om durch das durch die Belastung hervorgerufene Moment in Richtung der Belastung ge- 
bs u Die GréBe der Drehung ist von der Spannung und der Starke der Bindungen abhangig. 

er beschriebene Vorgang gehorcht dem Hookeschen Gesetz. Bei konstanter Last tritt nun ein 
als Kriechen bezeichneter zweiter Vorgang stark in Erscheinung; dieser auBert sich im Versuch 
durch eine geringere Ordnung gegeniiber Vorgang 1 bei entsprechender Dehnung; die Geschwin- 
digkeit mit der dieser Vorgang ablauft, kann als gering gegeniiber der des Vorganges 1 be- 
zeichnet werden. Erklart werden kann das Kriechen durch das Ausrichten von Ketten oder 
ganzen Verbanden in Richtung der Belastung. Der zur plastischen Dehnung gehérige Vorgang 


ergibt keine Ordnung der optischen Grundelemente; er kann durch eine Translation von Ketten- 
gliedern dargestellt werden. 


Zip Versuche mit Dekorit verschiedener Hirtungsgrade. Die bisherigen Untersuchungen be- 
pe sich auf ein Dekoritmaterial bestimmter Hartung. Folgend soll das rein elastische Verhalten 
(Hooke) verschieden stark geharteter Dekoritmaterialien verglichen werden. Diese Materialien 


verlangen zur Steigerung der Ordnung um eine Einheit ein bestimmtes Ao und bestimmtes Ae, ; 
ihr Produkt ist bei gleicher Dicke d eine Konstante. 


6_, = aAo Ae, = k (bei d = konstant) . (9) 


Die GréBe des Quotienten Ao/Ae, ist ein MaB fiir die Starke der Hartung. 
Bei Beriicksichtigung der Modelldicke d ist 


6_.= ado Ae, = k- (10) 


Die Gleichungen (9) und (10) ergeben sich aus Versuchen, die bei Zimmertemperatur durchgefiihrt 
werden. Versuche bei Temperaturen bis zu ungefahr 70° C zeigen die gleichen Abhangigkeiten 
und auch praktisch die gleiche GréBe der Konstanten a. 

Zur Erklarung der Gleichungen (9) und (10) werden zwei neue Begriffe eingefiihrt. 

1) Fess) ist der effektive Materialquerschnitt, d.h. der Querschnitt, der zum Tragen der Last 
vorhanden ist; also der Querschnitt des gebundenen Materials. 

2) digg) ist die effektive oder optisch wirksame Starke des Materials. Je mehr Material 
gebunden ist, je gréBer ist die optisch wirksame Starke des Probestiickes. 

Bei einem Zugstab, der z. B. ungehartet eine effektive Starke von 2 Einheiten und einen 
effektiven Querschnitt von 6 Einheiten besitzt, soll durch eine Hartung die effektive Material- 
starke verdoppelt werden; bei einer linearen VergréBerung der effektiven Starke andert sich der 
effektive Querschnitt quadratisch; er betragt also nach der Hartung 24 Einheiten. 

Durch die Hartung wurde 
die optisch wirksame Material- “%r/artung : LEE AE LY 


starke verdoppelt, d.h. zur ashranies| (E YY), 2 eee 
Erhéhung der Ordnung um LLL LE Z 


durch Hartung 
eine Einheit ist nur noch die 
halbe Dehnung gegeniiber dem Richtung der Durchstrahling Richtung der Durohstrahling 
ungeharteten Material erfor- Abb. 4. Die effektiven MaterialgréBen. 
derlich; der lastiibertragende 
Querschnitt wurde vervierfacht; die 4uBeren Abmafe blicben ungeandert. Die effektiven Material- 
gréBen fiihren zwangslaufig mit den iiblichen Formeln zu den Gleichungen (9) und (10). 


8. DieBeziehung zwischen der spannungsoptischen Konstanten’ S und dem Elastizitatsmodul EF. 
Die spannungsoptische Konstante S wurde als die Hauptspannungsdifferenz je Isochromatenord- 
nung bei 1 em Modelldicke — beim Arbeiten mit Natriumlicht — definiert. Mit (9) kann mit 
Hilfe des Hookeschen Gesetzes die Beziehung 56, — @(Ao)2/E aufgestellt werden. Beim Fort- 
schreiten um eine Ordnung kann nach der Definition fiir Ao nun S gesetzt werden, wenn die Ver- 
gleichsstabe immer eine Starke von 1 cm haben, damit Witd 02, =a 97... 

Werden Probestiicke verschiedener Hartungsgrade verglichen, so gilt 

Pine | 25h hy (11) 
E, E, Sy So Ss 


1 FE. Ménch, Ing.-Arch. 16 (1948), S. 269. 


oder 


7 * 
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Gleichung (11) sagt aus, daB beim Arbeiten im elastischen Bereich (Hooke) zwischen dem Elasti- 
zititsmodul und der spannungsoptischen Konstanten bei Stiicken mit verschiedenen Hartungs- 
graden ein Zusammenhang hesteht. s 

9. Der Einfluf8 des Kriechens auf die Konstante x, Bei der folgenden Betrachtung soll % 
eine durch einen Versuch ermittelte Konstante sein, die nicht nur den rein elastischen Anteil, 
sondern auch noch einen Kriechanteil enthalten kann. Der Quotient x/x—=y gibt Auskunft 
daritber, wie gro der Kriecheinflu& oder weitere Zeiteinfliisse im ermittelten %-Wert sind. Ist 

— 1, so wird der Versuch im elastischen Gebiet durchgefiihrt. Ist yp >1, so ist noch ein Zeit- 
einfluB vorhanden. 

Der Einflu®B des Kriechens auf die GréBe von y soll an einigen Beispielen erlautert werden; 
hierzu werden die Werte des ungeharteten Materials aus Tabelle 2 genommen. Dieses Material 
erfordert zur Erhéhung der Ordnung um eine Einheit durch eine Belastungsdehnung bei einem 
lem starken Stab eine Spannung von Ag — 10 kg/em®; nach Gleichung (8a) ist dann és 
= Ax/a= 180/42 - 10! 4,3/104. Um durch Zeiteinfliisse — keine Spannungserhéhung — eine Er- 
héhung der Ordnung um eine Einheit zu erhalten ist eine Dehnung Ae, — 180/b, - 10* erforderlich. 

Wird b, sehr hoch mit 20 eingesetzt, so ist Ae, = 9/101. 


Abb. 5 u. 6. Der Randeffekt bei zunehmender Belastung. 


Bei der folgenden Betrachtung bedeutet %), soviel wie: 1/10 einer Isochromatenordnung 
wird durch eine Belastungsdehnung bewirkt; die restlichen 9/10 miissen, um eine Steigerung der 
Ordnung um eine ganze Einheit zu erhalten, durch Zeitvorgange gedeckt werden. Bei %,, ist 
einzusetzen: Ag = 0,1-10 = 1 kg/em?; Ae = 0,1 Ae, 4- 0,9 Ae, = 8,5 - 10-4. Damit wird 

= Raf Lei koran: 

Hiya 34,5 Ve) em? / Roe em 
Mit steigender Beanspruchung wird die Becinflussung des x-Wertes durch Zeitvorginge immer 
stirker, d. h. die Steigerung der Ordnung wird durch kleinere Spannungsdifferenzen bei gréBeren 
Zeitdehnungen erreicht; die x- und w-Werte steigen an. 


Tabelle 3. Der Einflufs der Zeitdehnungen auf die Gréfe von x und wy. 


X 
| 


= 15,2 | yw, = 1,00 (aus Belastungsdehnung) 


|. = 

(elgg == L Tsk | Wyse == 105 
1X lo, i 25,0 Wor. = 1,64 
%\01 = 34,9 Rig = OP 


Ist bei der Eichung y = 1, so kann ohne Bedenken das Material bis zur Eichspannung im 
Versuch belastet werden; die Werte E und S gelten dann fiir jede Spannung, die die Eichspannung 
nicht iiberschreitet. 

Ist bei der Kichung y >1, so gelten die bei der Eichung ermittelten Werte Eund S nur fiir 
eine Spannungscharakteristik, die der des Eichstabes entspricht, d. h. der Spannungsverlauf in 
Abhangigkeit von der Zeit und Temperatur muB beim Eichstab und Modell der gleiche sein. 


10. Der Randeffekt des Kondensationsprodukts Dekorit, Ein randeffekt- und spannungs- 
freies Material erscheint in der spannungsoptischen Einrichtung in allen drei Blickrichtungen 
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dunkel. Eine Platte mit Randeffekt zeigt verschiedene Isochromatenordnungen von der Platten- 
mitte zum Rand hin. 

Ausgehend von einem Material mit der ersten und zweiten Isochromatenordnung wird zur 
Erreichung des Bildes von Abb. 5 die Last um 65kg gesteigert ; zur Erreichung des Bildes von Abb. 6 
um weitere 25 kg. Bei 65 kg Laststeigerung war eine Linienzunahme von etwas iiber 4, bei wei- 
teren 25 kg um rund 2 zu beobachten. 

Die Auswertungen der Abb. 5 und 6 lassen den Randeffekt als eine Zunahme der effektiven 
Materialstarke zur Plattenaufenschicht mit einer Uberlagerung eines Eigenspannungszustandes 
erkennen. Der Verlauf der effektiven Materialstirke ]aBt sich 
ebenfalls aus den Abbildungen ermitteln. 


11. Beschreibung einiger Bruchbilder, Ist die Belastungs- 
geschwindigkeit sehr hoch, so zeigt der zerrissene Stab keine 
oder nur geringe verbliebene optische Effekte; der Stab er- 
scheint auch nach dem Bruch, wie seine Umgebung, in der 
optischen Apparatur dunkel. Ein mit geringer Belastungs- 
geschwindigkeit zu Bruch gebrachter Stab zeigt in der Apparatur 
noch Zeitdehneinfliisse, deren Gréfe von der Hartung des Ma- 
terials abhangig sind. Der Stab der Abb. 7 ist aus einem Ma- 
terial gearbeitet, das nur leicht gehartet ist; er wird mit ge- 
ringer Geschwindigkeit bis zum Bruch belastet. Sehr gut sind 
die verbliebenen Zeitdehnungsanteile, des vor der Belastung 
vollkommen schwarz erscheinenden Stabes, zuerkennen. Abb. 8 
zeigt einen Zugstab nach der Entlastung; seine Last wurde 
langere Zeit konstant gehalten. Bei einem Celluloidstab, der 
langsam belastet wird, zeigt es sich, daf der Bruch von einem 
Zentrum mit mehr oder weniger groBer Geschwindigkeit aus- 
geht. Es steigt plétzlich an einer bestimmten Stelle im Bruch- 
querschnitt die Isochromatenzahl an; das gleiche Ansteigen der 
Isochromaten wird dann noch an weiteren Stellen, bis der Bruch 
eintritt, beobachtet. 


12, Zusammenfassung. Die Versuchsergebnisse lassen sich 
erklaren, wenn der spannungsoptische Effekt als eine Ordnung 
von optischen Grundelementen gedeutet wird. Die optischen 
Grundelemente besitzen eine gewisse Doppelbrechung, die aber 
durch die vollkommene Unordnung der Elemente im unbelaste- 
ten Material nicht nach auBen in Erscheinung treten kann; 
eine aufgebrachte Belastung ruft eine Ausrichtung hervor und uy. 7. 5. gugstitbe mit verschieden 
liBt die bisher verborgene Doppelbrechung nun nach aufen starken Zeiteinfliissen. 
sichtbar werden. 

Im belasteten Material kénnen drei Dehnungsanteile unterschieden werden, die sich am mar- 
kantesten durch die GriéRe des Quotienten, der aus dem optischen Effekt und der zugehérigen 
Dehnung gebildet wird, unterscheiden. 

Versuche mit verschieden stark geharteten Dekoritmaterialien ergeben, daB durch die Har- 
tung eine Zunahme des effektiven Querschnitts bewirkt wird; der effektive Querschnitt ist der 
Querschnitt des gebundenen Materials; nur dieser ist befahigt, die eingeleitete Last weiter- 
zugeben. 

Fiir den rein elastischen Bereich (Hooke) wird fiir Materialien mit verschiedenen Hartungs- 
graden eine Konstante x gegeben, die auch praktisch bei Versuchen mit erhéhter Modelltempera- 
tur erhalten bleibt. AnschlieBend wird ein Faktor yp definiert, der Auskunft iiber die GréBe der 
rein elastischen Dehnung zur Zeitdehnung bei der Erhéhung der Isochromatenordnung um eine 


Kinheit gibt. 


(Eingegangen am 17. September 1951.) 


Anschrift des Verfassers: Dr.-Ing. A. Schafer, Miinchen, Technische Hochschule, 
Institut fiir Metallurgie und Metallkunde. 
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Uber zwei mit der Tragfliigeltheorie 
in Zusammenhang stehende Integralgleichungen. 


Von J. Dorr. 


1, Aufgabenstellung, Bei der Untersuchung der Strémung durch ein senkrechtes Fliigel- 
gitter! wurde die Integralgleichung 
+1 


fl=s¢ g(y) Cig = —”) ay mit h=h’+ ith’ (1) 


a1 
behandelt, die im Sonderfall h’’=0 den Zusammenhang zwischen Abwind und Wirbelstarke 
an einem Fliigel des senkrechten Fliigelgitters beschreibt. Fiir die Sonderfalle h’’—0 und 
h’=0, —1>h”’>-+1 wurde unmittelbar gezeigt, da® die Integralgleichung (1) den Ausdruck 
eee aed le 9 5 ay LESS 
kGin 7 Coj 5 f(y)|/ Sin? — Sin 5 
s()= — aie ee = aay 2) 
n \/ Gin’ - —Gine= kh / Gin? | —Gin? Cite peas 
=i 
als Lésung besitzt. Die in (1) und (2) auftretenden Integrale sind dabei wie itblich als Cauchysche 
Hauptwerte zu verstehen. Man kann zeigen, wie in der zitierten Arbeit erwahnt wurde, da (2) 
auch dann noch die Lésung von (1) darstellt, wenn h=h’+-ih”, h’+ 0, h’’+- 0 ist. Der Beweis 
dafiir wurde nicht gebracht, da er fiir die technisch interessierenden Lésungen (h’”=0 und 
gegebenenfalls h’=0, —1>h’”>--1) bedeutungslos ist. Der Beweis fiithrt indessen auf bemer- 
kenswerte mathematische Ergebnisse und soll deshalb an dieser Stelle nachgeholt werden. 


1 
chs 


2, Beweis der Lésung fiir komplexes h, In der zitierten Arbeit wird gezeigt, daB man mit 


eed OTT. eae Se: 
Ci = OW Eee Olle (3) 
die Integralgleichung (1) auf die folgenden beiden zuriickfiithren kann: 
+1 a5 
gee eh pera re 2() 1 [{ yi) 
a(t) =a) len iL eee (78 pale) = 3% | Fes dg. (4) 
aay Ri: 


Hier kennzeichnet der Zeiger 1 bzw. 2 jeweils den ungeraden bzw. geraden Funktionsteil. Unter 
der Voraussetzung h’’—=0 oder h’=0, —1>h’’>-+1 bildet sich durch die Transformation (3) 
das Stiick der reellen Achse von —1 bis +1 der y-Ebene auf das Stiick der reellen Achse von 
—1] bis +-1 der -Ebene ab und zwar eindeutig. Deshalb sind die Integralgleichungen (4) in diesen 
Sonderfallen identisch mit dem Typ, fiir den Betz in seiner Dissertation den lésenden Kern an- 
gegeben hat. 

Im allgemeinen Fall, wenn h komplex ist, bildet sich das Stiick der reellen Achse von —1 bis 
-+1 der y-Ebene nicht mehr auf die reelle Achse der £-Ebene ab. Das Bild dieser Strecke in der 
§-Ebene ist vielmehr eine von —1 iiber 0 nach +1 durch die komplexe Ebene laufende doppelt- 
punktfreie Kurve endlicher Linge, die iiberall eine Tangente hat und antisymmetrisch zur reellen 
Achse ist. 

Auch bei der Integration iiber diesen allgemeineren Weg ist der Cauchysche Hauptwert 
zu nehmen, was ohne weiteres méglich ist, da der Integrationsweg in der &-Ebene iiberall eine 
Tangente hat. 

Wir betrachten nun den Ausdruck 


10(m) = = | i Ae. (5) 


s 


Der Integrationsweg s soll dabei eine von —1 nach +1 beliebig durch die &-Ebene gefiihrte 
doppelpunktfreie Kurve endlicher Linge sein, die iiberall eine Tangente hat. Die y- Werte 
sollen immer auf dem Wege s liegen, und das Integral ist als Cauchyscher Hauptwert zu ver- 
stehen. Die Lésung der Integralgleichung (5) fiir den beschriebenen allgemeinen Integrations. 


1 J, Darr, Ing.-Arch. 19 (1951), S. 66. 


XX.Band1952. Dérr: Zwei mit der Tragfliigeltheorie in Zusammenhang stehende Integralgleichungen. 89 


weg s findet man am schnellsten, wenn man formal bekannte Uberlegungen iibernimmt. Da die 
so gefundene Lisung spater direkt verifiziert wird, wird der formale Lésungsweg ohne Kritik 
der Voraussetzungen nur kurz skizziert. 

Mit den Transformationen 


y= — cos 0 & = — cos £ 
erhalt man aus (5) die Integralgleichung 
Pee i, geo) sin Bb 
MA ome 3 | cos f —cosa wine (6a) 


iS 
Der Weg S ist hier eine von 6=0 nach {= 2 laufende doppelpunktfreie Kurve endlicher 
Lange, die ebenfalls itberall eine Tangente hat. Das Integral ist wieder als Cauchyscher Haupt- 
wert zu verstehen. Wir machen nun ebenso wie in dem bekannten Sonderfall, da®B S iiber die 
reelle Achse von 0 nach x gefithrt wird, folgenden Lisungsansatz: 
P(f) sin B= Y a, cos nf. (6b) 
n=0 

Man findet durch funktionentheoretische Betrachtungen, daf auch fiir den beschriebenen 

allgemeineren Weg S die Beziehungen 
~ cos np __ _sin no 
| cos B — COS & aa. sin a 


fit Te (alee ea, (6c) 
8 

giiltig bleiben, wenn a auf dem Integrationsweg S liegt und das Integral als Cauchyscher Haupt- 

wert verstanden wird. Daraus folgt noch auf Grund einer einfachen trigonometrischen Relation 


[oer dp so coenn fiir == by Me BY aorwe (6d) 


cos Pp —cos a 


iS 
Geht man mit dem Reihenansatz fiir P(/) sin / in die Integralgleichung (6a), so kann man 
mit (6c) alle Integrationen ausfiithren und erhalt 


I a : 
W(«) = ae ay CesiOe (6e) 


Nun kann man (6b) auf Grund von (6d) auch wie folgt schreiben: 
. ip ees y | 
Phy) i —— | a a, sinnf df . 
m | cos Bf — cos % = 
s 
wenn man die Vertauschbarkeit von Summation und Integration als erlaubt voraussetzt. Fiir 
den letzten Ausdruck kann man aber auf Grund von (6e) auch schreiben 


ee 5, 


FG) SDE oe x J cos B —cos a 
Ss 


Die hiermit gefundene Lésung der Integralgleichung (6a) ist nicht vollstandig, da auf Grund 
von (6c) fiir n—0 auch 


Pla) etnias ee re 
% mu J} cos fh —cosa 
iS 
eine Lisung ist, wobei die Konstante k beliebig gewahlt werden kann. Durch Ubergang auf die 
Koordinaten 7 und & erhalt man als formale Lésung der Integralgleichung (5) den Ausdruck 


k 1 fw /1—# 
Bee amie fae ro () 


Das Integral ist als Cauchyscher Hauptwert zu verstehen. Es diirfte schwierig sein, gestiitzt 
auf diese Herleitung allgemeinere Voraussetzungen anzugeben, unter denen die Lésung (7) giiltig 
ist. Die Konvergenzvoraussetzungen der Reihen (6b) und (6e) auf den Kurven Sin der komplexen 
Ebene sind namlich wesentlich enger als bei den Fourierreihen fiir Funktionen, die nur fiir eine 
reelle Variable definiert sind. Deshalb wird nun die Richtigkeit der Liésung (7) auf direktem 
Wege unter gewissen Voraussetzungen verifiziert. 
Durch Einsetzen von (7) in (5) folgt 
(4) = —-3 aoe [ oo dy dé (8) 
i= ; - 


= ae 


Ste 


Ss s 
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bei Beachtung der aus (6c) fir n=0 folgenden Beziehung 
d& 
\ (jos) | amen 


8 


si) (Cauchyscher Hauptwert). (9) 


Wegen (9) diirfen wir zum Integranden des zweiten Integrals in (8) einen beliebigen auf s inte- 
grierbaren Summanden hinzuaddieren, der nicht von € abhangt. Wir diirfen z.B. schreiben 


"= — 3 | ea aitee | Oy) pone 


s s 


Das ist aber identisch mit 


w(7)) = Ses i wONVI— oy ap te 


Gee) Ay 7) 


Nach einer kleinen Umformung wird daraus 
Lage 1 yl1—y? 
tat) Slee as) | Vie le =A y= ny ee 


- 10b 
ye JI wO)—w) gage. | OP) 
2 oe = <3 esate 
a ) Vl1—& Vee y= | 
Wir setzen nun voraus, auf dem Wege s seien 
w(y) und 5, (9) stiickweise stetig. (11a) 


Die Unstetigkeitsstellen von w(y) seien y,. Dann ist auf dem Wege s fiir eine beliebig kleines, 
aber festes € 
w(y) — w(n) 
oho 
wenn 7 festgehalten wird und y auf s variiert. Unter der Voraussetzung (11b) darf man die 


stiickweise stetig und stiickweise stetig-differenzierbar fiir |7—y,|> e¢, (11b) 


Reihenfolge der Integrationen des zweiten Doppelintegrals vertauschen. Man beweist das, | 


indem man fiir den Weg s in der £- und y-Ebene eine geeignete Parameterdarstellung wahlt. 

Ein brauchbarer Parameter ist z.B. die von einem beliebigen Wegpunkt aus gemessene Weg- 
lange. Dann kann man einen von L. Schwarz! angegebenen Beweis direkt auf unser Doppel- 
integral iibertragen. Nach Vertauschung der Reigenfolge der Integrationen sieht man auf Grund 
von (9), daB das zweite Doppelintegral in (10b) verschwindet. Es bleibt also nach einer kleinen 
Umformung des ersten Doppelintegrals: 


Mit (6c) findet man 
V1 —y? fe 
- dy == ae dy— ae. 
i io eee ee 


s s 


Also ist 


wn) = 3000) | ae = ln) 


Unter den Voraussetzungen (11) ist somit die Beziehung (8) richtig fiir alle 7-Werte, die auf dem | 


Wege s liegen mit Ausnahme der endlich vielen Werte 7=y,. Damit ist der Beweis fiir die Rich- 
tigkeit der Lésung (7) der Integralgleichung (5) unter den Voraussetzungen (11) erbracht. Die 
Integralgleichung (1) 1aBt sich bei komplexem Parameter h auf die Integralgleichung (5) zuriick- 
fiihren. Deshalb ist die Lésung (2) von (1) auch bei komplexem Parameter h richtig, wenn die 
vorgegebene Funktion f(x) auf dem Wege —1<x<-1 1 stiickweise stetig und stiickweise 
stetig-differenzierbar ist. 

(Eingegangen am 10. Oktober 1951.) 


Anschrift des Verfassers: Dr.-Ing. Johannes Dérr, (16) Darmstadt, Technische Hochschule, Institut fiir 
praktische Mathematik. 


1 L, Schwarz, Luftfahrtforschung 17 (1940), S. 385. 
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Biege- und Torsionsschwingungen von Stiben 
mit beliebigen Querschnitten. 


Von H.-R. Lembcke. 


1, Einleitung. Die ,,elastische Achse“ eines Stabes ist diejenige in Stablangsrichtung ver- 
laufende Linie, in welcher alle duBeren Krafte angreifen miissen, um eine torsionsfreie Biegung 
hervorzurufen. Bei allen stabférmigen Gebilden, deren elastische Achse nicht mit der Schwer- 
achse zusammenfallt, entsteht eine Kopplung zwischen den Biegeschwingungen beziiglich der 
beiden Haupttragheitsachsen und der Torsionsschwingung. Diese Erscheinung ist bekannt und 
im Rahmen der Berechnung von Flugzeugtragfliigelschwingungen teilweise behandelt worden!. 

Samtliche bekannten Untersuchungen sind spezieller Natur und beschranken sich auf die 
angendherte Berechnung der Koppelschwingungszahlen und -formen in Fallen, bei denen sich 
nur zwei Schwingungsarten, eine Biege- und eine Torsionsschwingung, gegenseitig beeinflussen. 
Hierbei liegen die elastische Achse und eine der beiden Haupttragheitsachsen des Querschnittes 
in einer Ebene, so daB die Biegeschwingungen um die andere Haupttragheitsachse fast ungestért 
bleiben. Der Querschnitt ist dann meistens nur einfach symmetrisch mit einer Symmetrieachse, 
die mit der zuerst genannten Haupttragheitsachse identisch ist. 

Zwecks Unterscheidung wird die Kopplung zwischen einer der beiden Biegeschwingungen 
und der Torsionsschwingung mit Zweierkopplung und diejenige zwischen beiden Biege- 
schwingungen und der Torsionsschwingung mit Dreierkopplung bezeichnet. 

Die vorliegende Arbeit? befaBt sich mit der Erweiterung der teils bekannten Theorie der 
Zweierkopplung auf den allgemeineren Fall der Dreierkopplung bei beliebigen unverwundenen, 
stabférmigen Gebilden. 


2. Herleitung der Bewegungsgleichungen und des Variationsproblems, Die nachfolgenden 
Ausfiithrungen gelten unter Zugrundelegung der ,,linearisierten‘‘ Elastizitatstheorie. Samtliche 
Dampfungseinfliisse, die Querkraftverformungen und die Rotationstragheit der Stabelemente 
beziiglich der z- und y-Achse werden vernachlassigt. Die Massen- und Tragheitsmomenten- 
verteilung sei stetig. 

Aus Griinden der Anschaulichkeit wird als Beispiel fiir Stabe mit beliebigen, in der EKin- 
leitung gekennzeichneten Querschnitten ein einseitig eingespanntes, ungleichschenkeliges Winkel- 
eisen betrachtet, welches Biege- und Torsionsschwingungen ausfiihrt (Abb.1). Wird mit U die 
potentielle und mit T die kinetische Energie des schwingenden Stabes bezeichnet, so verlangt 
das Hamiltonsche Prinzip, daB das Zeitintegral 


ee 


stationar wird. Wahrend eines beliebigen Schwingungszustandes betragt die potentielle Energie 


1 2 2 22 (x, t)\? Aa(x, t) \? 
Un = + { [Er (Se) + erm MS) + ore (HEY) ae. 
0 


Die Rotation der Stabelemente erfolgt um die elastische Achse. Diese ist identisch mit der 
Verbindungslinie der Querkraftmittelpunkte aller Stabquerschnitte. Mit den beiden Kom- 
ponenten der Schwerpunktgeschwindigkeit 


(ie ——— ee 4- s(x) da Ee =e (x, ))) - | : 
PEG 5 felon Seo) | - 


1 H. Wittmeyer, Zur Berechnung der Eigenschwingungszahlen und -formen yon Fliigeln. Jb. dtsch. 
Luftf. 1 (1939), 5S. 505. ; 

2 Gekiirzte Fassung einer von der Technischen Hochschule Hannover genehmigten Dissertation 1949 
(D 89); Hauptberichter: Prof. Dr.-Ing. O. Flachsbart; Mitberichter: Prof. Dr.-Ing. habil. A. Pfliiger. 
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ergibt sich fiir die gesamte kinetische Energie 


1 
Tal ie(ap (22 air es Cleat S26 1) : 


0 


L O,(x) oo ce 


Im Bereich kleiner Auslenkungen verschwindet die GréBe « (x,t) in den beiden trigonometrischen 
Gliedern; denn das Produkt «(x,t) da(x,t)/dt fiihrt bei einem harmonischen Ansatz a= dy singt 
auf 1/20? m sin 2@t und darf somit als Ausdruck zweiter Ordnung vernachlassigt werden. Beide 
Energieformen, in das Hamiltonsche Integral eingesetzt, ergeben 


Gala, A[a (x, t) sin (B —oa (x, t))]\? 
-( ee s(x) Eat one oe 5 


ta 


, 


el J {ee (e®) [love (os #))® $2 8%) eet) te (24) $- (2 (a t))® — 2 8% (x) 20 (4) be (258)] 


wars? Og (2) (ot (2,1)? — EJ? (x) (Yun (50)? — E I(x) (2x2 (%0))? — | 
— GT (x) (a, (x, t))?\ dxdt= stat. 


(2.2) 


ausgelenkte Lage 


Q-Querkratimitte/gunkt \ 
1S — Schwerpunkt 


Abb. la. Eingespannter Stab. Abb. 1b. Schwingender Stab. 


Hierin bedeuten die Fufzeiger x und ¢ partielle Ableitungen nach x und t; ferner ist 
s*(%) = s(x)cos fh, s¥(x)—s(x)sinf und Op (x)= O,(x) + w(x) s?(x). 

Das Variationsintegral (2,2) enthalt drei Koordinatenfunktionen y (x,t), 2(x,t) und a(x,t), 
die einzeln so beschaffen sein miissen, daB die Forderung J = stationadr gewahrleistet ist. Die 
zugehorigen Eulerschen Differentialgleichungen! lauten 


[EB J* (2%) Yu (% t) Jew Fo (%) Ye, t) 4+ u(x) 87 (x) On (x, t) = 0, (2,3) 
[Ei J¥(%) 2x (%, t) lex + 6 (%) 24 (%, t) — u(x) 8%(%) &u(x, t) = 0, (2,4) 
[G P(x) a. (%, t)}e — Og (2) ou (%, t) — p(x) 8* (x) Yul %, t) + u(x) 89(x) zu (x,t) = 0. (2,5) 


Da Dampfungseinfliisse unberiicksichtigt bleiben, schwingen alle Teilchen des Stabes synchion 
und das Hamiltonsche Integral kann nach Einfithrung der Ansitze fiir die harmonische Schwin- 
gung nach der Zeit t ausgewertet werden. 


Mit der dimensionslosen Langskoordinate € = = , den Amplitudenfunktionen Y(é), Zt) 


und a*(&) und der Kreisfrequenz @ lauten die Ansatze 
y(€,t) = Y(E)sin wt, 2(€,t) = Z(€)sinwt, a(&,t1) = a*(&) sin wt. 


1 R. Courant und D. Hilbert, Methoden der mathem., IPSs Mh Sp LOL, D Auflage. Berlin 1931, 
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Nach ittblichem Rechnungsgang ergibt sich das Variationsintegral 
1 
Tay | fo We) (VO) +28 EYE a* + 2G —26G Zu) + Oo) a4] — 
oe B+(&) (Y”’(é))? — By(€) (2 (€))? — D(&) (a*’(€))*} dé = stat. (2,6) 
mit den entsprechenden Eulerschen Differentialgleichungen 

[Be (6) Y"(8)]” — wu(€)¥ (8) — wu (E) s* (8) o*(E) = 0, (2,7) 
[Bx (E) 2” (E)\"” — wu (E) Z(E) -+ ww (€) s7(€) a* (6) = 0, (2,8) 
[D(E) a*’(E)\’ + @?O9(E) a* (E) + w?u(E) s*(£) ¥(E) — w? w(&) s(E)Z(E)= 0. (2,9) 


Hierin bedeuten die GréBen B() = EJy(&)/l4, B*(€) = EJ*(&)/l! reduzierte Biegesteifig- 
keiten beziiglich der zur y- bzw. 2-Achse parallelen Haupttrigheitsachsen, D(&)— GT (é)/I? 
reduzierte Torsionssteifigkeit und Striche Ableitungen nach &. 


3. Lésung des Variationsproblems mit Hilfe des Ritzschen Verfahrens, Stimmt die elastische 
Achse mit der Schwerachse iiberein, so verschwinden im Variationsintegral und in den Differen- 
tialgleichungen die sog. Kopplungsglieder. Die Schwingungen der drei Freiheitsgrade y, z, a 
werden unabhangig voneinander und kénnen fiir sich berechnet werden. Sie heifSen ,,ent- 
koppelte* oder ,,freie** Biege- bzw. Torsionseigenschwingungen und die zugehorigen Eigen- 
funktionen ,,freie‘* Eigenfunktionen. Entsprechendes gilt fiir die Eigenfrequenzen. 


Es liegt die Vermutung nahe, daf diese ,,.entkoppelten Schwingungsformen besonders 
gute Naherungsansatze fiir die Ritzsche Rechnung darstellen, da die Kopplungseffekte nur eine 
mehr oder weniger starke Stérung der freien Schwingungen bedeuten. Deshalb werden an Stelle 
der vorderhand unbekannten Lisungsfunktionen der Variationsaufgabe fiir alle drei Funktionen 
Y(&), Z(&), a*(€) Linearentwicklungen nach freien HKigenfunktionen der ersten, zweiten, drit- 
ten usw. Ordnung angesetzt. Bereitet im Falle veranderlicher Massenbelegung sowie verander- 
licher Biege- und Torsionssteifigkeiten eine exakte Berechnung der freien Schwingungsformen 
zu groBe Schwierigkeiten, so sind rechnerische oder graphische Naherungsverfahren heran- 
zuziehen!}. 

Die Naherungsansatze enthalten frei verfiigbare Parameter, die durch die Rechnung so 
bestimmt werden, da die Forderung (2,6) erfiillt wird. Hierdurch entsteht ein lineares homo- 
genes Gleichungssystem, dessen Koeffizientendeterminante verschwinden muf. Diese Bedingung 
liefert obere Schranken fiir die exakten Eigenwerte. Zwecks Unterscheidung werden die Nahe- 
rungswerte iiberstrichen. 


a) Allgemeiner Fall des stetig veranderlichen und beliebigen Stabquer- 
schnittes. Fir die Amplitudenfunktionen Y(&), Z(€) und a*(&) werden zweiparametrige Nahe- 
rungsfunktionen angesetzt. Diese brauchen vorerst weder orthogonalisiert noch normiert zu 
sein, sondern miissen nach H. Wittmeyer! lediglich die wesentlichen Randbedingungen erfiillen: 


VY=>a,¥,%6aY,, 2=),24,--),2,, a =cak-+ 0% - (21) 


Mit diesen Ansidtzen werden die einzelnen Anteile des Variationsintegrals (2,6) berechnet. Aus 


den Bedingungen 
dS (a,b,c) _ 9 aJ (a, b, ¢) ie ON ase dJ (a, b,c Die 


da, E¢ 1 ae 0b, Ob» 
0 J (a, Dy @) dJ (a, b, ¢ b,c — () 
Bag 0, nae , 


entstehen sechs lineare homogene Gleichungen mit der Koeffizientendeterminanten-Gleichung I. 


1 H. Wittmeyer, Zur Berechnung der Eigenschwingungszahlen und -formen von Fligeln. Jb. dtsch. 


Luftf. 1 (1939), $.505; L. Collatz, ‘Eigenwertprobleme und ihre numerische Behandlung, 2. Auflage, 
Leipzig 1949. 
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Koeffizientendeterminante I, 


1 il 
wo? | w¥?dé w?[wY, Y,dé a aes 
fr : {! 0 0 -+- CaN +o Mee 
—[B(Y/pdé —[BYTY"dé - Yyatdé  Yyatdé 
0 


0 


1 1 
@ PES ELE o Javad ‘ P o{us* oP | js? 


A 5 “ OF oe 0 
—[BrY{yidé —/B(Y,)dé Yaoi dé * Yoo,d6 
0 0 
1 1 
w?|uZ2dé w?| uZ,Z,dé ys ie see 
eS i cc ee 
0 0 
—fBy(Z’)2d— —[BYZ{Zidé Zyo05 dé Zio; de 
A= 0 0 
if 1 
w2{ uZ,Z,dé w*| uZed& he bee 
0 0 {i Ba di 5 oe in a 
—| BYZ{ Zi dé — | BY (Z))?d&é Line de -Z, 03 dE 
0 f 
1 J 
1 x 1 ee 1 a 1 @ (6) *2 qd me} (Q) *. *d 
+ ous? + of us — | ws? — of us af Cor g Gea g 
0 0 0 0 
-Y, a*dé -Y,a*d& Lyon dé -Ziatd— —{D(af)*d& —f{Dof’atdé 
0 0 
1 1 
1 ar L i ae Pete ‘a2 #27 £ 
+ w? | us* = oo” | s — @*{ us? — w*| us as [Ooo Sam | % se: 
0 0 0 0 1 1 
‘Yyatdé "-Yyo8dE —’ -Z,abdf = Zab —JDat'at'dé —fD(at)*dé 
0 0 


Bei 7/-, J- und ¢-parametrigen Funktionsansatzen ergibt sich die Determinanten-Gleichung II 
von (y+ 0-+ C)-facher Ordnung. Sie liefert 7--0-+€ reelle und positive Eigenwerte und kann 
als allgemeine Liésung des vorliegenden Schwingungsproblems angesehen werden. Bei ent- 
sprechenden Ansiatzen lassen sich alle Kigenwerte und Eigenfunktionen mit beliebiger Genauig- 
keit ermitteln, wobei allerdings mit steigender Ordnung der Rechenaufwand rasch anwichst. 
Die Bedeutung der einzelnen Glieder ist durch Vergleich mit den Elementen der Determinante I 
fir 7=0=C=2 leicht zu erkennen. 


b) Sonderfall der prismatischen Stabe mit konstantem, aber beliebigem 
Querschnitt. 1. Herleitung der freien Eigenfunktionen fiir die Biege- und 
Torsionsschwingungen (kurzer AbriB der Theorie). Die Biegeschwingungen genannter 
Stabe werden bekanntlich durch die partielle Differentialgleichung vierter Ordnung 


jy eG ia CN 5 


Oxt CAT (3,2) 
und die Torsionsschwingungen durch die partielle Differentialgleichung zweiter Ordnung 
07x (x, t) 07x (x, t) 
GT - A rere 0, — Aaa 0 (3,3) 
beschrieben. Mit 
fe : Pay ap x)2 @, 
f=7. (y= SEE, yp ee (3.4) 
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Koef fizientendeterminante IT. 


0 OS mee sen CRN ree SONNETS OTe Paste 
Vi Wh -—— OB, ee, ae Bi Me ere rm pie 
| = | | oe | | chN 

| 
| | < | | | XS | | | x | 
| | el | | mt | | SGT 
Toews 8 0 2 0S 0 UT Pe aH 
0 i Vi, Wi, SOP Lala) Take = he 
0 Ds A Na Pee rarer ge! IR OI Sy Nea gg en ie 
| (ae | | bo | | eri | 
| | \ | | | X | | | NS | 
| | Soe ass | | he | | Ney 
0 ee eg = 0 SE ee 
ee eel Pa eT PA RO Re ee 
ee eee ee PS = i PR = = ee 
SS | | \ | | | X | 
| | x | | | NS | | | ‘ | 
| | | | Kl | t al 
ee eee ee Le A, HPS 


und den Zeigern n bzw. m zur Kennzeichnung der laufenden Ordnung lauten die allgemeinen 
Lésungen in einer fiir die Koeffizientenermittlung giinstigen Form 


¥o(E)n = Cy (cos kp & + Coj kr &) + C,(cos kr E — Cof kh) + C, (sin kpé + Gin kZé) | (3,5) 
+ C,(sin hE — Sin kX), | 
Ope) sin kee + D, cos | (3.6) 
Die Koeffizienten ergeben sich durch Anpassung dieser allgemeinen Kigenfunktionen an die 
Randbedingungen. Durch eine zweckentsprechende Wahl des Bezugskoeffizienten werden die 


Eigenfunktionen normiert, so da die einzelnen Integrale der Determinante I einfache Werte 
annehmen. Wegen der Orthogonalitat ist 


[Y,(é)¥i(€) dé = 0, [Yn (E)¥7(E) dé = 0 fitrae= i, 


0 0 


1 1 
fm (E) on (£) dé = 0 , [om (&) on (§) dé = 0 fiir m ~k 
0 0 
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und nach S. Timoshenko! weiterhin 


é 6 (3,7) 
[(a@) 48 = sh 5 [(H%)® (a (E))® + (a8 (1))2] 


Alle Integrale der Gestalt 
1 
[rate Ey be eS {[Zn(&) om (E dé =e nm 
0 


haben die Bedeutung von Kopplungsfaktoren zwischen Biege- und Torsionsschwingungen. 


Fiir die ersten drei untersuchten Lagerungsfialle (s. Tabelle 1) zeigen die Eigenfunktionen 
der Biege- und Torsionsschwingungen gleichen Aufbau. Folglich kann hierfiir eine gemeinsame 
Gleichung der Kopplungsfaktoren angeschrieben werden: 

[ * F i " = k* cos k¥ sin in + k* sin k¥ cos k* 

Me — [” = pr a 
n(&) Om(E) dé B (ayy — (43) 
) 


vee vos EATEN y a Ch ne Hx! y Mugen Nee a ne 
k* Gj kX sin kX — k* Sin k> cos Kirn ) ki cos k* cos ky + k¥ sink? sin k¥ — kh” 


Te | [a y\2 / 2 (3,8) 
(KY + (Ke) (2)? — (Ke) | 


= ke cos k” Goj kX + kX sin ki, Sin k* Aa ke 


fir ky A ke. 
(er)? + (kay 


Fiir den vierten Einspannungsfall ,,unterstiitzt—unterstiitzt gilt 


— k? cos es cos Pee a: sin k* sin ke tas k> 
USe= ik Y: =e? fir ky ~ ke 


und 


fir ki =k, 


a 
Ho. 
uw 
S 
\o 
ee 


Die Herleitung dieser Gleichungen ist trivial und kann im Original der Arbeit nachgelesen | 
werden. 


Da die Kigenwerte k, und k;, der beiden freien Biegeschwingungen gleich sind, gilt [%,, =I™ m- 
Kine Unterscheidung durch die beiden Zeiger y, z ist nicht unbedingt notwendig, aber wegen 
der guten Ubersicht zweckmaBig. 


Die in Tabelle 1 angegebenen numerischen Werte der Kopplungsfaktoren umfassen nur die 
Kombinationen zwischen Biege- und Torsionsschwingungen der ersten beiden Ordnungen. Wird 
fiir spezielle Untersuchungen eine gréBere Genauigkeit des Ritgschen Verfahrens gewiinscht, 
so lassen sich ohne gréBere Mithe auch die Kopplungsfaktoren héherer Schwingungsordnungen 
aus (3,8) bestimmen. In den meisten Fallen jedoch werden die Angaben der Tabelle 1 ausreichen. 


lS. Timoshenko, Schwingungsprobleme der Technik, S. 250. Berlin 1932. 
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2. Dreierkopplung. Fir die Amplitudenfunktion Z(&) gelten die analogen Aussagen wie 
fiir Y(E). Bei konstantem Querschnitt, unveranderlicher Massenbelegung und den in der 
Tabelle 1 angegebenen Lésungen der Integrale (3,7) ergibt sich aus der Determinante I nach 
einigen Umformungen die Determinante III. Durch geeignete Spalten- und Zeilenadditionen 


Koeffizientendeterminante ITI. 


1 
— y y 
cos” B B(k})s 0 0 0 + Ih + IF 
~ @u cos? B 
i 
0 cos® p BF (KY) 0 0 +Dn +DPn 
~ @2u cos? B 
if 
- 2 Zz Zz 
0 0 sin® p BY (ki)! 0 = 195 ah £18 
Oy sin® B 
1 
0 0 0 sin® B BY(ke)s Peis —I22 
~ o2wsin2B 
OQ 
Ih +P ae = ae ; 
= 2 w* py s® 
@ 
; ¥: z z 2 
+14, + Ia = 10% — Ws. 0 BIE D (kz)? 


gelingt eine iibersichtliche Darstellung der Determinanten-Gleichung, die sich fiir die Zahlen- 
rechnung besonders gut eignet: 


o\ pa | Ghar aa Wea Cia )e 
Oi =A B.D ee vig 6 a C D 
Prete de be APG URL 
pp — Sp pF aa ae) (. | Gam 
Py Ty, Pu Tn 4 Thy Tye Talal _ 9 
Aaa ee CEL RD Ay 
Hierin bedeuten 
met BY(k})" 1 _ Bey) Bape en AY 
Wcos2 so OA icose im cos? B= @2u cos? B ” “sin? B  @2usin? B’ 
EE) ee Oo DE tre 
sin™p * | @°jfsm p* 2 js? DCOe Ss ae 2s? IH? ws? 


mit Im =m und ky—=k,. Die Gleichung (3,10) gilt fiir die ersten drei Lagerungsfalle 
,,eingespannt—frei“, ,,eingespannt—eingespannt“, ,,eingespannt—unterstiitzt‘*, wahrend fiir den 
vierten und letzten Fall ,,unterstiitzt—unterstiitzt‘‘ nur die halben Werte der Glieder A bis D 
in Rechnung zu setzen sind. Ihre sechs Nullstellen, die am einfachsten durch Probieren und 
Interpolieren gefunden werden, bilden obere Schranken fiir die unteren sechs Kopplungseigen- 
frequenzen. 
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Die Ermittlung der Ritz-Parameter erfolgt nach bekannter Methode. Mit z.B. a, als Bezugs- 
gréBe ergibt sich fiir die restlichen Parameter 


On, 
== bi Sb c — c 


Cor? ‘L_——) 2 
ON ON ON ON ON 


Die Nennerdeterminante #y entsteht aus der Frequenzdeterminante III durch Fortfall der ersten 
Spalte und ersten Zeile. Aus ihr ergeben sich die Zahlerdeterminanten Dies Ds, Orient) ee 
wenn die entsprechende Spalte der Nennerdeterminante gegen die mit a, multiplizierte fort- 


'gefallene erste Spalte (ohne das erste Glied) vertauscht wird. Nach elementaren Zwischen- 
rechnungen wird 


Oy = BCD cos?f sintf s! if ST ray a E why Tiny" a | 


(3,11) 
LEME oa ay ay 
B G Dae 
a, = a, cos*f sin'f st CD [EI Ty, + Fy, Ty], (3,12) 
Bi, = — a, costf sin? st BD] E(Ty,)* + PUY)? — 25 (ry Py, — 71, Pay, (3,13) 
Oy, = — a, cos*f sin®f st BC [ETy, Fy, + FIy, Ty), (3,14) 
OQ. . | / Je ; eA Lee 
eseercoeyeem 6s) BCD) 1%, r— a oe @.15) 
apg aaa els ci 
= r,| 7 21 4 a c “ | : 
- : { an tit) eae. 
}., = a, cos*f sintf s? BC D\— Ty, | E— ~2 - = 
2 1 B | 1 B C D (3,16) 
Tol LE es PI. 
ry,| F 22 | % al 4 D | ; | 


3. Zweierkopplung. Bei speziellen Querschnittsformen, wo die Verbindungsgerade vom 
Schwerpunkt und Querkraftmittelpunkt mit einer der beiden Haupttragheitsachsen zusammen- 
fallt, ist eine exakte Lésung der Differentialgleichungen (2,3), (2,4), (2,5) méglich. Hierbei 
findet eine ,,Entkopplung*: der Biegeschwingungen in Richtung dieser einen Achse statt. Ist 
z.B. die z-Achse eine Haupttragheitsachse und identisch mit der oben definierten Verbindungs- 
geraden, so verschwindet der Winkel (, und die Biegeschwingung in z-Richtung wird entkoppelt. 
Um einen Vergleich der Naherungsrechnung mit der exakten Rechnung zu ermoglichen, wird 
fiir diesen Fall die Frequenzdeterminante IV angegeben, die sich aus der Determinante III bei 
z.B. vier-parametrigem Ansatz von freien EKigenfunktionen entwickelt. 


Koeffizientendeterminante IV. 


Zweierkopplung, vier-parametriger Ansatz von freien Eigenfunktionen. 


ey 4p Geen By ee OT 1.0, = Co O* 
ia 0 zoo Fiesty 

oO? wu ‘ 

0 Eells meee 410 
ca Pte —0 
Oe fo) D (ki)? . 

Q (ky) 
iter +P ition 
Og _ Dik? 


. Ph A Ps 2 2a ise 20% ws? 

4, Untersuchung der Differentialgleichungen (2,3), (2,4) und (2,5) bei Dreierkopplung und 
Lésung derselben fiir den Sonderfall der Zweierkopplung. Fiir Stabe mit der Linge nach ver- 
anderlichen, beliebigen Querschnitten ist es im allgemeinen nicht mdoglich, eine strenge Lisung 

8 
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in geschlossener Form zu finden. Auch bei unveranderlichen Querschnitten beliebiger Form er- 
scheint eine Integration praktisch undurchfithrbar. Nur bei Zweierkopplung gelingt die Losung. 
Die Bedeutung dieser exakten Rechnung liegt weniger in ihrer Anwendung fiir die Praxis (hier- 
fiir ist eine numerische Durchrechnung zu zeitraubend), als vielmehr in der Méglichkeit, an einem 
Sonderproblem die Genauigkeit der Ritzschen Naherungen zu diskutieren. 

a) Dreierkopplung £40. Ist der Stabquerschnitt der Linge nach konstant, so lauten die 
drei gekoppelten Differentialgleichungen (2,3), (2,4) und (2,5) 


ES? Yxxxx(%st) + ¥u(%,t) + ws cos Boy(x,t) = 0, 
EJ? Bexex(%t) + [b24(%,t) — wssin Boy(x,t)—0, (4,1 
GT dxx (x,t) — Ogdu(x,t) — ws cosh yu(%,t) + jes sin Bze(x,t) == 9 


Durch geeignetes Eliminieren und Differentieren entstehen hieraus die drei entkoppelten Diffe- 
rentialgleichungen 


010 z (x, t) OIA 62519) 08 z (x, t) 0°z (x, t) 0° z(x, t) OP z(x, 2) } 
Aisrage 1 Aa Gah get OB paso epee ge ore Se Sao | 
810 y (x, 1) a0 (x, t) A y(x,t) | 4 OP y(x,t) a y (x, t) ay(x.t) «| 
iS Bale A, Ox O82 A, “ax® Ot2 fF A, ox* ott A; Ox? att A, ars =, (4,2) 
a OC, t) 019n (x, t) a(x, t) | Oa (x,t) | a(x, t) | a(x, t) 0 
core 2 Ox? o# A A MR ae 2) MELE G2 Se ah 


Sie sind linear, homogen und von zehnter Ordnung mit den gemeinsamen konstanten Koeffi- 
zienten 


ys aEFES CT, = —ZESED Oo, A, =~ OT (EF EJ»), 


A= 8° [EJ? cost + EJ*sin?6]— “ (EX LEI), Ap=S6T, Ae ope 


Die Integration der Gleichungen scheitert an der geschlossenen Formulierung der Wurzeln der 
charakteristischen Gleichung fiinften Grades. 


b) Zweierkopplung f=—0. Die z-Achse sei identisch mit der Hauptachse in z-Richtung. 
Der Winkel 6 verschwindet, und die Stérung der Biegeschwingungen in z-Richtung durch die 
Schwingungen beziiglich der beiden anderen Freiheitsgrade entfallt. Es besteht demnach nur 
noch eine Kopplung zwischen der Drehbewegung und der Translationsbewegung in y-Richtung. 
Die Differentialgleichungen zehnter Ordnung der Funktionen y(x,t) und a(x,t) reduzieren sich 
auf solche sechster Ordnung, deren Behandlung méglich ist. 

Mit sin S—=0 und cos6—1 wird aus den Gleichungen (2,7) und (2,9) 


B  Y' (6) — op ¥ (6) — ot ws a*(€)=0, (4.3) 
Doa*"' (€) + w? Oga* (€) + w? us Y(E)=0. (4,4) 


Die beiden zugehérigen ungekoppelten Gleichungen sind lineare homogene Differentialgleichungen 
sechster Ordnung mit konstanten Koeffizienten. 


YE (E) 4 C8" yar (gy — ne YG) te ea, (4,5) 
oXVI(é) + tne oXIV (£) — am a"(é) ba ee a*(E)— 0. (4,6) 


Beide Gleichungen ergeben dieselben EKigenschwingungszahlen. Wegen der einfacheren Rand- 
bedingungen wird Gleichung (4,5) weiter untersucht. 
Fiir den meistens zutreffenden Fall 


o*u MiSs, w* Oo\3 
spB (9 zus) > (Sp) 
lautet ihre allgemeine Lisung in trigonometrischer Gestalt 


Y(é) = A Oink, E+ Boj k,F + Csin k,é + D cos k,&+ Esin kgE-+ Feosk,&. (4,7) 
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Hierin bedeuten 


Oa 2 Ih Upeiees w? Og\3 
peli ae ver 4, BD Re 3H") —-(a) 
ND es IE COS) eee sae TT Aaa Vay CNS eg 7 > 
e G97" @ ae 
3D) + 3B 
ha ee a eae Tie 
PNY /(@? @o\2 | wo? nee Oo 
m= 2/(S5° ei eae ake OE YR 
rs : jw? Oo\2 wo ‘ a o? Og 
k= | A (3D Be ree ee) locates 
Mae 9 (w2 Oo\2 ho wo? Og 
= |/ 2 V (pe) + BE - cos (60°49) + 282 


Als Sonderfall wird der einseitig eingespannte und anderseitig frei schwingende Stab aus der 
Vielzahl der Einspannungsméglichkeiten herausgegriffen. 
Die Randbedingungen heifen 


Y(0)=0, Y(0)=0, Y"(1)=0, Y(1I)=0, o*(0)=0, oF (1) =O. 


Die vorletzte Bedingung ist identisch mit Y!’(£)—0 und die letzte mit Y’(1) — sa YL) ==0 

Uo 
entsprechend den Forderungen der Differentialgleichung (4,3). Mit diesen sechs Randbedingungen 
und den einzelnen Ableitungsfunktionen entstehen sechs lineare homogene Gleichungen, deren 
Koeffizientendeterminante V verschwinden mu. 


Abb. 2. Gleichschenkeliges Winkeleisen L 3660 100 x 100 14, 
einseitig eingespannt, anderseitig frei (Zweierkopplung). 


Die Querschnittsdaten und Materialkonstanten werden mit 
6 stelliger Genauigkeit angenommen. 


s = 3,50000[ cm] 
= 2,10.000: 10-4 [kg -'s? - cm* >] 


GT = 186 340 - 10? [kg - cm?] 


EJ* = 744000 - 10 [kg - cm?] 
99 = 0, + ps? = 6,35000 10-2 [kg « s?}. 


os 
Koeffizientendeterminante V. 
0 1 0 il 0 1 
key 0 ke 0 ks 0 
ki Sin k, k2Coj ky —kesink, —kcos k, —kgsin k, —kjcos kg "3 


e k? Gof ky kB Gin k, —k3cosk, —kgsin k, —kjcosk, W—k3sin k, 
0 ka 0 ka 0 ks 


kg, Co) ky kyg,Gin k, kg,cosk, —h,g,sink, k3g,cosk, —k,g,sink, 


* kt "ke B* ks 4 
mit a=—1—— g=1— =, oat |e io far chktse kee he 

5. Beispiel fiir eine Zweierkopplung. Um einen Uberblick iiber die Fehler des Ritzschen 
Verfahrens zu erlangen, werden die Eigenschwingungszahlen eines gleichschenkeligen Winkel- 
eisens (Abb. 2) nach der Differentialgleichungs- und der Ritzschen Methode durchgerechnet. Die 
Ergebnisse (Tabelle2) zeigen, daB schon bei einem vier-parametrigen Ansatz die Fehler des Ritz- 
schen Verfahrens praktisch verschwinden. Deshalb sind auch bei Dreierkopplung und gleich- 


artigen Ansatzen nur Abweichungen dieser GréSenordnung zu erwarten. 
8* 


102 


Lembcke: Biege- und Torsionsschwingungen von Staben usw. Ingenieur-Archiv 


Tabelle 2. Eigenschwingungszahlen eines gleichschenkeligen Winkeleisens (Abb. 2). 


Bezeichnung Oo; [s— 2] Ory [s—2] ry fs] 
Genaue Frequenzen nach Determinante V. . 48,948 269,10 336,80 
Niherungsfrequenzen nach Determinante IV . 48,963 269,30 337,47 
Fehler des Ritzschen Verfahrens + 0,030% + 0,074% + 0,20% 
Freie Eigenfrequenzen nach den Beziehungen Ol = Op oo 
(3,4) und Tabelle 1 49,388 301,38 309,51 
EinfluB der Kopplung —(,89% |—10,7% + 8.8% 


Abb. 3. Querschnitt des experimentell untersuchten Stabes. 


b, = 5,05. mm], 
b, = 4,85[ mm], 
a, = 52,00[ mm}, 
a, = 30,90) mm], 
Ta 8,00 mm], 
DD = 3,25 [mm], 
s = 16,0 [mm], 
c, = 8,05[ mm], 
co = 17,5 [mm], 
SOLES 
Qu—e23,27% 
Querschnittsflache: qi = 4,260[ cm?], 
Ex sig? 
Stabmasse je Langeneinheit: pe = 3,410 - 10-5 [=| 
cm 
Rotationstragheitsmoment, 

bezogen auf die elastische Achse: Og = 2,131: 10-4[ kg - s?] 

bezogen auf die Schwerachse: O, = 1,258 - 10-4[ kg - s?] 
Axiales Tragheitsmoment bezogen auf die 

Hauptachsen parallel zur J* = 2,536[ cm‘), 

z- baw. y-Achse: i j= 13,20[ em‘), 
Torsionssteifigkeit nachF.Schleichert: T = 0,4050[ em4], 
Elastizitatsmodul: E = 2,100 - 10°[ kg/cm?] 
Gleitmodul: G =0,8080+108[ kg/cm?] 


1 F. Schleicher: Handbuch fiir Bauingenieure. S. 173, 
Berlin 1943. 


6. Beispiel fiir eine Dreierkopplung. Zur Bestatigung der vorliegenden Theorie werden die 


gekoppelten Biege- und Torsionsschwingungen eines ungleichschenkeligen Winkeleisens (Abb. 3) 
als Beispiel fiir den allgemeinen Fall der Dreierkopplung bei prismatischen Staben mit willkiir- 
lichem Querschnitt experimentell und rechnerisch untersucht. 


Abb. 4. Der Versuchsstand. 


a) Experimentelle Untersuchung. Abb. 4 
und 5 zeigen die Versuchseinrichtung. Durch einen 
Wechselstrommagneten mit harmonisch pulsieren- 
dem Feld wird der Stab am freien Ende zu er- 
zwungenen Schwingungen angeregt, die im Reso- 


Abb. 5. Die Anordnung des Erreger-Magneten. 


nanzfall angenahert in stationare Higenschwingungen iibergehen. Der Magnetisierungsstrom 
wird in einem Leonardaggregat erzeugt, das in einfachster Weise der Forderung nach weit- 
gehender und feinstufiger Drehzahlregelung gerecht wird. Bei Resonanz stehen Kigenfrequenz 
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und. Erregerfrequenz im Verhaltnis ganzer Zahlen, wenn Dampfungseinfliisse keine Rolle spielen. 
Folglich lat sich aus der Drehzahl des Leonardgenerators die Eigenschwingungszahl berechnen. 


Zur Darstellung der Schwingungsfiguren des Stabendes dienen vier polierte Stahlkugeln, die 
an den beiden Schenkelenden, im Schwerpunkt und im Querkraftmittelpunkt des Stirnquer- 
schnittes aufgelitet sind (Abb. 5). Die Stahlkugeln wirken als konvexe Spiegel und lassen das 
Bild einer nahen Lichtquelle als Reflexpunkt erscheinen. Erfolgt die Bewegung der Kugeln 
schnell genug, so erzeugen diese Punkte im zeitlichen Ablauf helle, scharfe Kurven, die nahezu 
fehlerfrei die wahren Schwingungsfiguren der gewahlten Querschnittspunkte wiedergeben. 
Durch Vergleich aller vier Bilder kénnen dann Riickschliisse itber die Formen der Kigenfunktionen 
am Stabende gezogen werden. Die Mefergebnisse sind in Tabelle 3 unter ,,gemessen‘ angegeben. 


b) Rechnerische Untersuchung. Durch Aufsuchen der Nullstellen der Frequenz- 
determinante II bzw. der Gleichung (3,10) ergeben sich die in Tabelle 3 unter ,,gerechnet* ver- 
zeichneten gekoppelten Eigenfrequenzen. 


Die Ermittlung der Parameter fiir die angenaherten Eigenschwingungsfunktionen ist zwar 
langwierig, bereitet prinzipiell aber keine Schwierigkeiten. Deshalb wird die Betrachtung des Be- 
wegungszustandes nur bei der Stablinge 1—207,0 [em] und @yy = 377,3 [s] durchgefithrt. 


Tabelle 3. Gerechnete und gemessene Eigenschwingungszahlen eines ungleichschenkeligen Winkeleisens 
(Dreierkopplung) sowie Angaben iiber den Kopplungseinfluf. 


Stablinge 1 cm 130,5 Teel | 207,0 | 244,7 
Meore. Kigentreguenz:. °) . 2°. 5. eT a ks ow, Sau 81,59 
Gekoppelte SCKCC Le tne ene Oy » | ~ 81,34 | 
Rigentrequenz \ypemessen. . (92... .. Wy + 80,97 | nicht gemessen 
Abweichung der Rechnung......... (ewe e025 
opplanyseinilnbs <0. 2 onl se I ys eevor iia 0.8 
EeciesPisentrequenzimiiena sare sea a. 6 wi Set 186,15 101,07 73,98 
Gekoppelte [EINES So 5 oe BoB oe Oy | » 183,5 100,32 713,58 
Hicentrequenz: \beemessenm . 9. . «2 « ..5 ; Ory 27 182,5 99,77 73,09 
Abweichung der Rechnung ........ 27> | +0,55 | + 0,6 st OL7 
op planes cutis aaeme ie eee err. % | —2,0 Bea ee Ke, 
Er elLem bce nie quenzunm rms kicurseiemcrren =) oy sol 511,3 277,6 203,2 145,4 
Gekoppelte Bere Cun tawwe-e tn New w os Myyy be 491,7 OBS 201,2 144,5 
igentrequemz: eemessen) = 9... «4 = - Ory si 491,0 DONS 201.0 144,2 
Abweichung der Rechnung ........ % = =I +0,3 | +0,1 + 0,2 
KopplungsermtluG meme emo ees sc) oe.) | % aA Zot eo ee |, Tete =). 
Preie-Nigentrequen7 G40; 6 6 se se ots st | 613,9 452,4 alee nas 
Gekoppelte ({ERRRTA MN ofS ob pmee LO 4c Ory > 630,5 448 8 gue SN 
Eigenfrequenz \gemessen ......... Ory 2 630,0 448 ,2 tae en 
Abweichung der Rechnung ........ vo SUS |) Sr eG =i Be oa 
KopplungseinfluB Roe, a er cpucl Levis xe + 2,5 as 8x, ? 
| 
Rrevestioentrequenzman mye) sate ies wo Se ee: | Dae pe oes 
Gekoppelte ETCCHINE UE ort say ese 3 Oy 1 fare ibe: ree Beas 
Eigenfrequenz |gemessen ........- Oy 2 107 ee i a = Ne 
Abweichung der Rechnung. ......% . vo i a ai ms af Ne a A 
KopplungseinfluB : : 70 NRK ae ? ; 
Mit B= 1,2854, C= 2,1476, D=— 1,1391, E= 0,46238 und F — — 5,6031 wird nach 

den Gleichungen (3,11) bis (3,16) 

= a qe Freman 

Dy —=+3,2251 [emt], d,, = + 14066 a, [em’), , = — 2,9926 a, [em], 

~ a 9. corn 

Oyo = + 11,7628 a,[cm4], 2, = + 3,9957 a, [em], Y., = + 0,10537 a, [em*] 
und mithin a,= + 0,43614a,, b,= — 0,92790 a,, 6,—= + 0,54660 a,, 


c, = + 1,2389 a, [em], ¢,= 0,032 671 a, lems ||. 


F : ; Sr genre 
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Hierdurch ist die Eigenschwingungsform vermége der drei Koordinatenfunktionen y(é, t), 2(&,t) 
und a(&,t) angenahert beschrieben: 


¥(E, thry = a, {+ Ay (sinky & — Gin ky €) + cosky § — Gof ky €+ | 
+ 0,43614 [+ By (sinky & — Ginky &) + cosky § — Wo) ky] sin @ryt, | 

2(E, ty = a {— 0,92790 [-+ BF (sink: £ — Sinks &) + coshz & — Gof ky €] 
1 0,54660 [+ fe (sin ke & — Sink £) + cosht & — Cofkzé]} sinoryt, | 

até, try = a, {+ 1,2389 sin k* € +- 0,032671 sin kh? é} sin @yyt [cm] (6,3) 


fiir a,in [em]. Die Konstanten (3, 6), und kx, k; sind der Tabelle 1 zuentnehmen. Mit Y(1 pee ae” 
(Pluszeichen fiir gerade und Minuszeichen fiir ungerade Ordnungen) nach der Beziehung (3,7) 


(6,1) 


(6.2) 


S 2 
g \of [(YE)n)? dé 
~ 0 1 
q/sec SY Ow HON = AY (1); ed 
600 © aoe 4 
N wt? SS d 
SS =. lauten die auf a, bezogenen 
500 x 7] I, : : 
<= i =a Maximalverschiebungen 
oe T a | des Stabendes 
oO ee es Stabende 
at Sa Sk i 
Le BY: ey Bo Y (l)ivy = — 11277 a,, 
ui ee See 
zn —— <a Z (Ir = + 2,9526 ay 
z Sa = = a (1)rv= + 1,2062 [cm Jay. 
100 = 
i = . . 
‘2 Der Faktor a, wird jetzt 
_ 
0 B05 199 2070 arom” ~~ so gewahlt, daB 
Abb. 6. Verlauf der freien und gekoppelten Eigenfrequenzen in Abhingigkeit 
von der Stablinge. y Y (12 Z(1)? 


mit der wirklichen Amplitude der reinen Transversalbewegung iibereinstimmt (Punkt | vgl. 


Abb. 8 und Abb. 9). Diese wird der Abb. 8 entnommen und betragt ~ 0,285 [em]. 


Somit wird 


0,285 
SS pe cm |= 0,094 [em 
1 Vi1a772 + 2.95262 al Lens 
und 
YQ) = —0,106 [em], ZA) == 0277 [ems oy = Od (6,4) 


Abb. 9 zeigt diese Verschiebungen. 


Abb. 7. Stabende im Ruhezustand von oben gesehen, Abb. 8. Stab im vierten Resonanzzustand. 


c) Vergleich der gerechneten mit den gemessenen Schwingungszahlen und 
Ermittlung des Kopplungseinflusses. In Tabelle 3 sind die Rechnungs- und MeBergeb- 
nisse zusammengestellt und in Abb. 6 kurvenmafig aufgetragen. Die mittleren Toleranzen bei 
der Angaben infolge MeBfehler und ungenauer Erfassung der Querschnittsabmessungen sind ge- 
ringer als + 0,5 % und damit durchaus tragbar. In allen Fallen liegen die gerechneten Werte 
etwas itber den gemessenen. Diese geringen Differenzen entstehen erstens durch den Verfahrens- 
fehler der Ritzschen Methode, zweitens durch die behinderte Querschnittsverwélbung im Ein- 
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spannbereich, die bei den kurzen Stablangen die Verdrehungssteifigkeit etwas erhéhen wiirde, 
und drittens, allerdings nur sehr wenig, durch die vernachlissigte Dampfung. 

Trotz dieser Unterschiede wird der KopplungseinfluB, d.h. der prozentuale Unterschied 
zwischen der Kopplungseigenfrequenz und der freien Eigenfrequenz nach Gleichung (3,4) mit be- 
friedigender Genauigkeit erfafBt. 


d) Vergleich der berechneten Schwingungsfiguren des Stabendes mit den 
wirklichen. Abb. 7 zeigt das obere Stabende von oben gesehen im Ruhestand. Die hellen 
Reflexpunkte auf den vier Kugeln sind deutlich zu erkennen. 
Die schwarzen Striche markieren die Trégheitshauptachsen. 
Der Kiirze halber wird nur der vierte Resonanzzustand mit 
@iy = 375,2 [s-!] wiedergegeben (Abb. 8). Die zugehérige 
freie Bewegungsform ist die erste Torsionsschwingung mit 
der Schwerachse als Drehachse. Die Abweichung durch 
den Kopplungseffekt ist unverkennbar. Abb. 9 zeigt die 
durch die Angaben (6,4) festgelegten Figuren fiir den obigen 
Zustand. Sie stimmen gut mit den wirklichen in Abb. 8 
iiberein und lassen deutlich erkennen, wie sich die Bewegung 
aus einer reinen Translation des Querkraftmittelpunktes 
und einer Rotation um ihn zusammensetzt. pe ee necro MIG ees ence 


des vierten Resonanzzustandes. 


7. Zusammenfassung. Es wurden die gekoppelten Eigenschwingungen von Staben mit belie- 
bigen Querschnitten untersucht. Bei ,,Zweierkopplung‘: (Definition siehe Einleitung) wurde 
eine exakte und bei ,,Dreierkopplung“‘ eine angenaherte Berechnungsmethode der gekoppelten 
Schwingungszahlen und Schwingungsformen hergeleitet. An einem ungleichschenkeligen Winkel- 
eisen als Beispiel fiir eine Dreierkopplung wurden die theoretischen Ergebnisse experimentell 
bestatigt. Der Verlauf der Eigenfrequenzen in Abhangigkeit von der Stablange 1aBt folgendes 
erkennen (Abb. 6): Die gegenseitige Stérung der Biege- und Torsionsschwingungen ist so gerich- 
tet, daB sich benachbarte freie Biege- und Torsionseigenfrequenzen auseinanderdriicken, und 
zwar um so starker, je naher sie zusammenliegen. Diese Erscheinung ist von den gekoppelten 
Pendelschwingungen her bekannt. Die GréBenordnung der Stérung durch die Kopplung (beim 
ersten Beispiel maximal 10% und beim zweiten maximal 8 %) zeigt, da® sie in manchen Fallen 
Beachtung verdient. Von groBem EjnfluB ist der Abstand zwischen elastischer Achse und Schwer- 
achse, da er quadratisch in die Rechnung eingeht. Bei kurzen, verdrehungsweichen Staben ist 
ein starkerer KopplungseinfluB zu erwarten, da hierbei die erste freie Torsionsfrequenz im Be- 
reich der unteren freien Biegefrequenzen liegt. Bei langen, verdrehungssteifen Staben dagegen 
unterscheiden sich die unteren, gekoppelten Biegeschwingungszahlen kaum von den jeweiligen 
freien Eigenschwingungszahlen, denn die erste Torsionsfrequenz liegt in diesen Fallen wesentlich 


hoher. 


(Eingegangen am 13. Oktober 1951.) 


Anschrift des Verfassers: 
Dr.-Ing. H.-R. Lembcke, Westerrénfeld bei Rendsburg, z. Zt. Kiel-Friedrichsort, MAK. 
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Uber die Stabilitat punktweise ausgesteifter Rechteckplatten. 
Von S, Woinowsky- Krieger. 


1, Einleitung, Bekanntlich lift sich die Beulsicherheit dinner, einem Druck oder Schub 
ausgesetzter Platten erheblich steigern, wenn man ihre Flache durch eine oder auch mehrere biege- 
steife Rippen unterteilt. Unter Umstanden kann es sich als zweckmaBig erweisen, die Platte nur 
punktweise festzuhalten und es erhebt sich dann die Frage nach der Wirksamkeit einer solchen 
Aussteifung. 

Die Platte 1aBt sich in einem innerhalb ihrer Begrenzung vorgegebenen Punkt entweder ge- 
lenkig oder auch in Verbindung mit einer Einspannung abstiitzen. Bei der ersten Anordnung 
leistet die Stiitze nur der Durchbiegung der Platte Widerstand, bei der zweiten zugleich ihrer Ver- 
drehung, sei es in einem vorgezeichneten oder in einem beliebigen Sinn. Die Abstiitzung und die 
Einspannung kénnen ebensogut starr als auch nachgiebig, etwa elastisch, sein. lm letzten Fall 
ist der vom Stiittzpunkt ausgeiibte Widerstand der Durchbiegung 
der Platte respektive ihrer Verdrehung proportional. 

Die folgende Untersuchung beziecht sich auf eine Rechteckplatte, 
die einer einfachen Art von Druckbelastung ausgesetzt ist, den 
Navierschen Randbedingungen geniigt und nur in ihrer Mitte 

punktweise abgestiitzt ist. Die angewendete Methode ist in- 
4 4 r-—g  dessen allgemein genug, um auch in wesentlich anders gearteten 
| Z Belastungs- und Stiittzungsfallen zum Ziele zu fiihren. 


2, Aufstellung der Beulbedingung, Die in Abb.1 dargestellte, 
ringsum gelenkig gelagerte und in ihrem Mittelpunkt unverschieb- 
saan | | lich gestiitze Platte sei einem allseitigen Druck von der GréBe p 
ae ae. je Langeneinheit der Randlinie ausgesetzt. Die Platte mége inihrem 


y é ; 
Mittelpunkt als eingespannt angenommen werden, doch gelten 
Abb. 1. Die Rechteckplatte ris is 6 5 pS 
ceeiccondiaacsar tt ca alle folgenden Uberlegungen auch fiir die gelenkige Stiitzung, sofern 


man nur von einer Beulflache mit doppelter Symmetrie ausgeht. 
Die Grundlage der Stabilitatsrechnung bildet die bekannte Differentialgleichung der gleich- 
formig zusammengedriickten Platte 


AAw + 4 Aw=0, (1) 


worin w die Auslenkung, D die als konstant angenommene Biegesteifigkeit der Platte und 
A = 0?/0x?+-0?/dy? den Laplaceschen Operator bezeichnet. Mit den Abkiirzungen 


P It LLacr 
2 Tp On = ey ? Br=Vor—x? (2) 


a9 


1a8t sich eine passende Lisung der Differentialgleichung fiir das Gebiet 0 < vex ay in folgender 
Form anschreiben: : 


w= = (A Ginony + B, Cof any + C, Sin By + D, Col Bry) sin on x . (3) 


; 1 
Im Gebiet 0> y > — 30 soll derselbe Ausdruck gelten, nur mit —y an Stelle von y. Da die 


elastischeFlache (3) die Navierschen Rundbedingungen w= 0, Aw=0 langs der Rander x= 0, xa 
wie ersichtlich erfiillt, bleibt es nur noch iibrig die Parameter 4,, B,, ... den Grenzbedingungen am 


J : oe : 
Rande y= 5 sowie den Ubergangsbedingungen langs der Mittellinie y=0 anzupassen. Indem 


man am Rande speziell die Navierschen Bedingungen einfiihrt und die Symmetrie in der Mitte 
beachtet, d. h. von den drei Bedingungen 


(= 0, 4) ao 8) oo (a 


Nie 
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Gebrauch macht, hat man die Méglichkeit, drei Integrationskonstanten, etwa A,, B, und D,, 
durch die vierte auszudriicken. Hiermit ]aBt sich der Ausdruck (3) auf die Form 


Rie 2G Bo( Sin On Y — Tq Coj on) —an{ Gi Buy — Tg! Gof buy) anes (5) 


2 


bringen, die fiir jedes n nur noch eine Integrationskonstante enthalt. Die elastische Fliche (5) 
liefert lings der Mittellinie der Platte die Querkraft 


: a 
qy= — D ls (Avo) = — Dx? DS! CB, sin a, %. (6) 
of v0 n=1 
DieSumme (6) muG nun identisch sein mit dem analytischen Ausdruck fiir die halbe Reaktion, 
die der Stiitzpunkt in der Mitte an die Plattenhalfte y>0 abgibt. Bezeichnet P die gesamte Reak- 
tion der Stiitze, so folgt aus der eben angestellten Uberlegung die Gleichung 


ce oo 
: Po.sonn. 
oe ee Ci Bs sin &, X= — » aoe SID OX, X, (7) 
n= 


n= 


aus der sich fiir == 1,3, 5>0.- 


n+1 
EPS 
SS mera (8) 


ergibt, wahrend fiir gerades n der Koeffizient verschwindet. Die Substitution von C, in die 
Reihe (6) liefert nun den Ausdruck fiir die verbeulte Plattenflache mit nur noch einem unbestimm- 


ten Faktor P. Im Folgenden wird allein der Biegepfeil in Plattenmitte x= ; a, y= Obendtigt und 


dieser ergibt sich auf die eben beschriebene Weise zu 


(9) 


pi On Bn 

Ware die Stitzung elastisch nachgiebig, so hatte man P= fw, setzen miissen, sofern f eine Fe- 
derkonstante bedeutet. Durch Substitution dieses Ausdruckes in die Gleichung (9) kame man 
dann auf die endgiiltige Bedingung zur Bestimmung des Beulwertes x?= p/D. Im eingangs an- 
genommenen Fall einer unnachgiebigen Stiitzung lautet diese letzte Bedingung einfach w—0. 
Dies hat zur Folge entweder P=0 bei einer unverdinderten Gleichgewichtslage der Platte oder 
aber das Verschwinden des Summenausdruckes in (9). Nach Beriicksichtigung der Abkiirzungen 
(2) und nach Einfithrung der Bezeichnung 


SS 
Dax Rly osree 


72s pe (10) 


De 
womit sich also k= ax/z ergibt, wird man schlieBlich auf die Gleichung 
bn 7 —— 
S s Qa Ig & Vn? — i) 
n=1,3,8,... re /n? — k? 
in k? gefiuhrt, deren kleinste Wurzel den gesuchten Beuldruck p, liefert. 

Es ist bisher stillschweigend von der Annahme n?>k? ausgegangen worden. In der weiter 
unten durchgefiibrten Auswertung der Beulgleichung ergab sich indessen 1<k<3, wonach das 
erste Summenglied in der Gleichung (11) einer besonderen Behandlung bedarf. Es ist nun, leicht 
nachzuweisen, da® der Sonderfall n?<k? durch einfaches Vertauschen von n und k bei gleich- 
zeitigem Ersatz der zugehérigen hyperbolischen Tangens-Funktionen durch die trigonometrische 
erledigt werden kann. Die Auswertung der Beulbedingung macht geringe Miihe, da fiir n>3 
simtliche Tg-Funktionen durch eine eins ersetzt werden diirfen; fiir n>11 wurden die Reihen- 
glieder, und zwar durch Entwicklung nach den Potenzen der Grobe k/n, auBerdem noch auf die 


Form 


= 0 (11) 


I 1 hk? (12) 


n /n2— kh? eo 2 n3 


gebracht, wonach der betreffende Rest der Summe sich zu 0,00086 k? ergeben hat. 
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3. Ergebnisse der Stabilititsberechnung. Die aus der Beulbedingung (11) resultierenden 
kleinsten Beulwerte k2 sind in der ersten Zeile der nachstehenden Tabelle enthalten. Sie bestim- 
men den tatsichlichenBeuldruck, falls die konstruktive Ausbildung des Stiitzpunktes auBer der Be- 
dingung w= 0 noch die Erfiillung der beiden weiteren Bedingungen dw/dx—= 0w/dy= 0 in Platten- 
mitte sichert. Ist die Stiitzung gelenkig und nur die erste dieser drei Bedingungen erfiillt, so mul 
mit dem Ausbeulen der Platte in zwei Halbwellen um eine Knotenlinie gerechnet werden, deren 
Verlauf bei b/a>1 durch x= 0 bestimmt ist. Die zugehdérigen Beulwerte sind in der zweiten Zeile 
der Tabelle zusammengestellt, wahrend die dritte Zeile die Beulwerte einer Rechteckplatte ohne 
jegliche Zwischenaussteifung angibt. 


Py a? 


Tabelle der Beulwerte k? = — ‘ 
22D 


Beulwerte 


Art der Abstiitzung Beulbedingung 


im Mittelpunkt US UES 5 | ug 2.0 | BU 0 
a | a a a 
Gelenklose Gleichung (11) 7,38 4,04 2,62 1,65 
Gelenkige k?=1-+-4a?/b? 5,00 Mesihe 2,00 1,44 
Keine k?=1-+-a?/b? 2,00 | 1,44 25 Jodlit 


Nach der Tabelle erweist sich die gelenklose Abstiitzung der Plattenmitte weit wirksamer als 
die Anordnung einer Querrippe, wahrend die aussteifende Wirkung der gelenkigen Punktstiitzung 
der Wirkung einer Querrippe gleichkommt. Bei langlichen Rechteckplatten ist eine Stiitze in 
der Mitte fast ebensowenig effektiv wie eine Querrippe. In solchen Fallen ware die Anordnung 
einer Langsrippe oder aber die Abstiitzung der Platte in mindestens zwei Punkten besser am 
Platze. 

4, SchluBbemerkungen, Das hier dargelegte Verfahren bedarf nur geringer Modifikation, 
um auch in andersgearteten Fallen von Belastung und Stiitzung anwendbar zu bleiben. Grund- 
sitzliche Schwierigkeiten sind eigentlich nur dann zu erwarten, wenn der Ansatz von Levy auch 
bei der Platte ohne Zwischenstiitzung versagt. 

Ist z. B. die Rechteckplatte einem immer noch gleichférmigen, aber in den Richtungen der 
beiden Achsen verschiedenen Druck ausgesetzt, so andert sich wohl die Differentialgleichung, 
nicht aber wesentlich die weitere Herleitung der Stabilitatsbedingung. Sind zwei gegeniiber- 
liegende Plattenrander gelenkig gelagert, so lassen sich beliebige Randbedingungen langs der 
beiden anderen Rander durch entsprechende Abanderung der Gleichungen (4) sofort befriedigen. 
Wird dabei die Symmetrie, etwa durch die auBermittige Lage der Zwischenstiitzung, hinfallig, 
so vermehrt sich mit der Anzahl] der zu erfiillenden Rand- und Ubergangsbedingungen ent- 
sprechend auch die Anzahl der verfiigbaren Parameter. Der Bedingung der Einspannung an 
einer Zwischenstiitze kann bei Asymmetrie der Beulflache nur durch die Einfithrung von min- 
destens einem reaktiven Einzelmoment geniigt werden. Man wird dabei, und ebenso bei mehr 
als einem Stiitzpunkt, auf ein System homogener Gleichungen fiir die unbekannten Reaktionen 
gefiihrt. Durch Nullsetzung der Nennerdeterminante dieses Gleichungssystems ergibt sich dann 
die gesuchte Beulbedingung. Man wird allerdings gezwungen, Einzelmomente an einer Strecke 
oder einer Flache von endlichen Abmessungen angreifen zu lassen; punktférmig konzentrierte 
Kinzelmomente fithren in der Tat auf eine ins unendliche wachsende Plattenneigung im 
Angriffspunkt, die zur Bestimmung der GréBe des erzeugenden Momentes ungeeignet ist: 


(Eingegangen am 15. Oktober 1951.) 


Anschrift des Verfassers: Prof. Dr.-Ing. S. Woinowsky-Krieger, Quebec (Canada) 200 Av. Marguerite- 
Bourgeoys Apt. 8. 
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Uber einige Versuche an Stromungen in Sand. 


Von K, Wieghardt. 


1, Einleitung, Im Folgenden werden einige einfache Versuche mit Sand beschrieben, die 
1947 im Max Planck-Institut fiir Strémungsforschung in Gottingen durchgefiihrt wurden. Zu- 
nachst war ganz allgemein beabsichtigt, Mischungen von Sand und Wasser — von trockenem 
Sand bis zu klarem Wasser — vom rheologischen Standpunkt aus zu untersuchen. Aber sehr 
bald erwies sich schon der trockene Sand als ein Strémungsmedium, das bereits geniigend 
kompliziert fiir eine Sonderuntersuchung war. 

Wahrend in der Erdmechanik das Hauptinteresse in der Statik liegt, handelt es sich hier um 
dynamische Versuche. Die ,,Theorie“‘ dazu besteht im wesentlichen nur aus dimensionsanalytisch 
gewonnenen Formeln, deren Hauptmerkmale jedoch von den Versuchen gut bestatigt werden, 
obwohl eine sicherlich wichtige Eigenschaft des Sandes, namlich die Veranderlichkeit des Poren- 
volumens, die sogenannte Dilatanz!, nicht beriicksichtigt wurde. 


2, Das Versuchsmaterial, Um sicher zu stellen, daf der Sand méglichst gleichférmig war, 
ohne ihn selber noch einmal sieben zu miissen, wurde er von einer Sandpapierfabrik bezogen und 
bestand daher aus gemahlenem Flint. Die Kérner der beiden benutzten Sorten (feiner und 
grober Sand) erwiesen sich unter dem Mikroskop als scharfkantige Prismen. Im Gegensatz 
dazu sind bei natiirlichem Sand, wie FluB- und Seesand, die Kérner nahezu zu Kugeln abge- 
schliffen. Daraus erklart sich, daB die unten angegebenen Materialkonstanten teilweise ab- 
weichen von denen fitr gewohnlich benutzten Sand. Die Porositat z.B., das ist das Verhaltnis 
des Porenvolumen zwischen den Kérnern zum Gesamtvolumen, kann bei Kugelpackungen nur 
zwischen 0,26, und 0,47, variieren; bei der dichtesten Packung bilden dabei die Mittelpunkte 
von je vier Kugeln ein Tetraeder, bei der lockersten Packung liegen die Mittelpunkte von je 
acht Kugeln auf den Ecken eines Wiirfels. Bei nadelférmigen Elementen ist dagegen rein geo- 
metrisch beinahe jede Porositat zwischen 0 und | méglich. 

Die folgende Tabelle enthalt einige Materialkonstanten, alle fiir lose Schiittung. Der Bé- 
schungswinkel wurde dadurch ermittelt, daB eine waagerechte Kreisscheibe in einem mit Sand 
gefiillten Eimer von aufen festgehalten wurde und der Eimer ohne Erschiitterungen langsam 


gesenkt wurde, bis ein kegelférmiger Sandhaufen auf der Scheibe iibrig blieb. 


Feiner Sand Grober Sand 


Lineare Abmessungen . 0,15—0,25 mm| 0,5—1,5 mm 


VE Pemiiise Besm Ge WAC Litters at une ie ie ce) Meee 1,2 g/cm? 1,3 g/em® 
n = Porositat = Poren-/Gesamivolumen, ... . | 0,51, 0-01, 
i) (eae OLeUugi ter aalm errr t uh Gate sls detent oa bs00 1,08 
y = Boschungswinkel iiber Kreisscheibe ... . 36,0° Sooo 
u = tang = Reibungskoeffizient Sand gegen Sand 0,73 0,79 


K.i,p =[(1+ sing)/cos ¢)]? = Koeffizient des 


aktiven ; Erddrucks 2. 
passiven 


0.26 0,24 
309) 4,3 

3. Messungen des Drehmoments von Reibungskriften, Ein mit Sand gefiillter Eimer wurde 
auf einer horizontalen Drehscheibe befestigt, die von einem Motor angetrieben wurde. In den 
Sand wurden zentrisch verschiedene Zylinderstabe und Kegel aus Stahl eingetaucht und das 
Drehmoment gemessen, das der drehende Sand auf sie ausiibte. Dieses Moment hing im Ver- 
suchsbereich von 20 bis 60 Umdrehungen in der Minute nicht von der Drehgeschwindigkeit ab. 
Wegen der relativ gro®en Streuung der Mefergebnisse konnte ferner nicht zwischen den beiden 
Sandsorten (fein und grob) unterschieden werden. 

a) Senkrechter Zylinder zentrisch eingetaucht. Um zu einer formelmaBigen 
Abschatzung der Reibungskrifte zu gelangen, soll angenommen werden, da die in der statischen 
Erdmechanik bis zum Bewegungsbeginn iiblichen Annahmen auch noch fir stationare Be- 
wegungszustande zutreffen. 


if 


1 Q. Reynolds, On the dilatancy of media composed of rigid particles in contact. Phil. Mag. 1885. 
2 K. von Terzaghi, Erdmechanik. 
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Fiir kleine Eindringungstiefen kann der Einflu8 der Eimerwande auf den Druck im Sand 
vernachlassigt werden, so da der vertikale Druck durch das Eigengewicht wird p,——vyy, mit 
y —spezifisches Sandgewicht und y =vertikale Koordinate (positiv nach oben) von der Sand- 
oberflache aus gemessen. Da der Sand sich hier passiv verhalt, ist fiir das Verhaltnis horizon- 
taler/vertikaler Druck der Koeffizient K, maBgebend: p,= K,p»- Die Schubspannung an der 
Zylinderoberflache ist t= 4,p, mit 1, als Reibungskoeffizient Sand-—Stahl, und somit wird das 
Reibungsmoment am ganzen Zylinder 


—H 
d ; 5 
M= [ iy Kpyy ads dy =F y @u,K, 
4) 
0 


mit d als Zylinderdurchmesser und H als Kintauchtiefe. Diese Formel stimmt mit den Versuchs- 
ergebnissen nach Abb. 1 iiberein, wenn man M/: 
Uy K,= 6 setzt, was mit K,~4u, ~1,5 be- cm FE tan® ce (cos oc+ sino) 
deuten wiirde. 

b) Senkrechter Kegel zentrisch ein- 
getaucht. Mit 2a als Offnungswinkel des 6 
Kegels und H als Eintauchtiefe wird der Halb- 


Sand: 


grob fen 0 
° x 65° 


messer des Kegels in der beliebigen Tiefe — y oe es is 
gleich r= (H-+-y) - tg a; das Oberflachenelement fo) x YBp° 


in dieser Tiefe ist eine Ringflache 2ar dy/cos a. 
Senkrecht darauf wirkt der Druck 


pn cos a+ p, sin a, 


; 


Ba Mahe 
I ay g2 
200-42" 
cmt Sand: s js 
grob fein a re 
©) x 71cm © 
7000 | ey ° 
800 $f 


200 


100 ee 
7 7A g CN GINTO, cm 20 7 2 Y 6 8 70 ctm2 
Abb. 1. Reibungsmoment an einem zentrisch Abb.2. Reibungsmoment an einem zentrisch 
eingetauchten Zylinder. eingetauchten Kegel. 


so daB mit denselben Annahmen iiber den Druck wie oben das Reibungselement am Kegel wird 
aan 


M=— — [ty (K, cos a+ sin a) y y 


0 


DY BBE 
r cos a dy =Fy tg?a (K, cosa + sin a)u,H*. (2) 


cos & 


Auch diese Forme] wird nach Abb. 2 bestatigt fiir Offnungswinkel 2~—33° bis 90° und K, ~4 
Fiir 4, witrde sich wieder 1,5 ergeben. eke 

Kin so hoher Reibungskoeffizient erscheint allerdings ziemlich unwahrscheinlich, so da® viel- 
leicht eher ein héherer Wert fiir K, bei Bewegung alsin Ruhe anzunehmen ist. Jedoch geben die 
Formeln die Abhangigkeit von der Eintauchtiefe, dem Zylinderdurchmesser bzw. dem Offnungs- 
winkel der Kegel richtig wieder. Zur Bestimmung der Zahlenfaktoren waren weitere Versuche 
etwa mit normalem, rundkérnigen Sand wiinschenswert. 


4, Bewegungswiderstand von Zylinderstaiben in Sand, Wahrend in den obigen Versuchen 
der Sand in Ruhe war relativ zum Eimer aufer in einer sehr diinnen Grenzschicht um den 
Widerstandskérper herum, wurde eine wirkliche Sandstrémung erst erzeugt durch exzentri- 
sches Kintauchen eines Zylinderstabes. Die Photographien Abb. 3 bis 5 illustrieren die Wirkung 
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auf die Oberflache des Sandbetts. Auf den ersten Blick erinnert die Strémung an die in einer 
sehr zahen Fliissigkeit, etwa bei einer Reynoldszahl von der GréBenordnung 1. Vor dem Stab 


ist eine Erhéhung und ein deutlich ausgeprigter Staupunkt (Abb. 3). Die Stromlinien sind 
auch im Nachlauf (Abb. 5) deutlich er- 


kennbar und glatt und stofen unter einem 
endlichen Winkel zusammen, — ein An- 
zeichen, da} dort auch in der dritten Ko- 
ordinate, der Tiefe, eine Strémung vor- 
handen ist. Kin wesentlicher Unterschied 
gegen die Strémung in einer normalzahen 
Fliissigkeit besteht darin, da®B das Tot- 
wasser scharf begrenzt ist. Danach ist 
der Sand, ahnlich wie gewisse nicht- 
Newtonische Fliissigkeiten ein Medium 
mit einer Mindest-FlieB-Grenze fiir die 
Schubspannung. Jedoch ist hier diese 
Mindestschubspannung proportional dem 
Normaldruck (Coulombsches Reibungs- 
gesetz). Schon in verhaltnismafig kleinem 
Abstand vom Zylinder, wo die Schub- 
spannung kleiner ist als dieser Mindest- 
wert, findet iiberhaupt keine Strémung 
statt. Gerade diese Tatsache erméglicht 
die ganze Versuchsanordnung, namlich 
einen Rundlauf mit sehr kleinem Radius. 
Ebenso wie das Reibungsmoment von 
zentrisch eingetauchten Rotationskérpern a 
erwies sich auch der Widerstand von ex- Boob Bo oeocae bani fiancncns Con de creraang 
(Zylinderdurchmesser 20 mm). 
zentrisch eingetauchten Kérpern unab- 
hangig von der Drehgeschwindigkeit, obwohl das natiirlich nicht zutraf fiir die Verformung der 
Sandoberflache. Da der Druck im Sand nicht isotrop ist wie in einer Fliissigkeit, und auch kein 
Unterdruck im Nachlauf entstehen kann, 
scheint hier der statische Erddruck an der 
Vorderseite des Kérpers allein maBgebend zu 
sein, wahrend der Impuls des bewegten Sandes 
dagegen vernachlassigbar klein ist. Dem ent- 
spricht auch eine sehr kleine Froudesche Zahl: 


Abb. 4. Seitenansicht. Abb. 5. Zusammenfluf hinter dem Zylinder. 


Geschwindigkeit/]/gd (g = Erdbeschleunigung, d— Zylinderdurchmesser) = 0,02 bis 0,3 1m 
Versuchsbereich. . ; 

Es soll deshalb angenommen werden, da der Widerstand proportional dem statischen Erd- 
druck auf die Widerstandsflache 
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W=c,| Kypyyddy= 5 K, yd H?, (3) 
0 
ist mit H als Eintauchtiefe. Nach Abb. 6 geben die Versuche tatsachlich 
W= 25 y d H?, 
so da hier c,K,— 50 wird. 

Als weiteres Beispiel wurde das Widerstandsmoment M, eines untergetauchten, waage- 
rechten Zylinders gemessen; das kleine Reibungsmoment der Haltestange wurde nach Abb. | 
abgeschitzt und abgezogen. Dieselben Annahmen wie oben geben 

1/2 
1 


5 Ma =o, | K, yy drdr=*c,K»ydPH, (4) 


ply 
2 Sie 
0 


mit 1 als Gesamtlange des Zylinders in einer 


met Tiefe H unter der Oberfliche. Nach Abb.7 
ce ist experimentell 
8 


jl 
Kes i 
Mal ya 20H, 
oder c, Ko = 20. 


T 


— 


PIx+00860 


t- 


70 
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Abb. 6. Widerstand eines exzentrisch eingetauchten Abb. 7. Moment eines horizontal eingetauchten 
Zylinders (Grober Sand). Zylinders (Grober Sand). 


In beiden Fallen geben die Formeln die Abhangigkeit vom Zylinderdurchmesser bzw. der 
Lange und der Eintauchtiefe richtig wieder, und nur die Zahlenkonstanten miissen den Messungen 
entnommen werden. In der Statik hingen solche Zahlenfaktoren meist stark von der jeweiligen 
Porositat ab, sie sind verschieden in loser Schiittung oder in festgestampftem Sand. Bei der 
obigen Versuchsanordnung kann das aber offenbar keine Rolle spielen, da ja bei jeder Um- 
drehung immer dieselben Sandmassen aufgewiihlt werden und wieder zur Ruhe kommen (im 
drehenden System), so daB stets der Fall der losen Schiittung vorliegt. 


5, Ausflu8 aus einem Rohr durch eine Lochblende, Ein anderes Beispiel einer Sandstrémung, 
so alt wie die Sanduhr, ist das Ausflie{en von Sand durch ein kleines Loch. Zum Zweck gut 
definierter Versuchsverhaltnisse wurde der Ausflu® aus einem Glasrohr gemessen, das unten 
durch ein diinnes Blech (0,25 mm) abgeschlossen war, in diessen Mitte ein Loch gebohrt war. 
Auch hier war es nicht nétig den Einfiillungszustand des Sandes naher zu beriicksichtigen; die 
Ausflu8menge ergab sich unabhiangig davon, ob man den Sand vorher in das Rohr durch einen 
Trichter hineinrieseln lie (lose Schiittung) oder schichtenweise in das Rohr einstampfte. Denn 
auch hier sind es im wesentlichen immer dieselben Sandmassen, die an der Bewegung teilnehmen 
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und dadurch sehr bald nach dem Beginn des AusflieBens aufgelockert werden, so daB wieder der 
Fall der losen Schiittung vorliegt. 

a) Senkrechtes Rohr. Zunichst sei an die bekannte statische Druckverteilung in 
einem mit Sand gefiillten vertikalen Kreisrohr erinnert. wenn der Boden des Rohres nur durch 
einen so kleinen Gegendruck gehalten wird, daB der Sand gerade anfangt wegzurutschen. Unter 
der vereinfachenden Annahme. da® die Drucke unabhangig von der radialen Koordinate sind, 


g/y' 


70 20 


Q= 0,2 (a-k) Moy 


cm 2 


ine I n B vy ine sh eine Lochblende 
7 o é +} ; shblende Abb. 9. AusfluB von feinem Sand durch eine 
Abb. 8. AusfluB von JU eipe ieee es n ein ochblende (Durchmesser d),. 


da8B wieder p, = K,p, und an der Wand‘ = (1 , pp, lautet die Gleichgewichtsbedingung an einem 
Zylinder der Héhe dh 


+ D? dpp= 7 Dy dh —x Dr dh == D*y dh —x Diy, K, py dh. 


Integration gibt den vertikalen Druck am Boden 


yD a. 
Pr (h= 0) = 3 7 K [1 — exp (— 44, K, H/D)|, (5) 


mit D als Rohrdurchmesser und H als Sandhéhe. Sowie die Sandhéhe zwei oder drei Rohr- 
durchmesser iibersteigt, ist also der Bodendruck unabhangig von der Hohe und zwar gleich 
y D/4 wu, K,. Da der Sand sich hier aktiv verhalt, ist der Koeffizient des aktiven Erddrucks K, 
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einzusetzen. Passiv wiirde sich der Sand verhalten, wenn etwa am oberen Ende des Rohres mit 
einem Kolben ein Druck nach unten ausgeiibt wiirde. 

Auf Grund dieser Silo-Gleichung (5) kann auch fiir den AusfluB von Sand aus einem solchen 
Rohr die Sandhéhe keine Rolle spielen, solange sie nicht zu klein ist. Wenn man diese fiir den 
Hydrodynamiker ungewohnte Tatsache weiB, ist es leicht, dimensionsanalytisch vorauszusagen, 
daB die mittlere Ausflu®geschwindigkeit v im AusfluBloch abhangen wird vom spezifischen 
Gewicht y, der Erdbeschleunigung g und dem Lochdurchmesser d. 

Die Beobachtung zeigt ferner, da8 der Sand nur in einem diinnen, schwach trichterférmigen 
Schlauch iiber dem Loch flieBt, wahrend der itbrige Sand rundherum nicht an der Bewegung 
teimimmt. Infolgedessen kann der Rohrdurchmesser D nur von untergeordneter Bedeutung 
sein, auBer etwa in sehr diinnen Rohren, wo der Schlauch mit flieBendem Sand die Rohrwandung 
schon kurz ttber dem Loch erreicht. Aus y, g und d kann nun nur auf eine Art eine Geschwindig- 


2/8) keit gebildet werden, namlich 


Y 


ee — eee (6) 


Es zeigt sich also, daB das spezifische Ge- 
wicht keine Rolle spielt. Es fallt in der- 


4 selben Weise aus der Rechnung heraus 

wie bei der entsprechenden Ableitung der 
ae AusfluBgeschwindigkeit von Wasser: 

ee 
i Die Ausflu8menge in der Zeiteinheit Q 
; in g/s wird danach 
Q= iy PvV dr. (7) 

4 A 


Da die Versuchsanordnung und die 
Messungen selbst (mit Waage und Stopp- 
02 = \ E uhr) so einfach sind, streuen die Ergebnisse 
| nur so wenig, daB eine Auftragung Q iiber 
| d noch eine kleine Korrektur k nahelegt, 
| _£ so daB Q~(d—k)>” wird. Theoretisch 
0 70 20 30 40 50 60° ‘ 3 

Abb. 10. AusfluB aus einem geneigten Rohr durch eine pera eS (ome ee IEG Ua 

Lochblende (feiner Sand). sache beritcksichtigt, daB der Sand kein 

Kontinuum ist sondern aus endlich groBen 

Kérnern besteht. So kann z. B. itberhaupt kein Sand ausflieBen, wenn das Loch zu klein 

ist: d< k. kist als ein mittlerer Wirkungsradius von der GréSenordnung der linearen Ab- 
messungen der Sandkérner aufzufassen, und Gl. (7) damit abzuandern in 


O=cs7y Ve (d— bP. (8) 


Fiir weite Rohre mit D > rd. 8cm sind die MeBergebnisse nach den Abb. 8 und 9 in cgs- 
Kinheiten fiir den groben Sand Q/y—16 (d—0,11)* und fiir den feinen Sand Q/y—17,6 
(d—0,04)°/*. Daraus folgt fiir den groben Sand c,— 0,65 und k= 0,11 cm, und fiir den feinen 
Sand c,= 0,71 und k= 0,04 cm. 

Bei engeren Rohren macht sich noch ein Einflu8 des Rohrdurchmessers auf c, und k bemerk- 
bar, wie in der Nebenfigur in den Abb, 8 und 9 angegeben. Und zwar wird in einem engeren 
Rohr der Ausflu8 durch dieselbe Lochblende kleiner als in einem weiten Rohr, obwohl der Schlauch 
mit flieBendem Sand die Rohrwandung eher erreicht und die Reibung zwischen Glas und Sand 
wohl kleiner ist als zwischen Sand und Sand. Andrerseits ist dort, wo der Schlauch die Rohrwand 
erreicht der Druck kleiner im engeren Rohr als im weiten. Denn nach (5) ist der Druck, abgesehen 
von sehr kleinen Sandhihen, proportional dem Rohrdurchmesser. 

Kine gewisse Bestatigung der Gl. (8) wurde nachtriglich in einer Arbeit von E. C. Bingham 
und R. W. Wikoff! gefunden, die den Ausflu8 von Seesand durch Glaskapillaren (D= 0316 bs 
0,36 cm, Lange 20 cm) gemessen haben. Das wiirde hier dem Grenzfall Lochdurchmesser — Rohr- 


Ke Tet Sai and R. W. Wikoff, The flow of dry sand through capillary tubes. J. Rheology 2 (1931), 
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durchmesser entsprechen. Sie faBten die Ergebnisse in einer empirischen Formel zusammen, 
nach der Q ~ D®*. Vermutlich wire auch hier eine Darstellung durch Q ~ (D — k)*5 
treffender. 

b)AusfluBauseinem geneigtenRohr. Das Verhaltnis des Ausflusses aus einem 
geneigten Rohr (Winkel / gegen die Senkrechte) zu dem aus dem gleichen, aber senkrecht ge- 
stellten Rohr ist praktisch unabhangig vom Lochdurchmesser und hangt nur wenig vom Rohr- 


durchmesser ab. Nach der obigen Abschatzung ist v ~ |g d; wenn man jetzt g ersetzt durch 
gcos /, so ware zu erwarten Q(()/Q(0) = Jcos 6. In Wirklichkeit ist jedoch nach Abb. 10 iiber- 
raschenderweise der AusfluB bei kKleinen Winkeln / etwas gré®er als aus dem senkrechten Rohr. 
Offenbar ist hier die kleinere Reibung zwischen Glas und Sand entscheidend, da beim geneigten 
Rohr ein gréRerer Teil des Schlauches mit flieSendem Sand von der Rohrwandung gebildet wird, 
— vor allem in einem engen Rohr. Es kénnte aber auch eine Folge des vergréBerten Porenvolumens 
sein, durch die der Ausflu8 durch kleineren Druck erleichtert wiirde; in diesem Fall wire ein 
kleinerer Effekt zu erwarten bei Seesand mit runden Kérnern und geringerer Dilatanz. 

Bei starker Neigung des Rohres, etwa bei 6 = 90°—  (m = Béschungswinkel) hért der Aus- 
fluB natiirlich ginzlich auf. Davor tritt bei etwa § => 40° eine Art ,,Kavitation’ auf. Im oberen 
Teil des Rohres rieselt der Sand dann so langsam nach, dafB weiter unten das Volumenverhaltnis 
Sand zu Luft zu stark unter den Normalwert fallt und schlieBlich eine Luftblase entsteht und 
hochsteigt, was an einer wellenférmigen Bewegung des Sandes zu erkennen ist. (Es sind dies 
nicht etwa Luftblasen, die am unteren Rohrende durch das Loch in das Rohr gekommen waren.) 


6. Versuche, das Geschwindigkeitsfeld auszumessen, Es wurde schlieBlich noch versucht, 
das Geschwindigkeitsfeld zu ermitteln durch Ausmessen der Verformung von Schichten ge- 
farbten Sandes. Der Ausflu8 aus einem weiten Glasrohr (D = 17,6 cm) durch ein Loch wurde 
dazu etwa aller halben oder ganzen Minute unterbrochen. Der Schlauch mit rieselndem Sand ist 
unten, iiber dem Loch, kegelférmig (gesamter Offnungswinkel: grober Sand 25,2°, feiner Sand 
16,8°.) Wenn die Sandhéhe groB genug ist, ist der Schlauchdurchmesser oben nahezu konstant 
(vor allem in feinem Sand). Dementsprechend nimmt die Sandgeschwindigkeit auf der Mittel- 
linie nach oben in Lochnahe quadratisch ab mit der Hohe, und ist im oberen Teil nahezu konstant. 

Um mehr Einzelheiten zu ermitteln waren weitere Versuche dieser Art noétig, die jedoch durch 
den Wegzug des Verfassers unterbrochen wurden. 


(Eingegangen am 24. Oktober 1951.) 
Anschrift des Verfassers: Dr. Wieghardt, ,, Uplands‘: Tilford, Surrey (England). 
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Der Einflu® der Querschnittsform auf das Verhalten 
von Bourdonfedern. 


Von W, Wuest. 


1, Einfiihrung, Kreisformig gekriimmte Réhrenfedern von ovaler Querschnittsform, sog. 
Bourdonfedern, werden in der Technik als DruckmeSelemente benutzt und in groBen Mengen 
erzeugt. Wahrend man urspriinglich eine ellipsenahnliche Ovalform bevorzugte, ist man in den 
letzten Jahrzehnten, wohl vorwiegend aus fertigungstechnischen Griinden, immer mehr zur 
Flachbogenform (vgl. Abb.1) iibergegangen. Bei gegebener Querschnittsform (z.B. Ellipse, 
Flachbogen u.a.) ist das Verhalten diinnwandiger Réhrenfedern vor allem durch das Achsen- 
verhaltnis « und den Parameter A= a?/sr (a groBe Querschnittshalbachse, s Wandstarke, r Kriim- 
mungshalbmesser) bestimmt. Fiir Querschnittsformen, die abschnittsweise aus Kreishégen oder 
Geradenstiicken zusammengesetzt sind (Korbbogen 
mit den Sonderfallen Spitzbogen und Flachbogen, 
vgl. Abb.1) sind Berechnungsgrundlagen von M. Tueda' 
entwickelt worden. Eine erhebliche Vereinfachung 
dieser Theorie ist vom Vf. vorgenommen worden und 
die Rechnung fiir den praktisch fast ausschlieBlich 
interessierenden Flachbogen numerisch durchgefiihrt 
und durch zahlreiche Messungen nachgepriift worden ®. 
. Neuerdings haben R.A.Clark und E. Reissner® eine 
Theorie der Réhrenfedern von elliptischem Querschnitt 
entwickelt. Die Ergebnisse dieser Theorie stehen je- 
doch in Widerspruch zu Messungen, die von einem 
fritheren Mitarbeiter des Verfassers mit Réhrenfedern 
von elliptischem Querschnitt durchgefiihrt worden sind. Die bei den Messungen verwandten 
Réhrenfedern sind zunachst mit elliptischen Ziehdornen gezogen worden. Doch besteht immer- 
hin die Méglichkeit, daB bei dem ,,Rollvorgang‘‘, bei dem das urspritnglich gerade Rohr eine 
kreisférmige Kriimmung erhalt, die elliptische Form nicht streng erhalten blieb. Fiir eine 
Entscheidung iiber die Giiltigkeit der Reissnerschen Theorie wird man also weitere sorgfaltige 
Messungen abwarten miissen. Solche Messungen sind auch deswegen von praktischem Interesse, 
weil die Reissnersche Theorie einen unerwartet hohen Einflu8 der Querschnittsform anzeigt. 
Méglicherweise kann man bei diinnwandigen Réhrenfedern Querschnittsformen mit giinstigeren 
Kigenschaften finden. 


/ 


Korbtogen  ‘Spitebogen  Hachbagen Elijpse 


Abb. 1. Verschiedene Querschnittsformen 
von Rohrenfedern. 


Kine erhebliche Vereinfachung der Theorie ist in folgenden beiden Fallen méglich: 
1. A klein (etwa < 1), a beliebig, 
2. A beliebig, a klein. 


Der erste Grenzfall entspricht der ,,Mitteldruckfeder‘‘, d.h. bei tiblichen Abmessungen ist diese 
Vereinfachung zulassig fiir Réhrenfedern, deren Nenndruck (GréBtwert des Innendruckes) 
groBenordnungsmabig 50 at itberschreitet, und bis zu einigen 100 at betragen kann. Die Mittel- 
druckfedern liegen also im Grenzgebiet zwischen den diinnwandigen Niederdruckfedern und den 
dickwandigen Hochdruckfedern. Die Berechnung von Mitteldruckfedern ist in vielen Fallen 
analytisch méglich*. Fiir beliebige Querschnittsformen haben C. B.Biezeno und J.J. Koch® ein 
numerisches Differenzenverfahren entwickelt, dessen Giiltigkeitsbereich durch mehrfachelteration 
bis etwa J=3 ausgedehnt werden kann. 


+ M. Tueda, Mem. Coll. Engin. Kyoto (Japan) (1934/35), H. 8 und (1935/37), H. 9. Die dort abge- 
leiteten Gleichungen enthalten jedoch einen Fehler. 
a ecb nlose der Rechnung s. W. Wuest, Technik 3 (1948), S. 23; vgl. auch Bild 4 der vorliegenden 
rbeit. 

8 R.A. Clark und E. Reissner, J. appl. Physics 21 (1950), S. 1340. R.A. Clark, T. I. Gilroy und 
E, Reissner, J. appl. Mech. Paper No. 51—A—11 (1951). 

* Einige Berechnungsbeispiele finden sich in Technik 3 (1948), S. 24, Bild 3. 

° C. B. Biezeno und J. J. Koch, Proc. Amsterdam 44 (1940), S. 779 u. 914. 


XX. Band 1952. Wuest: EinfluB der Querschnittsform auf das Verhalten von Bourdonfedern. 1s 


Auch im zweiten Grenzfall, namlich fiir sehr kleines Achsenverhaltnis des Querschnitts tritt 
eine erhebliche Vereinfachung ein. Dieser Grenzfall soll im folgenden naher untersucht und auf 
die Berechnung beliebiger Querschnittsformen angewandt werden. VeranlaBt wurde diese Unter- 
suchung durch die Unstimmigkeiten zwischen der Reissnerschen Theorie und den oben genannten 
Messungen. Aus diesem Grunde werden auch vor allem ellipsenihnliche Querschnittsformen 


untersucht. Die Ergebnisse dieser Naherungsrechnung bestatigen eher die Messungen als die 
Reissnersche Theorie. 


2, Ableitung der Grundgleichungen und der Randbedingungen, Diimnwandige Réhrenfedern 
fallen vom Standpunkt der Festigkeitslehre unter den Begriff Umdrehungsschalen. Jedoch ist 
der Ring bei der Réhrenfeder nicht geschlossen, sondern weist gewdhnlich einen Windungs- 
winkel von etwa 270° auf. Die Réhrenfeder ist am einen Ende im Federtrager (Abb. 2a) 
befestigt, der zugleich die .— én ~ 
Druckzufithrung enthalt, 
wahrend das freie Ende der 
Feder durch ein Federend- 
stiick verschlossen ist. Wir 
setzen voraus, dah die 
Réhrenfeder nur durch 
einen gleichmaBbigen Innen- 
druck p und allenfalls durch 
ein am freien Ende ein- 
wirkendes Moment M be- 
lastet ist. Die Belastung 
ist also kreissymmetrisch 
und daher auch die Span- 
nungsverteilung abgesehen 
von den beiden Enden, 
unabhangig vom Winkel y 
(Abb.2a). Der Querschnitt AE EE ae sa 
(Abb.2b) wird als doppelt- 
symmetrisch vorausgesetzt. (Diese Voraussetzung ist fiir die folgende Rechnung zwar nicht not- 
wendig, doch ist sie in fast allen praktischen Fallen erfiillt.) Als Achsenlangen a, b, der Mittel- 
flache bezeichnen wir die Achsenlange des umschriebenen Rechtecks. Die Mittelflache des 
doppelt-symmetrischen Querschnitts sei durch folgende Funktion gegeben: 


6 fe eee es bf (x). 


wobei f(x) eine gerade Funktion von x ist. Das Achsenverhaltnis a= b/a wird als sehr klein 
vorausgesetzt und soll im Grenzfall gegen 0 gehen. Wegen dieser Annahme fallen die tangen- 
tialen Verschiebungen wu und die normalen v (vgl. Abb.2b) im Grenzfall in die Achsenrichtungen, 
wodurch sich erhebliche Vereinfachungen ergeben. Wegen der vorausgesetzten doppelten Sym- 
metrie der Querschnittsform kénnen wir die Rechnung auf einen Quadranten beschranken. 
Bezeichnen wir den Abstand von der Mittelflache mit z (—h<z<h), so erhalt man fiir die 
tangentielle Dehnung ¢, und die Ringdehnung «¢, im Grenziibergang «—>0 


zur Krimmungsachse 


du d?v 
— — 7 —_, 1 
Cea ie At ge (1) 
e v 
fy = —— a. (2) 


Hierbei ist w)—= Ay /) die relative Winkelanderung, die als positiv angenommen wird, wenn 
die Feder sich aufbiegt, wie es normalerweise unter der Wirkung des Innendruckes geschieht. 
Die tangentiale Spannung o, und die Ringspannung o, ergibt sich dann zu 


On == is (x a [Es) ” 


7 


Oy — 1 aay, (é, + ia Ex) . 


Durch Integration iiber die Wanddicke erhalt man daraus die Zugkrafte T, und T,, sowie das 
ga 
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Moment M,: 
__ 2hE [du Bib 3 
iN ater er poet r on)| (2) 
2hE [| v y du\ if 
Gis ere (2 — 7, Mo + mu s) (4) 
BRE de 
Me 5 dak (5) 
Im Grenzfall ~— 0 geht auch T,— 0. Daraus folgt 
d v 
4 (* — ay). (6) 
Bezeichnet man mit N die Querkraft, so erhalt man als Kraftegleichgewicht in der y-Richtung 
dN : 
Pop le ph (7) 
und als Momentenbedingung : 
IM a (8) 
dx 
Durch Elimination von N erha]t man aus (7) und (8): 
Mx 


a eg Ly pr 
und durch Einsetzen von (4), (5) und (6) schlieBlich die gesuchte Differentialgleichung 


d4y 3 3 (1 — p2) 
oP spa = 0") h2 2 (o—YO)e ORE 


Fiihrt man noch die dimensionslosen GréBen §= x/a und w= v/a ein, geht diese Gleichung 
itber in 


wl¥(8) + 4 ntw = 4 nt [ae f(8) — A] (9) 
mit den Abkiirzungen 
Ree Oe 2 (jl — Wace 
4 nt = (1 — p?) Phe? A= — = od e ; (9a, b) 
Die Randbedingungen lauten 
bee b= we 0 
Aus (4) und (6) findet man die Ringkraft 
2hE 
7, = AE (yy ay). (10) 


Das Moment der Ringkrafte mu aber gleich dem von aufen auf die Feder einwirkenden Moment 
sein, d.h. 


1 
a SRE 
M=4/ T,ydx=8a°"™ | (w—apf) fae - (11) 
; Q 
Bei einer nur durch den Innendruck belasteten Feder ist M0, so daB man die relative Winkel- 
anderung w, berechnet zu 


1 
Wy) = —- ; (12) 


3. Allgemeine Lésung und Anwendung auf Sonderfille. Die allgemeine Lisung der Differential- 
gleichung (9) ist 
w= A(1+ C, oj n§ cosnE+ C, Gin né sin né + naw, jw (& — &’) f(&’) dé’ , | 
| (13) 
mit 


w* (£) = Cpj né sin nE — Gin né cos nE. 
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Hierbei ist schon beriicksichtigt, da die Lésung eine gerade Funktion in & sein muB. Die beiden 
Integrationskonstanten C, und C, ergeben sich dann aus den Randbedingungen w= w’—0 fiir 
§ = 1. In vielen Fallen kann aber die Lésung explizit angegeben werden, insbesondere in folgenden : 

a) f(&) durch trigonometrische oder Exponentialfunktionen gegeben, 

b) f(€) durch ein Polynom in & gegeben. 

Weiterhin soll insbesondere der Fall b) naher untersucht werden, und zwar werden wir uns 
auf ein Polynom sechsten Grades beschrinken. Die Funktion f(&) sei also gegeben durch 

F(E) = d+ 0,6" + ay&4 + a8. (14) 

Die Lésung der Differentialgleichung (9) ist dann 


w= A(1-+ C, Gof n& cos né + C, Gin né sin n&) 4 (Fr) . a, a8) 6.06%? (15) 


Die Integrale in (11) und (12) werden 
1 


eh, A 2 2 I 2 2 
| fede a3 3 4% % 5% UG 7 40% 5 4 7 Oy % bo Me Me 


9 
3 (16) 
ad Una, 2 c188 | 
74 cE hak oes 13 “6? 
1 1 
“A 6 
[wae oo [pea + (4% a, a0) (a, a | a ss fp 4 (C.F CF) | (17) 
0 0 ; 
90 a a, a a, 
_ a, a0 ( 3 oo a oi | 
wobei 


F, (n) = te (Sin n cos n+ Cpj nsin n) (a) + a,-+ a, + a,) 


= = Sin n sin n (a, + 2a, + 3a,) (Sin n cos n— pj nsin n) (a,/2 + 3.a,-++ 15/2 a) 
+- = oj n cos n (6a, +- 30a.) — = (Sin n cos n + Cj n sin n) (3.4, + 45 ag) 
Le tee : 
+ >e« Sin n sin n 90a, + 7 (Sin n cos n — oj nm sin n) 45 a, (17a) 
F,(n)= — = (Sin n cos n — Cpj n sin n) (a, + a, + a, + ag) 


+ fy Cof m cos n (ay + 2ay-+ 3a) — 1, (Ginn cos n+ Gof sin n) (ay]2-+ 3a, 4 15/249) 


= & Sin n sin n (6a, + 30a,) + = (Gin n cos n — oj n sin n) (3 a, + 45 ag) 


— se Co} n cos n 90a, + is (Sin ncos n+ Coj n sin n) 45 a, , (17b) 
ist. oe 
4, Kenngrofen der Réhrenfeder. Bevor wir zur zahlenmafigen Berechnung von Beispielen 
iibergehen, sollen aus den im vorigen Abschnitt entwickelten Beziehungen die fiir die praktische 
Anwendung wichtigen Kenngréen der Réhrenfeder abgeleitet werden. Ausgangspunkt ist 
einerseits die Gleichung (11), die in folgender Form geschrieben werden kann: 
oy = opp + euM (18) 
Andrerseits ist fiir die Beanspruchung vor allem die gréBte Biegespannung! mabgebend, die an 
den Enden der grofen Achsen erreicht wird, namlich 
h E 


OR w’’ (1) a.t—,- 


1 Fiir groBe Werte 2 kann die Zugspannung am Ende der groBen Achse wesentlich werden. Sie ist 


1 II WE 
offenbar nach Gl]. (5) und (8) gegeben durch oF INTL) are ((19) warmer, 
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Diese Beziehung laBt sich auf folgende Form bringen: 
Op= d, p+ 4, Wy | 
= (d, + d,¢,) p+ dycuM . | 
Die Koeffizienten c,, cy, d,, d, kénnen aus der in Ziff. 3 entwickelten Lésung errechnet und 


P w 
auf folgende Form gebracht werden: 


my Open A Sle ee 
= (1 77) ee aa Leas (20) 
Hierbei ist s—2h die Wanddicke und die Funktion F', auBer von der gewahlten Querschnittsform 
nur vom Parameter A= a?/rs beziehungsweise von n= 3 (1 — ww?) 2 abhangig. Der Koeffi- 
zient c, ist maBgebend fiir die ,.Druckempfindlichkeit‘+ der Réhrenfeder. Der Praktiker inter- 
essiert sich statt dessen mehr fiir den vom Federende zuriickgelegten Weg f, der daraus berechnet 
wird zu 


(19) 


f= ep YokP - (21) 
Hierbei ist y) der Windungswinkel der Feder und k der Abstand des Federendes vom ,,Polar- 
punkt*: der Feder!. Bei dem haufig angewandten Windungswinkel y= 3/27 1st _ k= 125 7, 
Die GroBe 1/cey wp ist ein Ma® fiir die Federsteife und soll als spezifisches Richtmoment 
bezeichnet werden. Man kann diese GréBe auf folgende Form bringen: 


(22) 


wobei das geometrische Tragheitsmoment J sich errechnet zu 
1 


dj OS | peas ; 
0 
Der Koeffizient F, gibt also Auskunft dariiber, inwieweit die Steifigkeit eines Hohlquerschnitts 
gegeniiber einem Vollquerschnitt von gleichem geometrischen Tragheitsmoment vermindert wird. 
Fiir kleine Werte von / nahert sich F., dem Wert 1. 
Den Koeffizienten d,-+ d,,c, fiir die gréBte Biegespannung unter der Wirkung des Innen- 
druckes schreiben wir in der Form 


2 
a 
d, +d, ¢, = | ") ee (23) 
SchlieBlich fithren wir noch ein 
d,/ (dp -- d,.¢,) = F,. (24) 
Die Funktion F, hat eine wichtige Bedeutung fiir die Berechnung der Uberdrucksicherheit der 
Réhrenfeder durch Abfangen und der Belastungsfahigkeit von dynamischen Differenzdruck- 
messern *. Die Gleichung (19) fiir die maximale Biegespannung kann namlich auf folgende Form 
gebracht werden: 
ai ( Sk 
on=(*) F,(F, p+ Ss * Wo) « (25) 
Bei diinnwandigen Réhrenfedern wird eine beschrankte Uberdrucksicherheit dadurch erzielt, da8 
die Feder bei einem bestimmten Enddruck abgefangen wird, also an einer weiteren Kriimmungs- 
anderung gehindert wird. Wenn F,< 1, wie es bei diinnwandigen Réhrenfedern immer der Fall 
ist, steigt die Beanspruchung nach dem Abfangen mit steigendem Druck geringer an, als es ohne 


Abfangen der Fallsein wiirde. Findet das Abfangen bei einem Druck p* statt, so ist die maximale 
Biegespannung fiir p > p* 


On = Bre Pe +F,(p—p*). (26) 


Dieser Vorgang ist in Abb. 3 dargestellt. Als Beispiel ist eine Feder gewahlt, die bei einem Nenn- 
druck von 2,5 kg/cm? eine maximale Biegespannung von 25 kg/mm? erreichend mége. F, wird zu 
0,2 angenommen. Die gleiche Feder werde jetzt bei 1,6 kg/em? Innendruck abgefangen. Der 


Druck kann jetzt bis auf 6,1 kg/em® erhéht werden, ehe die gleiche Beanspruchungsgrenze er- 
reicht wird (Kurve a’ in Abb. 3). 


' W. Wuest, Z. Instrumentenkunde 63 (1943), S. 416. 


* W. Wuest, Arch. Techn. Messen V 1343=5 (Juni 1948). 


XX. Band 1952. 


Wuest: Kinflu8 der Querschnittsform auf das Verhalten von Bourdonfedern. 121 


Eine weitere Anwendung bezieht sich auf dynamische Differenzdruckmesser, bei denen zwei 
gleiche Federn fest miteinander verbunden sind, so da eine Kriimmungsanderung nur bei einer 
Differenz der Richtmomente eintritt, wahrend sie sich bei gleichem Innendruck gegenseitig in der 
Bewegung hindern. Bei ungleichem Druck betragt natiirlich der Federweg nur die Halfte des 
Weges einer entsprechenden freien Feder. 
Bezeichnet man den niedrigeren Druck mit 
Pi den héheren mit p,, so ist die gréBte Biege- 
spannung bei der héher beanspruchten Feder 
offenbar 


max. Bregespannung 


os=(*) F, [Fmt i""(p.— py, 20 


ivemyD 
os=(“) FFP, 5) “(Ps P) (27a) 


Baut man mit den gleichen Federn wie im 
obigen Beispiel einen Differenzdruckmesser 
fiir 2,5 kg/cm? Differenzdruck, so darf, wie 
man aus Abb. 3 ersieht, p, bei gleicher zu- 
lassiger Beanspruchung den Wert po Ug Abb. 3. Erlauterungsbeispiel fiir die Uberdrucksicherheit von 


: 4 Roéhrenfedern und das Verhalten von dynamischen Diff - 
kg/cm? erreichen. (Kurve a—b’ in Abb. ale 5 Laie ee aa ay ynamischen Differenz 


5. Berechnungsbeispiele. a) a)—=1, a, —a,—a,—0. Dies entspricht also dem Grenzfall des 
Rechteck- oder Flachbogenquerschnittes von kleinem Achsenverhaltnis. Man erhalt hier fiir die 
Funktionen F,, F,, F, und F, das einfache Ergebnis 


Fa=3tF™—!} F 1 _ 3 Gin2n—sin2n 
i 


Lae ; 1S n F(n)’ 2 n Coj 2n—cos 2n’ 


mit 


Gin2n+sin2n 
Om pena pyraa 


Of 02 03 OF 7 E85 GW 2030 50 0 Gf O203 OF 7 é 5 10 2 30 50 100 
i 


NA 
Gf G203 05 1 2@ 8 5 OW 238 50 100 Wh WEE GE PB 5 0 2030 50 0 


Abb. 4a, b, ec, d, Hilfsfunktionen zur Berechnung yon Réhrenfedern mit Flachbogenquerschnitt. 


In Abb. 4 sind diese Funktionen in Abhangigkeit von 
— & ee ere os == 
au! y 31 — p?) 


aufgetragen. Aus einer unverdffentlichten Arbeit des Vf. sind diese Funktionen fiir den Flach- 
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0 G2 gy 46 O8 40 


Abb. 5. Querschnittsformen, die Polynomen zweiten Grades ents prechen 
(f(—) = 1—f & nach Rechenbeispiel bj, 


bogenquerschnitt auch bei endlichen 
Achsenverhaltnissen entnommen 
worden!. Mankann aus dieser Gegen- 
iiberstellung einen Anhalt gewinnen, 
wie Weit die Grenzkurven fiir « = 0 
ein qualitatives Bild des allgemeinen 
Falles ergeben. Die F,-Kurven ver- 
schiedener Achsenverhaltnisse wei- 
chen voneinander weniger ab als die 


1 Die vier Diagramme sind bereits 
in einem friiheren Aufsatz [Technik 3 
(1948), S. 23] veréffentlicht worden. 
Wegen eines dort vorhandenen Be- 
zeichnungsfehlers werden sie hier erneut 


abgebildet. 


L(A) 
G5 


a OF 


G3 
G2 = 


Mehwerte 
Elipt. Querstinitt 
x @=979 
° o=O1 


G01|———-+}~- ——}-—-}--}- 


i] 2 eA lb) 20 40 60 100 7 a 


bo Pe GD) 20 40 60 100 


008 


006 


Of lea N 
7 ‘e 4 x rep fd) 20 ¥O 60 80100 7 2 ¥ OVE. 20 40 60 80700 
Abb. 6. a, b, c, d. Hilfsfunktionen zur Berechnung von Réhrenfedern verschwindend kleinen Achsenverhaltnisses 


mit Querschnittsformen nach Abb. 5, 
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iibrigen Kurvenscharen. Es ist sogar zu erwarten, da sie bei einem elliptischen Querschnitt 
noch enger zusammenliegen, wie man aus dem Verhalten bei kleinen J-Werten schlieBen kann. 
Auch die Ergebnisse der Reifnerschen Theorie, die allerdings nicht gesichert sind, weisen in 
diese Richtung (vgl. Abb. 6a). 

b) ag=1, a,= —f, a,= a, =0. In Abb. 5 sind verschiedene Querschnittsformen fiir K=O Oe: 
0,5; 0,75 und 1 aufgetragen. Der Wert / — 1 entspricht angenahert einem spitzbogenférmigen 
Querschnitt, = 0,5 im mittleren Bereich einer Ellipse (gestrichelt eingezeichnet). Die Funk- 
tionen F’,, F,, F’, und F’, nehmen eine einfachere Gestalt an, wenn n = mz (m ganze Zahl) namlich 


hs 1 — 6/3 — B/n* — (1 — 8B) Tg n/n 
Hil") = 2 app) a ECL By" o? Sg hE D fin Cig hee 
_ 260 — 8) n*? + (1 —8)? nn? Tg n + 2 B2 n Ctg n — f? 
Fo) Paes Ry 
F,(n)= 3/n®? + (2n8 Ctg n-+ n? (1 — 6) — £) F,, 
F, | 


Fam) = 1/{1-+ | anf Gig- nmr —p— py]! 


Wenn Tg nund Ctgngegenseitig vertauscht werden, gelten diese Formeln auch fiir n= (2 m+ 1)z7/2. 
Fiir n= 0 ergeben sich folgende Grenzwerte: 


ee Afl5 Gh — 6/7) eee te 2} 
be 1—23/3p +R’ Poe bal t=. 
In Abb. 6a, b, c, d sind die vier Kurvenscharen in Abhangigkeit von J und (3 aufgetragen. (6 = 0,5 
stellt eine grobe Naherung fiir die Ellipse dar. Es ist daher in Abb. 6a eine bessere Naherung dar- 
gestellt, bei der die Ellipse durch ein Polynom sechsten Grades angenahert wird. Ferner sind 
noch die in der Einleitung erwahnten Messungen eingetragen, die in der Nahe dieser Grenzkurve 
liegen. Die gleichsfalls eingetragenen Kurven nach E. Reifner liegen dagegen etwa um den 
Faktor 2 héher. In Abb. 6c sind ebenfalls die Werte nach Reissner eingetragen. Allgemein lassen 
die Kurven erkennen, daB bei Wahl einer mehr spitzbogenférmigen Querschnittsform die Feder- 
bewegung bei sonst gleichen Abmessungsverhdltnissen bedeutend vergréBert werden kann. 
Andrerseits sinkt aber auch die Steifigkeit der Feder und wachst die Beanspruchung, so daB der 
Gewinn im ganzen doch zweifelhaft ist. 

c) Als drittes Beispiel sollen Federformen untersucht werden, die nach auben wulstartig ver- 
breitert sind (Abb. 7) und durch ein Polynom vierten Grades dargestellt werden sollen. Wahlt 
man 


A 4p __ 48 (8 —1) ee. Ap? _ 
SSATP Sl ae emma see oe ah ame TD es 
so ist die Federform gegeben durch 
FUP +e) (62) 


Aldo 
mit einem Gré®twert (1 + £)?/46. Vergleichsweise ist auch die schon frither behandelte Kurve 
mit Grenzwert dem {= 0 (vgl. Beispiel b) eingetragen. Die Rechnung liefert bei n= ma 
(m ganze Zahl) 


B elb pera lite 
F,(n) = — a te: : ee ae ele 
vty 4 
2(1+ pn tan + its Tgu—OF _y, B 
mit 
2 12 1+f 
N= 48(1+ 2)4 G2 fy 36) AG B) 
Ferner ist 
| 
3(1 +> 
F,(0) = | 21 (L+ 8? 
Wy canes 
Pig: 


124 Wuest: Einflu8 der Querschnittsform auf das Verhalten von Bourdonfedern. Ingenieur-Archiv 


Boe 
53 
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Abb. 7. Querschnittsformen, die Polynomen vierten Grades entsprechen. 


Cf (€) /f (0) = (1 — &) (14+ f &) nach Rechenbeispiel c]. 


vii 


O2 


G0 T + 


a 
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Abb. 8. Hilfsfunktion F, zur Berechnung der Druckempfindlichkeit von 
Réhrenfedern mit Querschnittsformen entsprechend Abb.7. 


s7(E) 


10 
| is 
Ob + + 
G2 \ 
4 


Abb. 9. Querschnittsformen, die Polynomen sechsten Grades entsprechen 
(Rechenbeispiel d). 


Die sich hieraus ergebenden Werte 
fir F, sind in Abb.8 in Abhangig- 


keit von 


A= n?//3 (1 — p?) 
aufgetragen. Bei VergréBerung des 
Wulstes nimmt die Druckempfind- 
lichkeit der Réhrenfeder ab. Bei 
Querschnittsformen mit schwachem 
Wulst ist die Druckempfindlichkeit 
gréBer als bei einer Ellipse mit 
sonst gleichen Abmessungen. 

d) Als viertes Beispiel sollen noch 
die in Abb. 9 dargestellten Quer- 
schnittsformen behandelt werden, 
die als Polynome sechsten Grades 
ausgedriickt werden kénnen. Die 
Kurve a stellt eine méglichst gute 
Annaherung an eine Ellipse dar und 
entspricht folgender Funktion: 

ace Ae 
BE ors: Si LOM 
Die ersten drei Glieder entsprechen 
der Reihenentwicklung einer Ellipse, 
das vierte Glied ist etwas vergréBert. 
Die Kurven 6 und c entsprechen 
folgenden Polynomen: 


&2 
P= 1 yas 
(0,48— 1,335 & + £8) 


bzw. 


C= 


1 2 
~ 0,1912 
(0,6912 — 1,5 £2 + £4). 


Die Funktion F, ist fiir diese drei 
Federformen in Abb. 10 aufgetragen. 
Form c) hat eine gute Druckempfind- 
lichkeit, die iiber derjenigen der 
Ellipse liegt. 


Wr 


6. Zusammenfassung. Die fiir 
die Beurteilung einer Réhrenfeder 
(Bourdonmanometer) wichtigen 
KenngréfBen (Druckempfindlichkeit, 
Steifigkeit, gréBte Biegespannung 
und Uberdrucksicherheit) lassen 
sich auf vier Funktionen zuriick- 
fiihren, die auBer von der Quer- 
schnittsform nur noch vom Achsen- 
verhaltnis und vom Parameter 


J 10 oot, 


(a groBe Halbachse, s Wandstarke, 
r Kriimmungsradius) abhangen. 
Allgemeine Berechnungen sind bis- 
her nur fiir den ,,Flachbogen“ 
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durchgefithrt worden und hier im Ergebnis mitgeteilt. Im Grenzfall kleiner Werte von J und 
beliebiger a ist fiir beliebige Querschnittsformen eine graphische oder numerische Lisung 
modglich. Ebenso tritt eine wesentliche Vereinfachung fiir kleine « und beliebige A ein. Dieser 
Grenzfall wird hier im wesentlichen behandelt und zwar fiir Querschnittsformen, die sich 


i 


G07 + 


— 
7 2 4 (eth dy 20 4O 60 


Abb. 10. Hilfsfunktion F, zur Berechnung der Druckempfindlichkeit 
von Rohrenfedern mit Querschnittsformen entsprechend Abb. 9. 


als Polynom sechsten Grades darstellen lassen. Die fiir eine Naherungsellipse durchgefiihrte 
Rechnung zeigt befriedigende Ubereinstimmung mit MeBergebnissen, wahrend die Reissnersche 
Theorie der elliptischen Réhrenfeder erheblich davon abweicht. Im Ganzen zeigen die durch- 
gefiihrten Rechnungen einen erheblichen EinfluB der Querschnittsform. 


(Eingegangen am 28. Oktober 1951.) 


Anschrift des Verfassers: Dr. W. Wuest, Géttingen, Bottingerstr. 6/8. 
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Die Formanderungen idealplastischer statisch bestimmter Balken. 
Von H. Craemer. 


In der Plastizititslehre bestehen teilweise Unklarheiten tiber die Formanderungen von Balken, 
besonders betreffend den allmahlichen Ubergang vom rein elastischen zum rein plastischen Zu- 
stand und die GréBe der Durchbiegungen und die Form der Biegelinie im vollplastischen End- 
zustand. Als Beispiel, an dem sich das Wesentliche zeigen 148t, behandeln wir einen Balken auf 
zwei Stiitzen mit Rechtecksquerschnitt und zwei symmetrisch angeordneten, gleich groBen Ein- 
zellasten. 

Wir legen dabei die bekannte idealisierte Spannungs-Dehnungsfunktion 


C= Bertir J0.50675 E SEs, | (1) 
o=o, fiir ¢,<&e <0 | 

zugrunde, in der g und ¢ Spannung und Dehnung, o, die Streckspannung und e, die Dehnung bei 

Streckbeginn bedeuten; wir nehmen ferner allmahliches, nicht plétzliches Eindringen des Flie8- 

vorganges in den Querschnitt an. Wird in einem Rechteck bh unter reiner Biegung die Rand- 

spannung 0,=0,, also die Randdehnung ¢, >¢,, so verbleibt ein elastischer Kern h, <h,, 


wahrend oben und unten je ein plastiziertes Randgebiet von der Dicke = (h — h,) und der 


konstanten Spannung go = o, entsteht. Das vom Querschnitt itbertragene Moment ist dann 
M=bo,| 7h +— (hh) (b+ h)). 
Da ferner 
7a ) 
ist, erhalt man hieraus 
Nie pn?(3 —%), 


insbesondere bei voller Plastizierung mit ¢,— oo, M, =o, b - ; 
Definiert man als Ausnutzung des Querschnitts die GréBe »— M/M,, so wird 
ee 3 (3) 
bei FlieBbeginn ist die Ausnutzung y= 2, am SchluB y= 1. 
In Abb. ] ist nun, stark umrandet, ein Teil des Balkens dargestellt. Ist 7) die Ausnutzung des 
Querschnitts unter der Last, so ist sie fiir einen beliebigen Querschnitt »,, entsprechend dem Ver- 


lauf der Momentenlinie; dies ist fiir verschiedene Laststufen <v) <1 oben links bei a ein- 


getragen. Die fiir jeden Querschnitt und jede Laststufe nach (3) hierzu errechneten GréBen 
é,/é, sind bei b zu sehen. Ferner sind nach (2) die Héhen der elastischen Kerne ermittelt und 
bei c eingetragen; man sieht, wie der FlieBvorgang sich einerseits nach der Tiefe, andererseits 
in der Langsrichtung ausbreitet. 
Die gegenseitige Verdrehung zweier benachbarter Querschnitte ist nun (siehe d) 
dx 

ihn = Er : 

eas | 
die Biegelinie ist dann bekanntlich die Momentenlinie des Balkens bei Belastung mit diesen sog. 
Winkelgewichten. Wir belasten nun statt dessen mit 

h t.) UF 

fe ae dan 

(siehe die schraffierten Differentiale in b), wobei dx = 1d& gesetzt wurde, und erhalten dadurch 
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die h/2 le,-fachen Winkelanderungen und h/2 /®¢,-fachen Durchbiegungen an beliebiger Stelle. 
Die so gewonnenen dimensionslosen Werte sind unabhangig von den Abmessungen (h, 1) und 
Baustoffeigenschaften (¢,) des Balkens und stellen den Vorgang als alleinige Funktion des Deh- 
nungsgrades ¢,/e, oder, wegen (3), der Ausnutzung dar. 


Bei dieser Auswertung, die rechnerisch vorgenommen wurde, ist das Verhaltnis Lastabstand: 
Spannweite mit « —= 0,01 angenommen worden. Bis etwa vy = 0,9 erhalt man praktisch dasselbe 
wie fiir eine Kinzellast in Feldmitte. Dann aber wachst der Einflu8 derjenigen Winkelgewichte, 
die zwischen den Lasten infolge der dort herrschenden groBen Randdehnungen entstehen, vgl. 
Skizze b. Bei voller Ausnutzung, yj — 1, hat die Flache der ¢, zwischen Auflager und Last zwar 
eine unendlich groBe Ordinate, aber, wie ich friiher auseinandergesetzt habe, einen endlichen 
Flacheninhalt. Mit «= 0 wiirde man also eine endliche Durchbiegung erhalten, und zwar eine 
solche, die weit geringer ist als sie bei Versuchen beobachtet wurde. Zieht man dagegen eine noch 


pon/es 


Abb. 1. 


so geringe, wenn auch endliche Spreizung der beiden Lasten oder, was im Ergebnis auf das Gleiche 
hinauskommt, eine Verteilung tiber eine geringe Lange in Betracht, so wird unter 7) = 1 das 
Winkelgewicht zwischen den Lasten und infolgedessen auch die Durchbiegung unendlich gro8. 

Freilich ist unter so groBen Verformungen die Beanspruchung des Balkens nicht mehr so wie sie 
hier zugrunde gelegt wurde und auch der Zusammenhang zwischen Winkelanderung und Durch- 
biegung ist anders, d.h. man miiBte eigentlich die Untersuchung durch Entwicklung einer 
Theorie zweiter Ordnung verbessern; doch wiirde dies am Grundsatzlichen wohl nicht viel 
andern. 

In Skizze e ist nun die Durchbiegung wy in Feldmitte als Vielfaches von ¢, /?/6 h in Abhangig- 
keit von der Ausnutzung %) aufgetragen. Bis zum FlieSbeginn, 1) = 2/3, ist die Beziehung linear. 
Dariiber hinaus wachst die Durchbiegung schneller als die Ausnutzung; z. B. erhalt man fiir 
¥9 = 0,99 einen Beiwert von 2,03 in der Figur, wahrend man bei Fortdauer der linearen Be- 
ziehung nur 1,48 erhalten hatte. 

Ferner sind in der Skizze die zu jeder Ausnutzungsstufe gehdrigen Dehnungsgrade ¢€,)/é, 
fiir den Querschnitt unter der Last angegeben. Die Dehnungen wachsen also schneller als die 
Qurchbiegungen, da letztere nicht nur von den grofBen Dehnungen in den meist beanspruchten 
Duerschnitten, sondern auch vom Verhalten der minder beanspruchten Schnitte beeinfluBt 
werden. Bei Versuchsauswertungen sind diese beiden Linien teilweise verwechselt worden. Auch 
ist zu betonen, daB® die Form dieser Kurven nicht eine Folge der meist vorhandenen Ausrun- 
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dung der (a, ¢)-Linie in der Gegend von ¢ = ¢, sind, die hier nach (1) ausdriicklich vernachlassigt 
wurde. 

Die Form der Biegelinien wird am besten durch die Verhaltnisse w,/wy (w, an beliebiger 
Stelle, wy unter der Last) gekennzeichnet (siche Abb. 1, f). Man sieht, wie die Kriimmung in 
dem elastisch beanspruchten Balkenteil immer mehr hinter der in den plastizierten Teilen zu- 
riickbleibt. Am Schlu8, unter ») = 1 und h, = 0, ist die Verformung in den Bereichen auBerhalb 
der Lasten, wo v,< 1, unendlich klein gegeniiber derjenigen zwischen diesen, wo ¥,—= = lI. 
Stehen die beiden Lasten so nahe beieinander, daf sie als eine einzige betrachtet werden kénnen, 
so bilden im Endzustand die gerade Linie der Skizze f und die entsprechende auf der linken Bal- 
kenhalfte einen Knick, d.h. das bekannte FlieBgelenk!. Im allgemeinen Falle bildet sich zwi- 
schen den Lasten eine Gelenkkette. 

Man kann die vorstehenden Untersuchungen auch so durchfithren, daB man als wirksam nur 
den elastischen Kern und den auf ihn entfallenden Teilbetrag des Moments ansieht und hierauf 
die Lehren der Elastosiatik anwendet. Unrichtig ist es dagegen, das Moment des ganzen Quer- 
schnitts auf den Kern wirken zu lassen, wie es in einem Buche iiber Plastizitat geschehen ist. 

Die vorstehenden Untersuchungen lieBen sich, an sich ahnlich, wenn auch miihevoller, fiir einen 
I-Trager durchfiihren. Wir beschranken uns auf einige Bemerkungen fiir den ,,idealen I-Trager“, 
d.h. einen solchen, bei dem der ganze Baustoff in den beiden aufersten Fasern zusammen- 
gedrangt ist. Es bleiben dann bis unmittelbar vor dem Bruch (v,<%< 1) alle Querschnitte 
elastisch. Beim Bruch hingegen wird zwischen den Lasten 1) = 1 und auBerhalb derselben wei- 
terhin 7,< 1, es werden also nur die Schnitte zwischen den Lasten plastisch und es gilt sowohl 
fiir wy wie fiir ¢,/e, die gestrichelte geknickte Linie in Skizze f; dies ist aber die gleiche Form 
wie die (a, ¢)-Linie nach (1). 

Wirkliche I-Trager unterscheiden sich nicht allzu stark von dem hier betrachteten Grenzfall, 
so daB der Unterschied im Verlauf der eben erwahnten drei Funktionen bei Versuchen oft kaum 
zur Geltung kommt, zumal die Nebeneinfliisse oft zu gréBeren Unterschieden fiihren. Dies sollte 
aber kein Grund sein, die Begriffe zu verwechseln, wie es gelegentlich geschehen ist. 

In Abb. 1, f wird ebenfalls der Ubergang von der rein elastischen Form der Biegelinie zur rein 
plastischen beim idealen I-Trager unstetig, d. h. die Zwischenformen entfallen. 


(Eingegangen am 11. Juni 1951.) 
Anschrift des Verfassers: Professor Dr. H. Craemer, Alexandria (Agypten), Bulkeley, Immeuble Binoor. 


5 Pina hierzu auch die Betrachtungen in einer friitheren Arbeit: H. Craemer, Ing.-Arch. 13 (1943), 
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Ausnutzungslinien als Darstellungsmittel in der Plastostatik. 


Von H, Craemer, 


1, Grundlagen und Aufgabe. Wahrend in der Elastostatik alle von einer Last hervorgerufenen 
Wirkungen — Spannungen, Stabkrifte und -momente usw. — dieser zufolge der Hookeschen 
Annahme unmittelbar verhiltnisgleich sind, trifft dies in der Plastostatik nicht zu; statt dessen 
haben wir es hier mit Kurven zu tun, die entweder krummlinig sind oder aus mehreren geraden 
Stiicken bestehen. Im Folgenden beschrinken wir uns auf Stabwerke unter Biegung, die der 
bekannten Annahme idealer Plastizitat 


Gee far oes Os oy, C= Of, 8 LULS 6p 6 2.00 (1) 


folgen, wobei o, und e, die Spannung und Dehnung bei Erreichen der Streckgrenze bedeuten; 
o, soll fiir Zug und Druck gleich sein. 

Die Gesetzmafigkeiten treten klarer hervor, wenn wir sie in dimensionslosen Groen aus- 
driicken. Wir definieren deshalb als Ausnutzung eines Elements 


Co 


b= (2) 
als Ausnutzung eines Querschnitts 
- (3) 
ane: 
und des Tragwerks als Ganzem: 
(0) === ia 3 


hierin ist M das in einem bestimmten Zustand wirkende Moment eines Querschnitts, P die vom 
System aufgenommene Belastung, M, und P, die entsprechenden Héchstwerte bei voller Her- 
anziehung der Plastizitat auf Grund von (1). Sind mehrere Lasten vorhanden, so wird an- 
genommen, daf sie alle im gleichen Verhaltnis zanehmen. Im folgenden sollen die Beziehungen 
zwischen den oben definierten drei Gré®en fiir einige charakteristische Falle untersucht und als 
Ausnutzungslinien dargestellt werden. 


2, Ausnutzung von Querschnitten, Wir betrachten einen I-Querschnitt nach Abb.1 unter 
Vernachlassigung der Spannungen im Steg. Solange die Randdehnungen ¢,<e, sind, d.h. 
im Hookeschen Bereich (Abb. 1a) ist 


os Vi 1 obM, (5) 
nach Ejinsetzen der Plastizierung (Abb.1b)! 
M= Caras 


Fiir die Elemente des elastischen Kerns ist nun o/o,= &/f, also 


*o% ee) bh? 
= es So 4 bl 6 
oe (3 o : Bom fale (6) 
und hieraus folgt im Falle o =o,, =a, d.h. bei voller Plastizierung 
2 
M, = (1—a) 0," (7) 


Dividiert man nun (5) und (6) je durch (7) und zieht die Begriffe nach (2) und (3) heran, so 
folgt im elastischen Zustand 


ale 08 
Sk Sa (8) 


1 H. Craemer, Ingenieur-Arch, 13 (1943), 5. 285. 
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bei teilweiser Plastizierung jedoch im elastischen Kern 


= 3 Be (9a) 
3 (1 —a?) 
sowie 
pean (9b) 


in den FlieBgebieten. 
Mit €=1, w=1 erhalt man aus (8) und (9) den gleichen Wert, namlich 
Belo) aaa 10 
rR | (10) 
dies ist diejenige Ausnutzung des Tragwerks, bei der gerade an einer Stelle die Streckgrenze 


erreicht wird; der Kehrwert 1/», wird teilweise als Biegefaktor bezeichnet und gibt an, um wie- 
viel das Biegemoment nach dem ersten FlieBen noch gesteigert werden kann. 


7 


Abb, 1. Spannungen im I-Querschnitt. Abb. 2. Ausnutzung der Elemente. 


Fiir eine bestimmte Querschnittsform, gekennzeichnet durch den Wert a, und ein bestimmtes 
Element, bezeichnet durch €, kann man aus Vorstehendem die Beziehungen zwischen 4 und y 
darstellen. Fiir a0, also das volle Rechteck, wird im elastischen Zustand 


2 
Pee (11) 
und fiir den elastischen Kern nach FlieBbeginn 
a 
ferner ist 
es =. (13) 


Abb.2 gibt den Verlauf der Funktionen (11) und (12) an; man erkennt das raschere Ansteigen 
der Ausnutzung der Elemente, nachdem der Rand begonnen hat zu flieBen. 

Der andere Grenzfall, a—1, betrifft den ,,idealen 1-Trager‘+, bei dem die ganze Masse in der 
auBersten Faser zusammengedrangt ist; er ist nicht weit entfernt von wirklichen I-Tragern, bei 
denen a meist etwa 0,9 betragt. Es kommen hier nur die Elemente £=1 in Frage; Gleichung (8) 
liefert hierzu zunachst einen unbestimmten Ausdruck, der nach den iiblichen Regeln zu 


9 (14) 
bestimmt werden kann. In ahnlicher Weise findet man aus (10) 
ye (15) 


Man sieht, da mit voller Ausnutzung der Randelemente hier auch die volle Ausnutzung des Quer- 
schnitts selbst gegeben ist, so daB (9) gegenstandslos wird. Es sollte mehr beachtet werden, daB 
dies nur fiir den idealen I-Trager gilt. Die Ausnutzungslinie ist hier eine gerade Linie (Abb. 3), 
d.h. alle Linien der Abb.2 sind in eine einzige zusammengedrangt. Es sei bemerkt, daB dieser 
Querschnitt statisch bestimmt ist, da es nur eine einzige Spannungsverteilung gibt, die den 
Gleichgewichtsbedingungen geniigt. 
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Fiir das gleichschenklige Dreieck mit um eine waagrechte Achse drehendem Moment beschran- 
ken wir uns auf einige wichtige Zustinde. Im elastischen Zustand (Abb.4a) ist die obere und 
untere Randspannung 


5 24M 12M 
(eo git eee (16a,b) 
in der Schwerlinie ist 
| o,—0 (16¢) 
und in den Abstinden y,=0,375 h und y,—0,293 h von der Grundlinie 
Moor M 
ae > PaO — 1,445. (16d,e) 


Das FlieBen setzt zunachst am oberen Rande ein und breitet sich von da nach unten aus; 
wenn im unteren Rande die Flie8spannung gerade erreicht ist, gilt Abb. 4b. Aus der Gleichheit 
der Zug- und Druckkraft findet man den Abstand der Nullinie zu ¥y—=3h/8 sowie aus dem 


Abb. 3. Ausnutzung des 
idealen I-Tragers. Abb. 4. Spannungsstufen im Dreiecksquerschnitt. 


Gleichgewicht gegen Drehen das wirkende Moment 
a 
M= 96 72 Oh? (17) 
die Spannungen sind 


N= O7, == O35 Oy = a O59 o,= 0, 0, = 0,22 6, (18) 


Bei voller Plastizierung schlieBlich (Abb. 4c) erhalt man auf gleiche Weise y, = f a = 12 | h 


= 0,293h; die im Schrifttum anzutreffende Angabe, dai hier die Nullinie durch den Schwer- 
punkt geht, ist voreilig, vielmehr ist sie bei vollplastischer reiner Biegung fiir jeden Quer- 
schnitt die Flachenhalbierende. Das Moment in diesem Zustand ist 


1 1 
M,=—(1— — 0, bh? = 00,0977 o, bh? . 19 
x ( ia) 4 


Driickt man wieder die Momente und Spannungen nach (2) und (3) aus, so wird im elastischen 


Zustand a 


My =— 2,347, M=1,1Ty, w=, o=— 0,159, w=0,14%, (20) 
weiter im Zustand } 
(0.103. fg, 1s gS OL fp = On fe 0,22 (1) 
und bei voller Plastizierung 
Po= bu=—Mao==1, MH, wl; (22) 
ferner ist 
v= oi = 0,428. (23) 


Die vorstehenden Angaben sind in Abb.5 dargestellt und, da die Zwischenzustande nicht be- 

handelt wurden, so gut es ging, erginzt. Die fiir das Element c schematisch angedeutete Ver- 

zweigung ergibt sich, je nachdem man einen Punkt unmittelbar unter oder iiber c betrachtet. 
10 


Ingenieur-Archiv 
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3. Ausnutzung von Tragwerken, Wir untersuchen einen voll eingespannten Balken mit 
konstantem Querschnitt und zwei Einzellasten P in den Drittelspunkten (Abb.6). Dann ist 


im elastischen Zustand 


ae 
M=—=2 pi, ) Meer Me Pie (24) 
a 9 9 m 9 9 
Der Stiitzenquerschnitt beginnt als erster zu flieBen. Wir legen im Folgenden ideale I-Trager 
zu Grunde, so daB M, bei FlieBbeginn sofort seine volle Ausnutzung besitzt und von da ab kon- 


stant bleibt; darnach ist 


Vie Mn= =z Pi— M,, Mo hee (25) 


Abb. 5. Ausnutzung verschiedener Fasern im Dreieck der Abb. 4. Abb. 6. Ausnutzung verschiedener Querschnitte. 


Dies dauert solange, bis auch zwischen den Lasten der FlieBzustand erreicht ist, das System somit 
voll ausgenutzt ist, d.h. 


Wie Wham Ie SAVE, ale (7— 3 (26) 
Fiihren wir nun die Ausnutzung nach (3) und (4) ein, so erhalten wir aus (24) 
4 x 4 
Vg = 2m = ZO 5 =| ian (27) 
Wir bezeichnen nun mit m, die Ausnutzung des Tragwerks, unter der an irgend einer Stelle zum 
ersten Male Fliefen eintritt und erhalten mit », = — 1 
Dimas (28) 


zu diesem Zeitpunkt ist demnach erst dreiviertel der theoretischen Tragfahigkeit ausgenutzt. 
Nach Bindung der FlieBgelenke an den Auflagern folgt aus (25) 
Yost i 10 


o—1 (29) 


und im Endzu stand 


a —— jel == (= -—— IL, (30) 
Abb .6 zeigt diese Beziehungen fiir die Schnitte a und m sowie x//=2/9 und x/l=1/6, d.h. fiir 
die Momentennullpunkte im rein elastischen und rein plastischen Zustand. Wahrend also in 
einem voll ausgenutzten Querschnitt (v=1) samtliche Elemente mit /4=1 ebenfalls voll aus- 
genutzt sind, trifft das Entsprechende fiir das Tragwerk im allgemeinen nicht zu; vielmehr 
erreichen bei voller Ausnutzung des letzteren (w=1) nur die Querschnitte am Auflager und 
zwischen den Lasten den Wert y—1. 

Das nachste Beispiel betrifft den gleichen Balken mit nur einer Einzellast, und zwar im 
linken Drittelspunkt. Im Hookeschen Bereich ist fiir die beiden Auflager a und b und den 
Schnitt c unter der Last 


M,—=— 0,1482 Pl, M,—=—0,0741 Pl, M,—0,0987 Pl. (31) 
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Nachdem der Schnitt a als erster flie&t, wirkt der Trager als einseitig bei b eingespannt mit 
konstantem Moment bei a und die elastostatische Untersuchung hierfiir liefert 


1 
M.=—M,, M,=>M,—0,1482 Pl, M,=— : M, + 0,1728 Pl. (32) 


Als nachster Querschnitt flieSt c und das System wirkt nun als Trager mit Einspannung bei 6 
und FlieBgelenk mit konstantem Moment bei c, ist also statisch bestimmt; es wird 


ANeaaMieeMe. M,=—+ Pl45M,, (33) 
bis schlieBlich auch Querschnitt 6 plastisch wird und 
2 eee ; PI (34) 


wird. 


Abb. 7. Ausnutzung verschiedener Schnitte. Abb. 8. EinfluB elastisch senkbarer Mittelstitze. 


Die Ausnutzungslinien fiir die verschiedenen Bereiche haben die Gleichungen 


4 2 8 

—— 3 Oo, 9 a 3 OM , i si 9 oO, (35) 

4 1 1 
y,=—1, y= — 3 (a) +e 9° ¥, = 1,555 w — Fon (36) 
== 7 se %=—6ow+5, (37) 
—) —) |, (ele (38) 

Die Giiltigkeitsgrenzen ergeben sich durch Gleichsetzen von (35) und (36) zu 
3 

0, = (39) 


(erstes FlieBen iitberhaupt, bei a) und aus (36) und (37) zu w,=0,965 (FlieBen bei c); die Aus- 
nutzungslinien sind in Abb. 7 dargestellt; die gestrichelten Linien dieser Figur werden spiater 
erldutert. 

Es gibt Falle, wo der Momentenausgleich bereits im elastischen Zustand vorhanden ist, z. B. 
der beiderseits eingespannte Balken mit Einzellast in der Mitte, wo M=-- P1/8. Hier sind die 
Ausnutzungslinien fiir die meist beanspruchten Schnitte Geraden wie in Abb. 3, wo hier @ statt v 
und y statt / zu setzen ware, und es wird w,—1, die Plastizitat kommt also gar nicht zur Geltung. 
Die Vorschriften, die bei Ermittlung der Spannungen nach der Elastostatik allgemein hohere 

Spannungen zulassen, um dem giinstigen Einfluf der Plastizitat Rechnung zu tragen, miissen 
-demnach als Notbehelf betrachtet werden; maSgebend ist hier vielmehr die GréBe «,. 

Abb. 8 behandelt den EinfluB einer elastisch senkbaren Mittelstiitze auf das Verhalten eines 
Lweifeldtragers mit der dargestellten Belastung. Im elastischen Zustand ist dabei 


M, = (—0,1785 + C) Pl, ~M,= (0,085 + 0,7) Pl, (40) 


wobei C eine die elastische Nachgiebigkeit der Stiitze kennzeichnende positive Konstante be- 

deutet; je nach ihrer GréBe wird der FlieBzustand entweder im Stiitz- oder im Lastquerschnitt 

zuerst erreicht. Der Endzustand ist von C unabhingig, da das System statisch bestimmt wird; 
10* 
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es ist 


M, Via Vin On 4eP ale (41) 
Die Ausnutzungslinien sind im elastischen Zustand durch 
y, = (—1,45+8,1C)m, v= (0,69 + 5,7 C) w (42) 


gegeben; ihr weiterer Verlauf ist leicht zeichnerisch zu finden, da die Endwerte bekannt sind. 
In Abb. 8 zeigt die volle Linie den Fall unnachgiebiger Auflager; bei C—0,03 kommt der Schnitt 
iiber der Stiitze spiter zum FlieBen als im vorgenannten Fall, bei C=0,0555 ist gerade der elasti- 
sche Ausgleich vorhanden, d.h. die Plastizitat wird gar nicht herangezogen und bei C=0,1 
schlieBlich flieBt der Querschnitt unter der Last zuerst. Mit wachsender elastischer Nachgiebig- 
keit der Stiitze drehen sich die fiir den elastischen Zustand giiltigen Teile der Ausnutzungslinien 
um den Ursprung. 

Neben den eben betrachteten, von der Belastung des Tragwerks abhangigen Auflagersen- 
kungen kénnen auch solche aus fremder Ursache eintreten; beieiner auf schweren Pfeilern ruhenden 
Briicke mit leichtem Uberbau kénnen z. B. durch das Eigengewicht der ersteren Senkungen ent- 
stehen, an denen der Uberbau ganz unbeteiligt ist, die in ihm aber, solange er elastisch bleibt, be- 
trachtliche Spannungen hervorrufen. Fiir den oben behandelten Lastfall ist dann 


Mi, =0.1785 PICK, M,= 0.0sser1- 4 Oke (43) 


fr 


=f} if 
Abb. 9. EinfluB erzwungener Senkungen. Abb. 10. Ausnutzung der Randfaser im Falle der Abb. 7. 


worin K eine den EinfluB einer einmaligen Setzung der Mittelstiitze ausdriickende GréBe be- 
deutet. Gleichung (41) bleibt bestehen und es wird im elastischen Zustand 


a 
Pil 


Vor Eintritt der Setzung ist natiirlich K=d6=0. Abb.9 gibt die Darstellung fiir 6— 0,02, 
6= 0,055 und d= 0,09, wobei angenommen ist, dafi die Setzung eintritt, nachdem das Tragwerk 
eine Ausnutzung w= 0,4 erreicht hat. Wiederum tritt mit VergréBerung der Setzung die Plasti- 
zierung iiber der Stiitze spater ein und schlieBlich flieBt, wenn 6 >0,055, der Lastquerschnitt als 
erster. Der Fall 6=0,055 entspricht wieder dem elastischen Ausgleich; er wird im Stahlbau 
manchmal durch kiinstliche Uberhéhung von Auflagern herangezogen, um einen Trager konstan- 
ten Querschnitts ohne Riickgriff auf die Plastizitat méglichst gleichmaBig auszunutzen. Die fiir 
den elastischen Zustand giiltigen Teile der Ausnutzungslinien verschieben sich hierbei parallel zu 
sich selbst. 

SchlieBlich soll noch das Verhalten des Tragwerks fiir den Fall erértert werden, daB die Stab- 
querschnitte allgemeine Form besitzen, so daB »,<1 ist. Im elastischen Zustand Andert sich 
hierdurch nichts, doch setzt das erste FlieSen jetzt frither ein. Bei dem Lastabfall Ab.7 z. B. wird 
die obere Giiltigkeitsgrenze der elastostatisch bedingten Gleichungen (31) und (35) von @, auf 
einen geringeren Wert w, herabgesetzt. Man erkennt leicht, daB allgemein 


Va —— 1,45 (60) - 8,1 ) 9 Ve = 0,69 ay -* Dall t) ° 0) = (44) 


/ 
Of) 


Ms | (45) 
ist. 

Im Gegensatz zum idealen I-Trager ist hiermit die Aufnahmefahigkeit des betreffenden Quer- 
schnitts nicht erschépft, sondern seine Ausnutzung nimmt, wenn auch jetzt langsamer, weiterhin 
zu; in den geringer beanspruchten Schnitten dagegen nimmt sie jetzt schneller zu als bisher. An 
die Stelle von (32) und (36) tritt demnach hier eine verwickelte Beziehung, die sich auf die geo- 
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metrische Vertraglichkeit der Forminderungen das allmahliche Eindringen des FlieBvorganges in 
den Querschnitt und seine Ausdehnung auf die Nachbarschnitte zu stiitzen hatte; fiir einen 
Sonderfall wurde dies in der oben genannten fritheren Arbeit des Verfassers behandelt. Das quali- 
tative Ergebnis ist gestrichelt in Abb.7 angegeben. Wesentlich ist hierbei, daB wegen der geome- 
trischen Vertraglichkeit volle Ausnutzung, y— 1, und damit unendliche Dehnung niemals in einem 
Querschnitt allein, sondern nur in mindestens drei Schnitten zugleich auftreten kann, im vor- 
liegenden Falle also in a, bund c. Volle Ausnutzung des Tragwerkes, w — 1, ist somit erst miglich, 
wenn die FlieBgelenke voll wirksam sind, d.h. wenn — », = —, =», = 1. 

Demgegeniiber wurde im Schrifttum teilweise die Ansicht vertreten, daB man zwar den Mo- 
mentenausgleich heranziehen diirfe, das Tragvermiégen des Querschnitts aber mit dem FlieSbeginn, 
vy = y,, als erschépft ansehen solle. Damit verwirrt man aber die Begriffe ; denn Abb.7 zeigt deut- 
lich, daB mit — y, = », der Ausgleich in den Querschnitten iiberhaupt erst beginnt und weit da- 
von entfernt ist, vollstaéndig zu sein. Auch das sog. Traglastverfahren ist in dieser Hinsicht logisch 
unbefriedigend, da es bei Untersuchung der Querschnitte eine Plastizitat tiberhaupt nicht zur 
Kenntnis nimmt; praktisch spielt dies zwar im Stahlbau, wegen y, ~ 1, keine groBe Rolle. Der 
Irrtum stammt daher, daB man den Momentenausgleich vielfach als nur den Gleichgewichts- 
bedingungen unterworfen angesehen und der neben diesen stets bestehenden Bedingung geo- 
metrischer Vertraglichkeit keine Beachtung geschenkt hat. 

Abb.10 stellt noch die Ausnutzung der Randfaser, ju,, in Abhangigkeit von @ fiir den Lastfall 
der Abb.7 dar; die vollen Linien beziehen sich wieder auf den idealen I-Trager und stimmen, da 
hierbei y = py, ist, mit denen der Abb.7 iiberein. Die gestrichelten Linien betreffen irgend einen 
symmetrischen Querschnitt. Fiir den Schnitt a wird uw,=1, also v= »,, wenn m= w,ist. Darnach 
bleibt u,, = 1, wahrend fiir die beiden anderen Schnitte, bund c, die Ausnutzung schneller zunimmt 
als bisher. Die Ausnutzungen, fiir die ,, = 1 und u,, = 1 wird, kann man ohne Schwierigkeit 
aus Abb.7 entnehmen, indem man dort — », =», bzw. v, = », setzt. 


(Eingegangen am 29, Oktober 1951.) 
Anschrift des Verfassers: Prof. Dr.-Ing. H. Craemer, Alexandria (Agypten), Bulkeley, Immeuble Binoor. 
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Berichtigung 
zu meiner Arbeit im Bd. XIX, S. 378 des Ingenieur-Archivs: 


»Tragwerksauslenkung unter bewegter Last. 


Von E. Pestel. 


Die Formeln der Arbeit sind nur giiltig, wenn fiir die ~,;(*) Formfunktionen eingesetzt 
werden, die von der Stellung der Last unabhangig sind. Sollen die Formeln in der dort ange- 
gebenen Allgemeinheit giiltig sein, so ist qn; (t;) Pni(*) durch q,(t) Pri x3 x;(t)) zu ersetzen und 
dann unter dem Summenzeichen in 2a) 


aa a 
CH Guo nF 


5 OB 
| a n 
r 2 dn Ox; 


Pn %; 
Xi 


sowie unter dem Summenzeichen in 2b) 


2 2 
| (2 qn se cao +2 5a ch +g me) 
hinzuzufigen. 
Dieselben Ergainzungen sind ferner in den zugehérigen Ausdriicken fiir die virtuellen Arbeiten 
der Tragheitskrafte und in (2.3) zu machen. 
Das in der Arbeit entwickelte Lisungsverfahren, das von (2.2) ausgeht, wird von der Be- 
richtigung nicht betroffen. 


(Eingegangen am 5. April 1952.) 
Anschrift des Verfassers: Dozent Dr.-Ing. E. Pestel, Hildesheim, Kaiser-FriedrichstraBe 15a. 
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‘Beobachten wu d Messen sind die grundlegenden Mittel zur Erreichung dieses Zieles. _Bauherr und oberste Bau- 
_ Teitung,- -denen die finanziellen Belange unterstehen, haben ein Interesse zu wissen, wie und was’ fiir meBtechnische 
% Mittel ihnen gestatten, Einblick in die tatsachliche Ausniitzung und Tragfahigkeit der Baustoffe und Bauteile zu 
geben. Der Biingeuicus, dem Projektierung und Ausfihrung unterstehen, ist mit der Theorie und Berechnung 
vertraut, Er wiinscht zu wissen, wie weit die Berechnungsgrundlagen, denen zablreiche Annahmen und Voraus- 

: - setzungen anhaften, mit dem tatsdchlichen “Verhalten” des Bauwerkes abereinstimmen. Auch fiir thn ist die MeS-— 
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. die Priifung der Bauwerke in Massenbeton. Von Dr. sc. techn. konsult. Ingenieur ay 
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Grundlagen der Aeromechanik und Flugmechanik. von p:-me 
Arthur Proll, Brotesort. R. ander Technischen Hochschule in Hannover. Mit 278 Text- 


abbildungen, XVI, 612 Seiten. 1951. sae Ganzleinen DM 48,— 


Ausden Besprechungen: Der Verfasser, von dem im Jahre 1919 ein Buch tiber Flugtechnik erschienen 
war, verdffentlicht auf dieser Grundlage ein ausfiibrliches Werk tiber Aeromechanik und Flugmechanik, Gefordert 
von dem Verlangen nach Verbesserungen in der Flugtechnik, hat die Aerodynamik in den letzten Jahrzehnten einen : 
bedeutenden Aufschwung genommen, Die schon.lange vorhandenen Erkenntnisse tiber ideale Flissigkeiten, Po- 
tentialstrémungen sind mit neueren Erkenntnissen tiber wirkliche Flissigkeiten zu einer fruchtbaren Synthege ver- 
bunden worden. Dieses Wissen hat besonders seit dem ersten Weltkrieg die Entwicklung der modernen Flugtechnik 
geférdert und diese zu groBer Vollkommenheit gefihrt. Der Verfasser war an dieser Entwicklung in hervorragendem 

" MaBe beteiligt. Sein Werk gibt als erstes nach dem zweiten Weltkrieg einen Uberblick tiber dieses ‘Wissensgebiet. .. 
Das Werk zeichnet sich in erster Linie durch die Fille des Gebotenen aus. Seien es aerodynamische oder flugtechnische 
Fragen, man wird, wenn sie nicht ausfiibrlich behandelt sind, zumindest einen Hinweis. zum eingehenden Studium 
finden. Dieser Vorteil in Verbindung mit zahlreichen anderen ergibt ein wirklich auegezeichnetes Buch, fiir das si 
dem Verfasser dankbar sein kann. Die Ausstattung des Buches ist sehr gut. 

», Maschinenbau und Warmewirtschaft* 


Die Welt der Vektoren. Einfiihrung in Theorie und Anwendung der Vektoren, 
Tensoren und Operatoren. Von Dr.-Ing., Dipl.-Ing. Franz Ollendorff, Professor der 
Elektrotechnik und Vorstand des Elektrotechnischen Laboratoriums, The Hebrew Tech- 
nical College, Haifa (Israel). Mit 68 Textabbildungen. VIII, 470 Seiten. 1950. 

v \ DM 37.50; Ganzleinen DM 40,— 


Aus den Besprech un gem: ... Der Verfasser fabrt seine Leser mit Leichtigkeit aus dem gewohnten 


euklidischen Raum in die abstrakten Mannigfaltigkeiten von Minkowski, Riemann und Hilpert, deren Metrik der 


Anschauung nicht mehr zuginglich ist. Die in eleganter Form gebrachten Anwendungen auf Weltvektoren und 


Materiewellen beweisen die Hinfachheit der Methodik, die die eingefiihrten Tensoren und Operatoren erméglichen. 


Der Verfasser wollte bestimmt kein Lehrbuch der Physik schreiben, aber er bietet durch seine duferst glickliche i 


Konzipierung gerade dem Physiker viel.” Alles in allem ein ausgezeichnetes Werk, fiir dessen Heratsgabe man Autor 


und Verlag nur danken kann. »Osterreichisches Ingenieur-Archiv‘ 


Lehrbuch der Funktionentheor 1@. Von Professor Dr. phil. H. Hornich, Graz. 
Mit 34 Textabbildungen. VII, 216 Seiten. 1950, DM 19.50; Ganzleinen DM 21.60 


e 
Das Buch gibt eine Einfihrung in die Funktionentheorie und will vor allem den Weg zum Verstandnis ihrer zahl- 


losen Anwendungen in Physik und Technik erschlieBen. An Kenntnissen werden nur die Anfangsgriinde der 


Differential- und Integralrechnung vorausgesetzt. Zuerst werden die auf einem Gebiet der Ebene eindeutigen 


Funktionen uhtersucht und erst mit der analytischen Fortsetzung auch die mehrdeutigen Funktionen behandelt, 


Jedem Abschnitt sind Beispiele zur Ubung und Erganzung des Stoffes beigegeben. 
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